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A PRACTICAL METHOD FOR CALCULATING 
THE COUPLING OPERATOR IN LINEAR APPROXIMATION 


L. N. Lipatov and Yu. P, Leonov 


( Moscow) 

Translated from Avtomatika i Telemekhanika, Vol, 22, No. 5, 
pp. 561-570, May, 1961 

Original article submitted August 1, 1960 


The problem is considered of obtaining a simple numerical method for determining the weight function 
of a linear object. As an example, the weight function for the coupling between the parameters of a 
continuously operating coil boiler is calculated, 


The problem of the statistical analysis of the relation between random functions has recently become of great 
importance, This is due to the fact that at a certain level in the precision of describing a complex system, it ceases 
to be clear how to write the differential equations of the system, In the operation of such a system, too many differ- 
ent factors affect the operation of the system (sometimes factors that are unknown to us), and it becomes difficult to 
take them all into account. This situation arises, in particular, in the solution of certain problems of automatic con- 
trol, when precision control of a certain object requires that it be described accurately. A statistical description of 
the coupling in this case can lead to the design of a more accurate control system. 


We will consider a simple method of obtaining a linear model of a coupling, which will be suitable for practi- 
cal calculations, We will use as an example the determination of the coupling between the parameters of a continu- 
ously operating coil boiler. 


1, The Problem 





Let there be an object relative to which two signals X,(t) and Y(t) are recorded by some device. The signal 
X;(t) can be assumed to be the input and the signal Y(t) the output of the system. It is desired to determine the coup- 
ling operator between X,(t) and Y(t) under the assumption that the problem is linear, In other words, we wish to con- 
struct a mathematical model of the object. 


2, Solution of the Problem 





To solve this problem, we can use the general method of finding the coupling operator described in[1]. In 
order to do this, we must seek the operator yielding the best approximation in the sense of that operator that gives 
the minimum of the mean square difference between the function Y(t) and the transform at X(t) of X(T), Le., the 
minimum of 

et = M {Y (t)— AX (x))*. 
(1) 

In agreement with the conditions of the problem, we consider operators At from the class of linear operators, 
Then Af X(r) can be written in the form 


AL X(t) = {we t) X,(t)dt 
As ” » Whee . (2) 


—oo 


From (1) we therefore have 
— 2 
ef = M IY (t)— \ w (t,t) X, (x) dx} : 


00 


(3) 


The usual procedure for minimizing this squared functional leads to the integral function 
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00 
\ w (t, t) Ux,x,(t, u)dt =Tyx,(t,t1) (—oe<ti<+o). (4) 


o 
—oo 


for the weight function, where x, x,(t;,%) andl’ yx,(t;,t) are the first moment of the function X,(t) and the mixed 
first moment of the functions Y(t) and X(t) respectively. In particular, if X,(t) and Y(t) are stationary random func- 
tions and have a stationary relation between them, and have zero mathematical expectation, then 


bi 
\ w(t)ky » (t— t)dt =k, (t) (-x<t<+o). (5) 


where kx, x; ( T) and Kyx, (7) are respectively the correlation of x,t) and the mutual correlation of the functions 
Y(t) and X;(t). 


The weight function corresponds in this case to a linear operator with constant parameters, The solution of the 
equation (5) is trivial if we calculate the two-sided Laplace transform of the left and the right hand sides. However, 
if it is required that w(t) satisfy the condition w(t) = 0 for t= 0, then we obtain a Wiener-Hopf equation, the solution 
of which is known from the theory of optimal systems [1]. If it is desired that a model of the object be realized by a 
physical device, then the condition w(t)= 0 for t = 0 must of necessity be satisfied. 


The problem of constructing a linear model of the coupling thus reduces to the solution of the integral equation 
(4) for the weight function of the linear operator, or, in the particular case, to the solution of (5). The first moments 
I x4 Xq(ty, t) and I yx,(t,t) are determined from a record of the random functions X<t) and Y(t) obtained in practice, 


Two cases can occur in the construction of 2 linear model of the coupling. In the first case, the linear model 
is approximate. This means that the operator A‘ ¢ the first approximation in (1) is nonlinear, but the best approxi- 
mation is sought in the class of linear operators. Ihe second case arises when the operator At’, for which the lower 
bound of ¢* is attained, is linear, The linear model of the object is then exact. 


As has already been shown, the equation (5) can be solved in closed form, at least for certain classes of corre- 
lation functions, In what follows, however, we will use the approximate solution. This is due to two circumstances, 
Firstly, the functions k,,,(t) and Ky,(t) are usually approximated by certain analytic expressions, and only after 
this approximation can the equation (5) be solved, The solution is then obtained by approximate methods, In addi- 
tion to the approximate solution of the equation, the question arises of the choice of the functions to use in approxi- 
mating the above functions, Secondly, the solution obtained by this procedure sometimes does not correspond to the 
problem. It is well known, in particular, that the solution of the Wiener-Hopf equation for the weight function to be 
used in a physically realizable system can contain delta-functions and their derivatives [2]. But a real object can 
not involve such functions. This difficulty is evidently connected with the fact that the functions ky x69 and kyx; (0), 
for which the equation is solved, must have Laplace transforms that are rational fractions, while the. teh correlation 
and mutual correlations functions do not necessarily possess this property in the general case. 


We finally mention the numerical solution of the equation (5) that does not require any approximation to the 
functions kx x0 and fo. In this case the integral equation (5) is reduced to a system of linear algebraic equa- 
tions. The number of these equations is usually rather large (equal to the number of values of the weighting function), 
and so a digital computer is required for the calculation. 


We consider below a simple graphical-analysis method for obtaining an approximate solution for the equation 
(5), suitable for practical calculations. It does not require the use of a digital computer. For the calculation, we 
must have a sufficient number of characteristics (moments) of the random processes at the output and input of the 
object, These moments are obtained either directly from the random processes themselves [3], or from the corres- 
ponding correlation functions. The proposed graphical-analysis method for obtaining the solution is naturally not 
the only one, but in the authors’ opinion it is suitable for practical calculations. 


We introduce the normed moments ©, of the function w(t): 


an = 


| 


(6) 
where 








+00 
eae \ w (t) tat. (6a) 


—oco 


We also introduce the central normed moments 
ie 
dng = —— \ (t— Gy)" w(t) de. (1) 


—oo 


Between the normed moments and the central normed moments there exist the following obvious relations: 


19 _ 0, Ooo i Og on (a1)?. 
Agq = Ay — Saga, + 2 (a,)%. 
gq = Oy — 4030, +- Gay (,)? —3 (a,)*. 


(8) 


se «© ££ WTI Gita ts: eae Se ae ee eee. oe Ss 


We will denote the normed moments for the functions ky x69 and kyx (9 by Yno and Ban respectively. 


It can easily be shown [3] that between Opp for a function w(t) satisfying the equation (5), ¥no and Bp, the 
relations 


Bo 
adi Yo 
a = bs 77 Tr: 


(9) 
no -- Rae no 
hold, 
In order to obtain an approximate solution of (5), it is convenient to have the form of the transfer function 
w(p), which depends on certain arbitrary parameters, After this, we can choose the m arbitrary parameters so that 
the m equations in (9) are satisfied, The moments Bn» and Yho are also assumed to be known, We note that the mo- 


ments can often be calculated from the Laplace transforms of the corresponding functions. Thus, if w(p) is the two- 
sided Laplace transform of the function w(t), then from (6a) we easily obtain the relation 


ty = (— 1)*w (p)|,_g- (10) 


It was shown in [3] how to determine directly values for the moments B np and Tho from the realizations x{ t) 
and y(t) of the random functions X¢t) and Y(t). 


We can, in particular, take w(p) to be of the form 


il (TP +T,, Pp +1) Ti (Top + 1ye— |} 
w(p) = kP* mt 


ii (TP +T,p+1) TL (Tp +1) 


— =] j= =] 


(11) 





Then, using the formula (10), we can easily express the moments of the weight function w(t) in terms of the 
parameters of the transfer function w(p). For the first four moments we have 


ty =k, 
MOIS) Ep en xomygs svi neg m, (12) 
= Hla + Ti — >) Tee — Sy Te +t. 

i=1 j=1 p=! q=1 





Re m, m, m 
due = Fi (TA 20) + STF — 3 (Th — 2p) — ST? 

fom j=l p=1 «=1 
ea ™ 's . i ™ 'g m, 3 m, (12) 
239 = > (T— 37:71) + > T; one ps (T 2p — 37 pT op) Bre > T%. 

im j=1 p=l q=1 

deo - 3 (T4—47uT? + 27,3) + s T;- 
3 


Fi (Tp — 40 pT? t+ 2Tt)— ST 


p=1 


It is somewhat arduous to solve the system of equations (12), It is therefore expedient to use an expression for 
the transfer function w(p) that may possibly be simpler, Here we will consider the graphical analytic solution for a 








function w(p) having the forms 
sits Tsp +1 e7? 
w (P) k apy oTep Ft ° (13) 
and 
k . 1 
w(P) =k ator t Die FD (13a) 


depending on the five parameters k, T, £, Ts andr. 
The transfer functions in a large number of industrial cases can be written with sufficient accuracy in either 
the form (13) or the form (13a). In these cases, the relations between the moments and parameters are 
O, = 2ET--+ +473. 
digg = ((28)° — 68] T? + T3. 
Og = ((2E)* — 16§* + 2) 74+ 73. 


The plus sign for Ts corresponds to the expression (13a), and the minus sign to (13). 
For the solution of the system of equations (14) it is convenient to write these relations in the form 








= k, 
a, = 2ET +4473. 
_ G0 F TS _ M0 F753 . 
~ FT? sp, (€) ~ seers ® (5). (15) 





T, = V |{(26)? — 2] T? — ago | 


In the equations (15) the minus sign corresponds to the expression (13a). The functions ¥(&) and ¥E€) are 


given by 
_ (2)?—2 (2E)? — 6 
WO) = am: ) — rere (16) 





These functions are shown graphically in Fig. 1. For more accurate calculations, the values of »¢E) and ¥<) 
and of the ratio ¥,/#, can be obtained from the table (see the appendix). 


The equations (15) can be solved by the method of successive approximations (either by iteration or by Newton's 
method [4]). By using the graphs of ¥¢£) and ¥{£) (or the tables of their values), we can solve the third equation 
of the system (15) and determine the values of T and E by successive approximations, In order to facilitate the use of 
the graph in Fig. 1, we can use the supplementary graph in Fig. 2. In the latter, we show the regions, relative to the 


quantities R, and R,, corresponding to various ranges of E. Here 
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deo +T3 


e—-. 
Gg + TS 


From (17), the third equation of the system (15) can be written in the form 
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Fig, 1 


(17) 


(18) 
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Region of zeros ® ® 
in the right half- : 
plane & 
Fig. 2 


It follows from (18), that for different ranges of the variable 
— , the regions of variation of R, and R, are sectors, These regions 
are numbered in Fig. 2 (the numbers in the circles) for the follow- 


ing ranges of variation of €: 


E=1— co for 
—E—0.92—1 for 
— = 0.86 —0.92 for 
—=0.7—0.86 fo 
—=0.38—0.7 fo 
—E = 0—0.38 for 


Oly. Oly. Og. Ag > 0: 

Gy. As. Ag. By > 0: 

Ay. Ae. as > 0, A, _ U0: (19) 
ae > 0. As. A, < 0: 

Oo. Og. ig < 0: 

Oe. Og << 0. a, > 0. 


(The signs of the moments are given to facilitate finding the in- 


terval of €,) 


If the graph in Fig. 2 is used, then for each approximate solution the values of R, and R, determine a rough 
value of €. From this rough value, a more accurate value can be obtained from the graph in Fig. 1 or from the table. 
If we now use the conditions (19) or Fig. 2, we can see that the solution of the system (15) we have obtained is unique. 


3, The Coupling Between the Parameters of a Continuously Operating Coil Boiler 





The method we have described for determining coupling operators (the operators of the object) was applied to 
the calculation of the coupling between the parameters of a continuously operating 67-2SP 230/100 boiler. This 
boiler was chosen for the object, since its characteristics are difficult to determine accurately by the usual methods 
(for example, by artificial perturbation methods), The nominal output of the experimental boiler was 220 t/hour of 
steam at a pressure of 100 atms. In its operation under a normal regime with a turbine, it was therefore highly un- 
desirable to introduce any perturbations into the system, Such perturbations could upset the whole normal regime of 


the boiler-turbine system. 


The nominal flow diagram of the boiler is shown in Fig. 3. One of the main elements of the boiler is a system 
of tubes into which water is directly injected. The tubes are heated by burning gases, and the water is heated andcon- 


verted into steam, In spite of the apparent simplicity of this continuous flow boiler, heat engineers have not yet 





obtained a satisfactory analytic relation between its parameters because of the complexity of the relation, which is 
described generally in terms of partial derivatives. In particular, no such quantitative relation has been obtained 
between the input parameters of the water and the output parameters of the steam, The aim of the present study was 
to find certain dynamic characteristics of the boiler. These characteristics were determined by two methods: 1) by 
the usual method of determining the frequency characteristics*(5); and 2) by the statistical method, 
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Fig. 3. Flow diagram for the 67-2SP 230/100 boiler: WEI) Water economizer, 
first stage; WEII) water economizer, second stage; LRS) low radiation section; 
HRS) high radiation section; TZ) transitional zone; RS) radiation superheater; 
CS) convection superheater. 


It is of great practical interest to obtain the coupling operators between the parameters of the steam and the 
water in a direct-flow boiler, This is the problem we consider below. In particular, a linear model was constructed 
of the coupling between the flow of water Q,,(t) (point A in Fig, 3) and the flow of steam Q{d (point B in Fig. 3). In 
order to do this, a realization of the random functions Q,(t) and Q(t) for the nominal regime (200 t/hour, p = 100 
atms) was recorded, and part of this record is shown in Fig. 4. It can be concluded from these realizations that the 
mathematical expectation is constant, and that the process can be assumed to be stationary with a sufficient degree 
of accuracy, 


Q;(9 
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Fig. 4 


In this case, the problem of obtaining a linear model of the coupling reduces to the determination of the 
weighting function w(t) from Eq. (5). For the approximate determination of w(p), the methods described in the previ- 
ous paragraph were used, It was assumed that w(p) had the form (13), and so the problem was reduced to that of find- 
ing the parameters of the function w(p). First of all, by using experimental results, the moments of the correlation 
and the mixed correlation functions were obtained from the realizations of q,(t) and q(t), For this calculation, the 


calculating device described in [3] was used, By an application of (9), the following values for the moments of the 
weighting function were obtained 
Gy = 1,08; a = 365; a9 = 98800; cao = 274-10; ayq = 8-10°. 
The parameters of the function w(p) were then obtained by successive approximation. 
The final values of the parameters were thus taken to be 


k = 1,08; 
t = 32; T = 205; — = 1,08; T; = 110. 


* The experimental frequency characteristics of the boiler were determined by V. Dement'ev (TsNIIKA). 
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For these values of the parameters, the frequency characteristic w(jw) of the object was constructed. It is shown 


as a continuous curve in Fig. 5. In this diagram we also show the frequency characteristic of the coupling (the dotted 
curve), determined experimentally by the usual methods of artificial perturbation, 


Im 





q O 
w-002 ‘ w=-002 rad/ sec 
rad/sec* \. 








w-0005 tad/ se=@ -Qn025tad/ sec 
Fig. 5 
A comparison of the two curves show that they agree in the region of low frequency, but diverge somewhat in 


the regions of intermediate and high frequency. This is to be expected, since the accuracy of the usual methods for 
determining frequency characteristics is appreciably lower for medium and high frequencies, 


Ky 6 Te Kygl T)s t/sec* 





Fig. 6 


Fig. 6 shows the correlation (the continuous line) and the mutual correlation (the dotted line) functions for Qw(t) 
and Q{t). In particular, Fig. 6 shows that the delay occurring in the system is 30 sec, and this agrees well with the 
value T = 32 sec that was obtained, 


The authors wish to thank Ya. I. Khurgin for several useful comments he made in reviewing the manuscript 
of the present article, 
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Table of Functions 






































g we 2 Vi/v2 z 7 he dibs 
0.00 +00 0 es) 0.74 | —0.1589) 0.610 —0.2604 
0.05 6.6555 | —0.1525 —43.643 0.75 —l.2222) 0.5806 —0.3827 
0.10 3.3108 | —0.3215 —10.298 0.77 | —0.3840} 0.5197 —0.7389 
0.15 2.1879 | —0.5314 —4.1172 || 0.80 | —0.7954 0.4150 —1.9166 
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1.6197 |—0.3199 
1.2727 |—1.2941 
1.0354 |—2.2970 
0.8594 |—6.2728 
0.8003 |-16.5735 
0.746 29.4227 
0.7203 | 12.5532 
0.6038 3.3756 
0.5000 , 2.0000 
0.4012 , 1.4311 
0.2991 1.1104 
0.1820 | 9.8945 
0.0961 | 0.7912 
0.0275 | .7286 
—0.0117 | 0.6984 
—0.0552 | 9.6688 
—0).1038 9.6395 
1.3891 | 1.9574 
1.2796 | 1.5617 
1.1905 | 1.3342 
1.1164 | 1.1854 
1.0538 | 1.0797 
1.0000 | 1.0000 
0.9122 | 0.8867 
0.8432 | 0.8086 
0.8139 | 0.7776 
0.7873 | 0.7504 
0.7409 | 0.7047 
0.7016 | 0.6674 
0.6677 | 0.6360 
0.6381 | 0.6090 
0.6246 | 0.5969 
0.6120 | 0.5855 
0.5886 | 0.5645 
0.5676 | 0.5457 
6.5486 | 0.5287 
0.5312 | 0.5131 
0.5231 0.5058 
0.5153 | 0.4988 
0.5005 | 0.4856 
0.4869 | 0.4732 
0.4742 | 0.4616 
0.4623 | 0.4509 
0.4567 | 0.4458 
0.4512 | 0.4408 
0.4408 | 0.4312 
0.4309 | 0.4221 
0.4246 | 0.4135 
0.4128 | 0.4053 
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0.3387 0.3354 
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The incomplete image of a weighting function is introduced, and its properties are given; it is 
shown that the concept is useful when determining the relation between the images of the corre- 
lation function of the signals at the input and the output of a linear dynamic system. 


Two-dimensional Laplace transforms have in recent years been applied to solve a number of problems, The 
existence of the image of the required function of two variables is a necessary condition for their use. 


A correlation function is a function of two variables and its Laplace image exists in most cases, It would, there 
fore, seem important to determine the relation between the images of correlation functions at the input and the out- 
put of a given dynamic system, 


For the sake of simplicity we assume a single input and a single output in our dynamic system. The correlation 
functions of the signals at the input and at the output are denoted by K,(t, t)), Ky (t, ty) respectively, When g(t, s) 
denotes the system's weighting function then [1] 


ty 


Ky (t. t)) = \et. 8) g(t. 8;) Kx (8. 8) ds ds. (1) 
0 


owe 


We multiply both sides of the formula (1) by e7 Pt at | and integrating with respect to t and t, from 0 to@ 
co co 


\ \K, (t. t,)e—?- a dt dt, = 
00 (2) 


co co it, 


= \ Jean | ( gt. 8) g (ty. 8) Kx (s. s,)ds ds, | dt dt. 


00 
The left-hand side of (2) represents the two-dimensional Laplace transform of the correlation function Ky{t,t) 
[2]. Denoting by Sy q)the image of the correlation function Kt, t)), and inverting the order of integration in 
the right-hand side of (2), we find that 


Dry (p. g) = ( \K. (s. s,)| g(t. s) edt g (ty. s,)e-* dt, |dsds,. (3) 
"See English translation, 0 0 8 % 
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We introduce the function 


G(p.s) = g(t. seat, (4) 


and shall call it the incomplete image of the weighting function, incomplete in view of its lower integration limit, 
The incomplete image of a weighting function is the outcome of an incomplete one-sided Laplace transformation of 
the original. 


One should’bear in mind that for physically realizable linear dynamic systems the incomplete image of its 
weighting function coincides with its one-sided image, as in this case we have g(t,s) = 0 when t < s, and therefore 


\ g(t. sje" dt = \ g (t. s) e—?' dt. 


vu 


By making use of the incomplete image of a weighting function together with its one-sided image it is possible 
to simplify the proofs of certain statements, as the region in which the function under the integral sign is not vanish- 
ing identically stands out more distinctly. 


The more important properties of the incomplete Laplace transform are similar in the general case (that is, 
when not only the weighting functions can be the originals) to the corresponding properties of the ordinary (complete) 
one-sided transforms [3,4] We cite the following: 


1, The limiting value theorem for the incomplete transform: ifs? 0, then 


lim G(p. s) = 0. (5) 
poo 


2, The theorem with regard to the product of an incomplete and a complete image: if G(p,s) is an incomplete 
image of the function g(t,s), and F(p) is the one-sided image of the function {(t), then the imcomplete image of the 
function 


t 
8 (t. 8) =\ et. s) /(¢—t)dt (6) 


is the product G(p,s) F(p). 


3. The theorem on the incomplete image of a product of functions: G(p,s) is the incomplete image of g(t,s) 
and ¢(t) a function which admits a Laurent expansion, then the incomplete image of the product g;(t,s) = ¢(t) g(t,s) 
is given by the formula 


G,(p. s) = ¢(—r) G(p. s) (7) 


wherer “5 ¢(-t) represents an operator acting on the function G(p,s). In order to find the value of the expression 


¢(-1)G(p,s) it is necessary to expand the operator ¢(-r) into a Laurent series; subsequently we operate with every 
term of this series on the function G(p,s). 


Using (4) we bring the relation (3) into the form 


90 60 (8) 
xy (p- 2) = \ | G(p. 8) GUg. 44) K,s, 51) dsds,. 
00 


The formula (8) shows that the image of the correlation function of a signal at the output of a dynamic sys- 
tem is given in terms of the signal's correlation function at its input and of the incomplete image of its weighting 
function, 
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We shall find the incomplete image of a weighting function in three cases which are of practical interest. It 
shall be assumed when determining the image that the initial values of the weighting function vanish together with 

the initial values of its relevant derivatives or integrals, Initial values of a weighting function or its derivatives and 
integrals are understood to be its values corresponding to the time moment t= s. 


1, The behavior of a dynamic system is described by the equation 
Q (p)z =/ (t) (9) 
where x(t) is the required response to the excitation f(t). Its weighting function is given by the equation 


Q (p) g(t. s) = b(t —s). (10) 


We multiply both sides of the equation (10) bye ~P* and integrate with respect to t from s to ©, In view of the 
rules governing the incomplete image of a weighting function, and those of the 5-function, we obtain 


e~Ps 
G(p. 8) = oy (11) 





where p is complex, It is obvious that om is the transfer function of our system, 
Substituting (11) into the formula (8) we get 
®,, (P» 9) 
@ a tL 
w(P ) = OOM (12) 


where jx (p,q) is the image of the signal correlation function at the input, Thus, in the case when the linear dyna- 
mic system is stationary, the image of the correlation function of the output signa! can be found by multiplying the 
image of the correlation function of the input signal first by the system transfer function with argument p, subsequently 
by the same transfer function but with argument q. The formula (12) has previously been obtained by V. L. Lebedev 
[5] and also independently by N. P. Sevost*yanov. 


The formula (12) serves as a basis for compiling the tables which simplify the determination of the correlation 
function of the output signal of a stationary dynamic system when its operation time is finite. The tables enable one 
to estimate the effect of separate elements on the accurate working of the whole system; they also enable one to de- 
tect the elements which matter most as well as those whose effect is slight, and thus reduce the analysis of the opera- 
tion of a given system to a more simple one, The tables, their derivation and also examples showing their application 
are the subject of a separate paper. 


2, Let the behavior of a dynamic system be described by the formula 


Q (p)ig (tz) = 7 () (13) 


where y(t) is a known function, x(t) is the required response to the excitation {t). The weighting function g(t,s) is 
given by the equation 


2 (P) [9 (t) g(t. s)] = 8(t—s). (14) 
We introduce the notation 


g: (t. s) = @ (t) g(t. 8). (15) 
Then, as in the formulas (10) and (11), the incomplete image of the function g¢t,s) is given by the formula 


e-P (16) 
Q(p)* 


G, (p. s) = 








But in accordance with the formulas (7) and (15), the incomplete images of the function g¢t,s) and g(t,s) are 


related by the formula 
1 
G(p.s)= en G, (p. 8) 


and therefore 








1 ep 7 i e~™ 
G(p. 8) = 9(—*) E ma G(q. &) = 9 (— 1) | Q wo) (17) 


where r = oe r= Sr . Substituting the expressions (17) into the formula (8) we obtain 





¥ { ©, (P- 9) (18) 
Div(P- 1) = Hee le (p)Q ca 


The square brackets in the formulas (17) and (18) are here to indicate that the operators $t=7) and Soh) ope- 
rate on the expressions within these brackets. . 


In the particular case of g(t)=t we have 


co co 
Duy (P- q) = \ \ ACTA CT benmaaan (19) 
P@ 


3, Let the behavior of the system be described by the equation 


Q; (p) {1 (t) 2] + Qo (p) bs (t) 2) = 7 (2) (20) 


where ¢ (t) and yt) are known functions, x(t) the required response to the excitation f(t). The weighting function 
is given by the equation 


Q: (P) [Px (¢) g (t. s)] + Qo (p) [Pa(t) g(t. s)] = 5(¢ —s). (21) 


We multiply both sides of the equation (21) by e~P* and integrate with respect to t from s to ©, We thus obtair 
the equation 


Q, (p) Pi (—r) G(p. 8) + Qa (P) P2(— 7) G(p. s) = e-”* (22) 


which determines the incomplete image G(p,s) of the weighting function. 


The equation (22) is essentially different from (21). In the equation (21) Q(p)and QXp) are operators and A(t. 
and ¢{ t) are known functions of the independent variable t, but in the equation (22) Q,(p) and Q{p) are known func- 
tions of a new independent variable p with ¢,(-r) and y-r) as operators, This signifies that for comparitively sim- 
ple ¢ Xt) and ¢t)functions, the equation (22) can be of a lower order than (21), and so the derivation of the equa- 
tion (22) seems useful. 


However, if even one of the functions ¢(t) or yt) is more complicated than ¢(p) or ¢{p), or if the weight- 
ing function has already been found, then using the equation (22) instead of the equation (21) is not justified. 


In particular, when ¢t) = t, yt) = 1, the incomplete image of the weighting function can be determined 
from a linear differential equation of the first order 


—Qr(p) SE” + Qe(p) G(p. 8) = e-™. (23) 


The solution of the above equation which satisfies the theorem on the limiting value of the incomplete image, 
can be expressed by the formula 

















du (24) 
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_— 


where the function T(u) is given by the relation 


‘ 2 25 
arr 4h = 





(a is fixed and can be selected in any suitable way). 


Similarly, 
oc 


i T (uw) 
G(q. 81) = ria) ies Qi (41) oy 





(26) 


Substituting the formulas (24) and (26) into the relation (8) and inverting the order of integration we obtain 


8 FE Tw) Tey (27) 
Dxy(P- q) = TOT) 0: (u) Qi (us) Oy, (u. uy) du du. 





Having compared the formula (27) with formulas( 19) and (12) we conclude that the. correlation function of the 
output signal of a dynamic system described by the equation 


Q,(p) (tz) + Qa(p)z = f(t) (28) 


is the same as the correlation function of an output signal of an open system consisting of three units in series: a unit 
with the transfer function T(p)/Q(p), a multiplier unit by 1/t, and a unit with transfer function 1/T(p). 


The above can be of use when determining the correlation function of the signal at the output of a nonstation- 
ary dynamic system. 


Example, White noise is fed to the input of a system described by the equation 
tx + (ct + 2)z = f(t). (29) 


We shall determine the correlation function of its output signal. 
Using the procedure worked out in [5], the equation (29) takes the form 


(p +¢) (tz) +2 =f (0). 


In the case under consideration we have 
QalipP=pte gith=t. Alp)=—1. galt)=1. 


In order to simplify, we take a= c-1, and therefore 
i 
T (u) = ute’ 


(30) 


i 
®,,. (u. u1) = 5 (t — ty) eth dt dts = ue : 


ou8 
ow 78 





O,,(P- q) = (p+c¢)(¢+¢) (u + c)® (uy + ¢)* (u + ) 


vo""8 


oo 
\ du du, 
q 


As far as the performance of the system is concerned, it is equivalent to one containing three units in series, 
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1 
1, A unit with its tansfer function equal to(,; jt. The image of the correlation function of the signa) at the 
output of such a unit is given by the expression 


®; (p. 9) = (p +c)? (¢q +)? (p+ q) ° 





Applying the inversion formula for two-dimensional Laplace transforms [2], we find thatthe correlation function 
Ky (t. t)) = a oO) 14 4  (t — ty) — eH 1 + ct 4 ety 4+ 2c%ty)) 
corresponds to the image @(p,q) for t= t;. 


2. A multiplier by 1/t. The correlation function of the signal at the output of the second unit is given by the 
expression 


Ky (t. ) = 





ae oN) 14 4 (tt — ty) — (1 + ct 4 cty-+ 2%). 


The function k{t,t,) has the image 





®; (p. 7) = (u + c)* (uy + c)* (u + uy) ~ 


ve 8 


co 
{ du du, 
q 


3. A unit with transfer function p+c. The formula (30) is now equivalent to 
® (p. q) = (¢ +¢) (Pp + ©) a (P- 9), 
therefore 


K, (t. 1) = (en +e) (Gr +e) K3(t. t)| 
By carrying out the operation we find for t = t; 


ett—t) 


K, (t. 4) = Gere 


(20%? — 2cty + 1 —e?), 


The value of a correlation function, when t = t, can be determined from its symmetry about the line t= t. 
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Block-diagrams of varying-parameter systems containing two kinds of generalized units are considered, 
viz, the lag and the phase-lead unit. Transformation rules are given of these block-diagrams obtained 


by the application of linear differential operators, The results are illustrated by an example, 


1, Introduction 





A block-diagram of a system represents usually, in a conventionally accepted manner, the interconnections of 
elements or units constituting the whole system, and being differentiated only by their dynamic properties, A system 
with variable parameters can also be represented by a block-diagram in the same way as the systems with constant 
parameters; the system will only contain four different types of units, called elementary units; the integrator, the dif- 
ferentiator, the adder, and the multiplier with a variable gain coefficient. In many cases, however, such a detailed 
breakdown may not prove appropriate, and it may prove sufficient to represent the system as consisting of bigger 
blocks. Reflecting the block-diagram structure there will be equations describing the motion of the system and es- 
tablishing the connections between the input and the output signals, Only a block-diagram, however, will clearly 
display the interconnections between the blocks as well as the propagation path of the signals in the system; this 


enables one to penetrate deeper and to get a true insight into the working either of an element or of the system as 
a whole, 


It is well known that the block-diagram method in automatic control systems with constant coefficients has 
been used for some time and that it has developed considerably (for example, see [1]). The difficulty in construct- 
ing a corresponding apparatus for varying-parameter systems lies primarily in the fact that it is not possible to alge- 
braize, in the general case, the differential equations which, of course, prevents one also from obtaining simple rules 
as was the case in systems with constant parameters. A way of tackling this problem was outlined in [2]. 


There are three main problems for which the transformation rules for block-diagrams representing varying- 
parameter systems can be used, 


First of all, by using the method and the rules of block-diagram transformations one is able to find, for prac- 
tically all feasible combinations, the sequence of operations leading from one system to the other one previously 
given. This indicates, that by taking into consideration the corresponding differential equations, an algorithm for 
transforming one system of equations into another one can be established. 


Further, insynthetizing linear optimem systems one has to construct a block- diagram of the optimum system 
such that it is actually realizable in spite of an unavoidable presence of certain elements and connections given be- 
forehand. This problem arises in the synthesis of compensating circuits. 


Finally, as an actual system represents a combination of units and elements, some of which may be complete 
assemblies, the necessity arises, when investigating the system, of transferring similar elements to other places leav 
ing their structure intact. 


The resulting block-diagram transformations are of course equivalent to the corresponding transformation re- 
-Sults on differential equations and so the result cannot be simple unless appropriate approximations are used. As the 
latter are only feasible when actually solving practical problems, we shall not be able to consider them. 


Transformations of such block-diagrams shall be considered below — those containing only two different types 
of generalized units: the lagging or the phase- leading one, 
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The general case of a linear system with variable parameters can be described by the equation 








s dz out d@zin 
2s ut —~— = 2 b(t) - 1 
root dt pao dt? (1) 


where xin is the input signal, x,,,, the output signal a(t) and b(t) are given functions of time t. 


We shall say that a block-diagram element is lagging if the relation 
between the input and the output signals is 





























in | Yo | 7 | Fout n i 
F =— d= out 
c > a(t) —= = m0 (2) 
iO dt 
Fig. 1 where x,,,,, is the output and x» the input signal of the element. 


These units will be denoted by > on the block-diagrams. 
An element of the block-diagram shall be said to be leading if the input and the output signals are related by 


a di 
y= >) o,(0—? (3) 


j 
t, dt 





where y is the output and x,_, is the input signal of the element. 
The notation adopted for these units on the block-diagrams is the symbol F. 


When y = Xo, the expressions (2) and (3) give jointly the equation(1), Therefore, a block-diagram of a varying- 
parameter system can be represented by the introduced units; the simplest way is to represent it by a leading unit 
whose output is connected in series with a lagging one (Fig. 1). 


The above introduced units can be found in block-diagrams representing automatic control] systems in various 
combinations and with differing parameters, in forward and feedback networks and are both of main and compensatory 
character, When opening the feedback network of a varying parameter system, we obtain a system of elements in 
series whose joint dynamic behavior depends on the point of opening, which is not true for a system with constant 
parameters, In other words, the units of a varying-parameter system are not commutative. This implies that with 
two feedback contours consisting of identical elements but connected in a different order, different transformation 


rules are necessary. 


2. Transforming Block-Diagrams with the Aid of Linear Operators 





Theequation(1) of a linear varying~parameter system can be written in its operational form as 
ViFou = V2% in (4) 


where 


n 


Vi =>dyai(t) pt. Vo =D} 5; (t) pi 
j=0 


i=0 


are linear differential operators and p = d/dt. 
In agreement with the previously introduced definitions, the equation of the lag element is written in the form 


LE our = Zo (5) 


and of the phase- leading element as 
y= Frip (6) 
where L and F are the respective linear differential operators. 











The transformation rules of block-diagrams considered below boil down in the end to finding some new opera- 
tors which are the results of transformations of the original ones, 
























Let there be two lag units with their equations 


La; = ZLo- Lee = TZ}. 


Required to find the operator representing their operation in tandem. Substituting the second into the first equa- 
tion, we have 


L, [LsXe] — Lo: 


The left-hand side of this equation should be regarded as operation L; acting on the operator L, but not as 
their product, That is why we introduce a special symbol 


Lo =LseL, 


to denote this operation. Thus the result shall be written in the form 
LoXe = [L, * Le} Le = Zo- 


An operator operating on another one constitutes the basic tool in block-diagram transformations, It is given 
by the formula 








Vi«V, - Sy 2y, 2 (7) 
a — ap’ , 
where 
VY. =Sa(p. Ve=Sb,()p 
= j=0 
= . : n+m 
Vo =V,eV,= >} Ag (t) p2 
“_ (8) 
where 
Aq(t) = > a; (t) Nig 
i=0 (9) 
@ aye, 
Ni, = Yi Ct pp? *** 8; (0). (10) 


j=0 
The differentiation symbol p in these formulas only relates to the functions which follow it directly. The for- 
mulas have of course all the properties resulting from the linearity of these operators, 


Without pausing to give simple examples of transformations which follow directly from the above formulas, we 
shall proceed to consider transformation rules for block-diagrams with feedbacks. 


Lead-lag block with a follow-up. From the block-diagram, Fig. 2a, we obtain 
Li, = Fe, e= To 4- Ze 
and hence 


[L— F | 2, — Fx. 
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Thus the obtained lead-lag block has its variable lagging part operator equal to L-F (Fig. 2, b). 
Lag block in network with feedback containing a lag block. From the block diagram Fig. 2, c, we have 
Lyg= 8. Let =%q. &@ = %y + 4 
and hence 
{[L,# L, — 1} = Leto. 
(12) 


Thus a feedback network of the above type can be represented as a lag- leading block with respective operators 
as in Fig. 2,.d, 
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Lag block in a network with feedback containing a leading block. From the block-diagram Fig. 2, e, we obtain 


Li,=t, 4 =Fr,, e=%+%; 
and hence 


[L — F] x, = 2p. (13) 


Consequently, the considered network can be made equivalent to a leading block with the operator L- F (Fig. 2,f), 


Lag-lead block in a network with feedback containing a leading block. From the block-diagram in Fig. 2,g, we 
obtain 


Lt = Fie. x, = F422. @€=%+ 7 


and hence 


[L, —_ F, ok F,) Le => F ,2p. (14) 


We thus obtain again a lag-lead block but with a different lag operator. 
Lag block in a network with feedback containing a lag-lead block. From the block diagram in Fig.2h we obtain 


[yt,=8. [yr,= Foyt. &€ = 2+ 2, 























[Le * L, <— F,| Ts = Leto. 


The studied network is therefore transformed into a lag-lead one with appropriate operators, 


The discussed block-diagram transformations illustrate adequately in what manner the linear differential opera- 
tors are used, Other relations could be obtained in a similar way. The resulting new operators are easily obtained by 
using the formula (7), One does easily notice, however, that in all our examples we have managed to eliminate the 
intermediate variables by performing appropriate substitutions, If the latter cannot be carried out, a transformation 
of a block-diagram by our worked out methods does not prove possible. For example, in the equations 


Lyx, dy. Let, a Te 
(16) 
one cannot eliminate the intermediate x, by our methods, 
In order to solve the problems in which more complicated block-diagrams occur, when it is required to elimi- 
nate such intermediate variables for which the equations are not solved, a method is used which shall be called the 


equalizing operators method, Operating on the left- and the right-hand sides of the expression (16) with some, still 
undetermined operators U, and U,, we have 


(U, eL,) 2, = U2. (U,*L,] x, = U2, 


and we demand that it be true identically 


U,#L, = U,#Lp. 


Then the left-hand sides of the expressions in(17) are equal to one another, and we obtain 
U ,%o = U 9%o. 


In this manner the intermediate variable has been eliminated, 


Now it is necessary to find the operators U, and U,, which would ensure the equality of the left-hand sides of 
the expression (17). We note, first of all, that for the operators to be equal it is necessary that their orders as well as 
the coefficients of the same arder be equal. Consequently, if in the expression (18) the operator L, is of the order n 
and the operator L, of the order m, the orders v andy of the operators U, and U, respectively should be such thatthe 
equality 

v+n=p+m=—r 


be satisfied where r is the common order of the expressions (18), As r= n+ m, therefore v= m and p = n, 


By comparing the coefficient of equal powers of-p in the expressions (18) of the operators U,* L, and U,* Ly, 
one obtains r+ 1 linear algebraic equations between the unknown coefficients of the operators U; and U;. We have 
shown that the operator U, must be of the order v = m and it must, therefore, contain m + 1 coefficients, and simi- 
larly that the operator U, of the order » = n will contain n+ 1 coefficients, The total number of unknowns is there- 
fore r+ 2, and one of the coefficients either of the operator U, of U, can be assigned an arbitrary value. 


Thus, let 


n 


L, = 5 api. Ly = \ bsp. 
j=0 


i=0 


U,;= > jp). U,= Bip' 


j=0 i=0 





where the coefficients a, and 6; are still unknown and to be determined. 
The operators U,* L, and U,* L, are expressed in the form 


U,«L,= >) An P*. U,*L, = >) Be p* (21) 
k=0 k=0 
By taking also the identity (18) into account, we obtain 


Ax = By (kan. 1. 2.065. r)e 


By applying the formulas (8) and (9) and from expression (19) and (20) we obtain 


m n 
| > Nig =) Bi Nia (= 0.4.2.....7) (22) 
Bi i=0 j=0 
r| where the functions 
if qa “BE 
Nig _ p> Chet tate in. Ne. al >» Ci atk i-atky (23) 
i k=0 k=-0 


are known because the coefficients aj and bj are known. 


By putting ©, = 1 in the equations (22) which is possible because the value of one unknown can be assigned 
arbitrarily, and by transposing the terms containing the unknowns to the left-hand side, a system of r+ 1 linear alge- 
braic equations is obtained 


aoNoq — BoVoq + aig — BN iq + + «+ + Om—1Vim—1q — BnVing = — Ning 
(q = 0. 4, 2..... r). (24) 


—————— 


t 
. 
} 


whose solution will provide the sought coefficients oj and 6j of the equalizing operators U, and Up». 
We shall now consider block-diagram transformations in which the above described method should be used. 
Moving a leading block behind a lag block. From the block diagram shown in Fig, 3, a, we obtain 


Lr = Tin . Tout = F 2. 


Operating with equalizing operators 


(U,#L,)%,=Uyx,,. (U.* Fy)x, = UxSut 
we obtain 


O28 oy = Uitiy (25) 
provided that 


U,#L, =U, #F,. (26) 


The transformed block-diagram is shown in Fig. 3.b. 


Lag-leading blocks with feedback containing another lag block, From the block-diagram Fig. 3,c, we obtain 
the equations 


L,t,2= Fyx9+ Fy, Lyx, = 22. 


By applying the equalizing operators we get 
(UO, # Ly) x = [Uy # Fy) 29 + (Uy * Fy). (U2*L5) 2, = Uare 











and in view of the identity 


U,*«F, =U,#L, (27) 
we finally obtain 


[U, #L, —U,_) xg = (U, * Fy) 2p. (28) 


The resulting block-diagram is shown in Fig. 3, d. 
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Lag- leading block in a network with feedback containing a lag-leading block, For the block-diagram shown 
in Fig. 3,e, the following equations are valid: 


LyXq = Fyto +- Fi. Lots = F325. 


By applying the equalizing operators and performing transformations similar to those in the preceding case, we 
obtain 


[U,#L,—U,*F;,| 2. = [U,#F | 2 (29) 
U,# F, =U, "Ly. (30) 


The resulting block-diagram takes the form shown in Fig. 3, f. 


The preceding examples lead to transformation rules for other types of diagrams. The derived formulas and 
transformations enable oneto “convolve" theelements of any given system so that it will be transformed to any 
previously given type of system; in particular, this can be done when a single equation connects the input and the 
output. In the case of a constant-parameter system, the respective operators representing the components are formally 
identical with complex Laplace transforms expressions. All our formulas and results of transformations are, of course, 
also valid in the constant-parameter case, but in addition the formula (7) then takes a particularly simple form 


V;" Vz = ViV2 = V2Vi- 
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When constructing automatic control systems with varying-parameter components it is useful sometimes to 
make the whole system into a filter with constant parameters. To achieve this, it is obviously necessary to find the 
corresponding matching complementary blocks together with their mode of connection in such a way that the re- 
quired result is obtained. Let, for example, the block-diagram of a servo be, as shown in Fig. 4,a with block F of a 
constant-parameter type, It is required to obtain additional elements in such a way that the complete new system 
becomes a constant-parameter filter. We join the lag block-- to a negative feedback containing a leading block F,4 
(Fig. 4,b). The operator of the new block within the dashed lines is equal to 


In = L-+- Ky 


Further, assuming that 


L= >) a (t) pt 
i=o 


then, by putting 


Pag > by (t) p* 


i=n 


we obtain 


Lg = 3} fas (t) + b; (0) pt. 


i=a 


When the condition 
a; (t) + b(t) = const 


is fulfilled, the operator Lg has constant coefficients as required, Such an approach is advantageous as one can ob- 
tain the values of the coefficients of the operators Lo by a suitable selection of the functions b(t) such that they be 
equal, for example, either to the mean values of the coefficients a,(t) or, say, to their maximum values. In this way 
the dynamic characteristics of the new block can be made to become close to those of the old one at certain fixed 
values of its parameters, It is obvious that with several blocks with variable parameters, the above method can be 
applied to each block independently; moreover, it is even possible to find one big feedback including all the blocks 
with variable parameters. In order to achieve this, we must perform a transformation in accordance with the rules 
previously described, thus reducing the system to one with one or two blocks, 


Example, Given a system described by the equations: 


ay (t) 8 + a(t) 25 = 24 
dz: 
Lg = by (t) — 4- by (t) x2 
Cg (t) = + cy (2) et + eg (1) 4 = Tin —2%3— %4 








dx d. 
d,; (t) oo + do (t) 2 oy, = 71- y= e4 (t) OME +e (t)2 


t out 


Required a transformation such that the connection between the input and the output signal be described by a 
single equation, 


By introducing the differential operators 
1 = cop? + aqp+eg Lgp=ap+a, Ig—=ap+d. Fe=bp+bh. Fs3=ep+ eo (31) 


we can construct the block-diagram shown in Fig. 5,a, The network within the dashed rectangle can easily be trans- 
formed with the aid of the equalizing operators. We have directly from the Fig. 5, a. 


Ly =@—2s. Lore = 2. 2g = Foro. 



























Introducing the equalizing operators we obtain 
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.. Ss. See —4 as 
Fy U yal, <= Ugek o. (33) 
90 ae ae ee ay 
rip e vif, |af 7 | tout This gives us the network shown in Fig.5.b. Transforming 
fe U, U,*U,*L, qT; its part within the dashed rectangle, we obtain 
b (. ee awe oe om oe oe oe oe , 
I, P = ({U; +} U geL,\eL,) T out 
J 
which is shown schematically in Fig.5,c, Then, we obtain 
z P r’ i Zout Use = (U1 + U gal, )el4) 7 out 
4 U; 1U,* Ft on = Fin —* Fin — ¥s% ow = &- 
Hence, finally, 
fy 
Fig . Uzi, = (UgeFs + [U, -f U qeL,\4Ls) zx out (34) 


Us, = Lexy, Lo UeeFst+ (Ur + Uselilels. (35) 
In order to find explicitly the operator expressions (35), it is necessary to determine, first of all, the operators 
U,; and U;. In view of (33) and (31), these operators must be of the form 
Uy=ap+a. Us=Bip+ By 
and the operators U,* L, and U,* F, are both of the order r= 2, 
From the (23) formulas, we obtain: 


when q= 0 
NS,=0. NE ab, 
Nv = Pao. Nf, = po: 
when q=1 
No = % Noy = b 
NY, = % + Pa NE = by + phi 
when q= 2 


NS, = 0 NE =0 


Nyve = %1 NE, = by. 


Applying the formulas (24) and putting @, = 1 gives: 


when q= 0 
Agoty — by By — pboBi = — Pao 
when q= 1 
4,2, — bsB, — (by + ph) Bi = — (a + pa) 
when q= 2 





Oat -+ OB, — bibs = — 41. 
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tet. b (aob — boas) {1 [apbq + bo (par +- 29) — bi pag) — boas (pbi+ bo)} 
By oo b ai by a) {a; [bipao + Ur (pb, + bo) — a,pb,} —_ gb; (pa, + a)} 
Bi = * . =i. 


Now applying the formulas (7) we obtain 


UF; = U2Fs + a Fp = Brerp* + (Biper + Bots + Bre1) Pp + Bot» + Bipeo 
Usely = Ugly + ie Lip = Breap? + (Bipes + Bors + Bucs) p* + 
+ (Biper + Boer + Bicg) P + Bipeo + Bory. 


We have further 
U1 + Ugely = aggp*® + agap* + Gp + A 


where 


Qg3 = Bice. Gen = Bipce + Boca + Bier. Gar = Bi per + Bors + Bico + O% 
4g = Bi pty + Bory + Mo. 


Applying again the formula (7) we obtain 


- = OE ee «0 
Leela = Lyla + 5% lap + So lat + Gp lak 


= by p* + bys p* -{- Deop* + bap + bao, 


where 


beg = M430; 
bas = 3agspd; + 7 god) + asd, 
baa = Basaply + 5 candy + ayid + aeapds + Says p*dy 
bar = Gandy + Gan py + Sagsp*dy + Aggy ++ ag pds + agop*dy + aggp"d; 
bag = Gate + G41 Pdy + M42 p*dy +- agsp*do. 
By substituting the expressions (37) and (39) into (35) we finally obtain 


Lg = basp* + basp® + (bea + Bier) p® + 
+ (bar + Biper + Boer + Bier) P + Deo + Boto + Bipeo. 


By solving these simultaneous equations we obtain the required coefficients of the operators U, and U,: 






















(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


The coefficients of the operators lg and U, can be computed with the aid of the expressions (36), (38) and (40). 


APPENDIX 


We shall now find a formula expressing the result of the operation of one operator on another 


Given operators 


n ; m 
Vi= >) op v= >} b; p’. 
i=0 j=0 








Consider the following differential equation 


n 


peV = ral DX a; P| , 


i=0 
By applying the Lei>nitz formula 
n 
“> a, p'|= > Soh (p* a,) p*& + — Yoryr V) p*— (42) 
i=0 i=0 A=0 A=0 
where 
<A 
r 
pV = >) (p*a,) p. 
=0 


Thus the operation on the operator V, with the operator V, is equivalent to successive differentiation up to the 
n-th order of the terms of the operator V,, and the multiplication by the coefficients a; 


n 


VieV2 => > a; (p'sV3) 
i=0 


as the elimination of the intermediate variable x, in the equations 


Via, = 2%. Vota=— 7% 
leads exactly to the relation 


[VieV 2] TT, = 
When applying the formula (42) to the operator V, we obtain 
i 
P'eVa= >) CP (p* Vs) p*. 


»4=0 


We multiply both sides of the last relation by at) and add 


n 
> a, (p'eV2) = VisV2 = S Sicha, (p* V2) p** 
i=0 i--0 A=0 


or, changing the summation index 


n ft 
VieV—, = >> > cia, p Vop'. (43) 
i=o l[=0 


Considering the sum 


we note that 
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and this gives us the formula (7): 


vale 3} a! a v, +. 


By substituting the value of the operator V, into (43), we obtain 


== VieV, = y > 3} cia, pi b; piti 


i=o l=0 j=0 


1 
The order of this operator is obviously determined by the exponent of p 4" and it is equal ton+m, This im- 
plics that the operator Vp contains n+ m+ 1 powers of p with corresponding coefficients, and it can be represented 
in the form 


n+m 


- >> A, (t) p*. 
q=-0 


The expressions for these coefficient can easily be written if it is taken into account that the power with the 
exponent q is obtained by varying the index? from zero to q and the index j from q to 0 in such a way thatl + j 
always equals q. Then 


n it+j m 


nn > > >) a; ct-ati pi-a+i b; p? 


i=0 q =0 j=0 


and 


A, (t) = y y a, chat pitti p.. 
i=0 j=0 


We introduce the notation 


q 
mn >> ci-ati pati b;. 
j=0 


Nia (44) 


A, (t) = a,N,. 
a 2 Mi (45) 


The expressions (44) and (45) determine fully the formulas (8) and (9). 
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The analysis is carried out of certain nonlinear and extremal pulse systems in the incremental phase 
plane, The incremental-phase-plane method is made more precise and further developed, The boun- 
daries of the switching regions are determined, and systems with delay and systems with nonlinear and 
pulse correction are investigated. 


Nonlinear pulse systems, and in particular, relay pulse systems, were investigated in [1-6]. 


For the investigation of nonlinear pulse systems with amplitude modulation, J. Aseltine in [6] proposed the ap- 
plication of the incremental phase plane, in which the coordinates were the variable itself and its first difference, By 
the application of this method, which is, in our opinion, the most illuminating and convenient of those at present 
available, we carry out below the analysis of certain types of extremal and relay pulse systems, We will also make 
more precise and develop to a certain extent some applications of this method, namely: the determination of the 
boundaries of the switching regions, the investigation of systems containing loops with a delay and circuits with non- 
linear and pulse corrections, and the investigation of pulse systems of the third order, 


1, The Investigation of the Simplest Extremal Pulse Systems 





The principle of operation of extremal systems, including step systems, is well known (see, for example, [7]), 
and so we turn directly to the consideration of certain very simple extremal pulse systems, We will analyse the tran- 
sitional processes in these systems for impulse perturbations (or for nonlinear initial conditions), assuming as in [8,9] 
that, as a result of the perturbation impulse, the extremal characteristic is displaced and occupies some new quasi- 
motionless position, which for simplicity of calculation will be taken to be the origin for the analysis, The structural 
scheme of the system is shown in Figure 1, a, 


We write the equations of a section of the system as 


UL = Oyu (1) 

Ta - © = Aol (2) 
Q = — Or". 

(3) 


With the assumption that the coordinate x can be measured directly, or in the opposite case with the assumption 
that a model is available of the linear part of the object of control, the equation for extremal control of the relay 
type* can be written in the form 


sign A@n—, Atn—1 for n<t<n+y 
e for nty<t<n+1 (4) 


- ¢t 
where t =, T= const, y < 1, 


—_— 


*In [8] a system with proportional control was studied, 
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When we solve (3) and (4) simultaneously, we obtain 










u, = — sign (Z_ + Ln—)- (5) 





Further, the simultaneous solution of (1), (2), and (5) yields the equation for the dynamics of the system in the 
form 
A: 


Atay +(1—e° *) Ary = — a(1— ©) wags (6) 


where 


. . ’ T1 , 
Uns, = — Sign (Ln41 + In) = O, Ay Os t=+> a= a7. 























Fig. 1. a, Block diagram of an extremal pulse system with direct coupling for the coordinate 

x (or with a model of the linear part of the object of control); b, block diagram of an extremal 
pulse system without a model of the linear part of the object of control; 1) executive drive; 

2) linear part of the object of control; 3) link with the extremal characteristic; 4) memory 
unit; 5) unit for calculation according to the formula (4); 6) relay element; 7) tigger unit, 
carrying out the calculation for the formula (13). 
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z= 2,—-¥—* +—_* [—tu|izict—y +i |izioty, 


We now introduce the new variable y= 4x. Then, as was shown in [6], the equation of the phase trajectories 
corresponding to (6) is } 
~ mY : - (7) 
i—e * (, —e * ) 
where xX) and yp are the initial conditions, 

If we start from only the difference equations of pulse systems, which determine the coordinates of the system 
only for discrete instants of time, then the equations of the phase trajectories, obtained by J. A. Aseltine in [6], gen- 
erally speaking also only determine discrete points in the phase plane. We denote by xp_¢ and yp.¢ the coordinates 
of the system displaced by an arbitrarytime 6(0<¢=1) relative to the discrete instants n. The value of x,__ can be 
determined for a given transfer function of the continuous part of the system, by using the D- transformation [10] or 
the modified Z-transformation [11]. When we substitute xp.¢ and yp-¢ in (7), it is easily seen that for any value of 
€, all the values of xp.¢ and y,__ will also lie on the trajectory defined by (7). Thus the equation obtained in [6] 
determines at every instant of time the position of the representative point in the plane. 








For constructing the phase trajectories of the system, it is sufficient to have the equation (7), but these calcu- 
lations are more easily performed, and the properties of the system become clearer, if we consider the boundaries 
(or at least one of the boundaries) in the phase plane. With the aim of simplifying the further exposition, we intro- 
duce a definition. We will call the expression after “sign” the switching function F. 


The equation of one boundary of the switching region can be obtained by equating the switching function to 
zero, since it is evident that switching can not occur before F passes through zero. Consequently, one boundary, in 
the system being considered, is given by the equation 


Int 1 Fn = 0. 





























But since x,,,,= Yp + Xp» the final equation of the boundary MN of the switching region (see Fig. 2, a) is 
obtained as 


y = — 22. (8) 


Another boundary of the switching region can be obtained as follows, We choose values of xp and y» lying on 
the boundary F = 0 of the switching region using the phase- trajectory equations (7), and then calculate the value y, 
corresponding to x; = X9 + yo. The equation (7) is transcendental, and we solve it graphically, Starting from values 
of Xo and yo, we locate point by point a second boundary EOF of the switching region (see Fig. 2, a), the configuration 
and size of which is determined not only by the form of the switching function F, but also by the parameters of the 
linear part of the system.* 
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Fig. 2. a, The phase representation of the system (6) with the control (4) for 
T = 1, Xo = 4,2, yo = -1.35; b. phase representation of the system (6) with the 
control (9) fort = 1, a=31, A=1, xo=4, yo=0. 


If we know the character of the variation of the coordinate x, and use the equation (3), we can obtain the 
variation of the output ¢ of the object in control. 


In Fig. 2, b, we show phase trajectories for a system with the switching function 


FP = 2,4, T, — Asign y, (9) 


i.e. for a system in which the nonlinear element is a relay element with hystersis loop 24, The equation of one pair 
of boundaries of the switching region is obtained by setting (9) equal to zero; 


—2r—A for y>0 
—2r+A for yz 


As can be seen from Fig. 2, b, oscillations with period T, = 8T occur in the system. 


* There is no difficulty in principle in obtaining the switching region, In this connection, the method of investigating 
the problem in which the switching points located inside a certain region uniquely determined by the parameters of 
the system are considered as random variables [5], unnecessarily complicates the work, and is not adequate for the 
present work, 





We will analyze further the transitional processes in extremal pulse systems for the case that is most interesting 
in practice, namely when the coordinate x is not measurable in the control system and there is no model of the linear 
part of the system (see Fig, 1, b). 


Let the continuous part of the system be given by the 
equations (1)-(3) as before. The equation for extremal con- 
trol can be written in the form 





u, = sign (A@,_, -|- A) up_y 











(10) 
where 4 takes into account the effect of friction in the trigger, 








We remark that extrema! pulse controls, tested under the 
conditions occurring in industry, and rather well developed 
structurally (see, far example [12,13])actually have this type 
of control equation. * 

















If we solve the equations (3) and (10) simultaneously, we 
obtain 











u,, = sign (22 — x si 6 A) unt" (11) 











It can easily be seen that the dynamics of a system with 



































» control (10) are described by the same difference equation (6), 
0 \ 5<0 while the switching function for such a system is 
l oc 
F 4=— F 1Un (12) 


Fig..3, Phase diagram of the system (6) with the 

control (10) for Tr = 1,a= 1, A= 2.5, x9= 4, 

Yo= 0. By using the above described method, we obtain from (12) 
the equation x = (A - y*)/2y for the boundary of the switching 
region, 


This boundary divides this whole x, y plane into three regions (Fig. 3). In two of these regions, F; > 0, and so 
in these two regions, by (11) or (12), switching occurs. In the third region, for which F, < 0, a change in sign of the 
controlling quantity u is not possible, 


where F, = Xn 41° X47 4, 


The diagram in Fig. 3 shows that the control system we are considering is also perfectly efficient; its stationary 
operating regime is self-oscillating. The oscillation period of the system is determined by its parameters and the 
initial conditions. The greatest difference between these systems lies in the fact that when the representative point 
moves downwards in the region where F;> 0, more than one switching can take place in a system with the control 
described by (11), but only one can occur in the case of (5), Since the sign of the controlling parameter for the sys- 
tem (11) is determined by the product of two quantities, multiple switching can occur parameters fixed relative to 
one another, 


The choice of parameters to ensure that no multiple switching will take place (the presence of which will entall 
losses in hunting) can be made in the following way. We start with values of xp and yp that can occur in practical 
conditions and that lie on the boundary of the switching region (F;= 0), and use the method of successive approxima- 
tion to find parameters (a, T) for the linear part of the system, such that the system leaves the region F,; > 0 after 
a change of sign in u corresponding to a single switching. 


2. Nonlinear Pulse Systems Containing a Unit with Delay 





We will now use the difference phase plane to analyze the behavior of nonlinear pulse systems with delay. In 
order to simplify the description of the method of investigating such systems, we will use a concrete example (see 
Fig. 4, a). On the basis of what has been discussed above, we can write the equations 


1/r 
* It is not difficult to prove from the equations (1)-(3) and (10) that for T- © (corresponding to e - 0) and for 
4= 0, the systems with controls described by (4) and (10) become identical. 
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1 


A@r,, -| ( —e 7) Ar, = a (1 —e ) sign rw (73) 
x = Inco (14) 


n 


As can be seen from (14), an analysis of a pulse system with delay in the difference phase plane requires the 
knowledge of the values of x and y at the time t= n- €, or, to be more precise, it is sufficient to know only the 
value of x,-¢ and y;~¢, since all the subsequent values of xp-¢ and yp. ,_ can be determined from the phase tra- 
jectory, For the example under consideration, during the interval 0 = t= 1 the change in the coordinate x(t), which 
is the reaction of the system with transfer function W(q) to an impulsive disturbance of +4, is given by the expression 


alt 


r@=4alt—e ) + sesh daha a), (15) 


With t = 1, the equation (15) can be used either for finding yo 
= X,~ Xo for given x, and X9, or conversely for determining X» for 
given Xp and yp. When we substitute t=1-€ in(15), we obtain 


sok + U W/q) ify , fx)" the value of x;.¢, and from x;.¢ , by using the phase trajectory of 
a 



































the system, we determine the value y;.¢ (see Fig, 4, b). The further 
z’ investigation of a system with delay does not differ in any essential 
way from that already described above. 





The method we have described of using the time designation 
n-€ in the difference phase plane is evidently the only possible accu- 
rate method, since attempts to use the difference phase plane for 
known methods in the theory of continuous systems [15] for determin- 
ing the time on the phase trajectory do not ensure sufficient accuracy, 
and require complex construction. 


It is pertinent here to make a remark concerning the properties 
of nonlinear systems with delay, In contrast to nonlinear continuous 
systems, with straight-line switching, nonlinear pulse systems have a 
switching region, Thus the presence of a relative delay «€ (for Fig. 
4, b, € = €") displacing a representative point in the phase plane 
within the limits of the switching region does not influence the pe- 
riod or “amplitude” of the self-oscillations, This peculiarity of non- 
linear pulse systems with delay, as far as we can determine, has not 
been previously noted, 





Fig 4, a, Block diagram of a nonlinear pulse 
system with delay; b. phase plane of a non- 
linear pulse system with delay. Circles 
denote the x-coordinates and dots with 


crosses the x"-coordinates for € = €*, while The advantage of investigating nonlinear pulse systems in the 
dots denots the x*-coordinates for € = €"; difference phase plane in comparison with investigating them in the 
the broken curve is the phase diagram of usual phase plane is particularly striking in systems with delay, since 
the system for € = €", in these systems we must have recourse to rather cumbersome methods 


of designating the time in the phase plane, 


3. Nonlinear Pulse Systems with Impulse and Nonlinear Corrections 





Both continuous correction and pulse-correction methods can be used in the case of relay pulse systems, De- 
spite differences in the schematic realization of these two correction methods, their action reduces essentially to that 
of altering the locating (and configuration) of the switching region in the x, y plane so that the dynamic properties 


of the pulse system are improved. The influence of continuous correction on a relay pulse system in a concrete case 
is described in [6]. 


In many cases, especially for systems with high inertia, the practical realization of continuous correction often 
meets with considerable difficulties, because of the necessity of carrying out the differentiation of slowly varying 
functions, In such cases it is therefore preferable to use the impulse correction method, the practical realization of 
which can be shown to be simpler because only a definite number of memory elements are necessary. 
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We now consider the influence of pulse correction on the extremal system considered above.* 


Let the continuous linear part of the system have the transfer function 
, a 
OO) = Tea FH (18) 


and let the nonlinear characteristic of the object of control be described by the equation (3). The equation of exjre- 
mal control can be written 


u, = sign [Agn—, + BAGn—2 + A] Un. (17) 
The simultaneous solution of (3) and (17) yields 
Up, = — sign [x* — (1 —B)a2_, — Ba2_, — A] Un-s. (18) 


The phase trajectories of the system under consideration are obtained from (7) with the use of (18). The only 
difference between the system we are considering and the system with the control (10) is that their transfer functions 
are different. Thus, for example, for parameters the same as those in the caption for Fig. 3 and for 6 = -0.55, the 
phase diagram of the system with the control (18) is identical to that shown in Fig. 2, b. The introduction of pulse 
correction thus permits a decrease in the amplitude of the self- oscillations (the loss in hunting) by almost 50%, 


From the theory of continuous nonlinear systems, it is known that the introduction into a system of nonlinear 
correcting links permits the improvement of the quality of a contol system. We will show by an example that this 
remains true for pulse systems. We use the extremal pulse system with proportional action analyzed in [9] (with 
pulse-width modultion). 


The extremal-control equation for this system is 
sign u, = sign (Agn_,; + A) un, 
junp alot elAgeal if S+0|Agal<c 
\c = const. if 6+ 6|Aga_4|>c. (19a) 


(19) 


The limitation on the absolute value of the control parameter uy, is introduced to prevent any great over-regu- 
lation in the system. 


A rather simple method of realizing such a control was described in [13]. 
By taking (3) into account, we rewrite (19) in the form 


signu,, = — sign(z2—a?2_,—A)u,_, (20) 
6+ b|a—a2_.| if 8+6| 2-2 |<c 
ih if 6+ b|2% —22_|>c. (20a) 


The dynamics of the system we are now considering are described by the equations (6) and(20), We have used 
(7) and (20) to construct the phase diagram of the system in Fig. 5. For comparison, the same diagram shows the 
transition process for a system with the control (10). The difference in the character of the two systems clearly shows 
the advantage of the system with the nonlinear correcting element. 


We will consider briefly the properties of the system in question for 5 = 0, We note at once that the practical 
realization of such systems for 4 = 0 is not particularly difficult. A characteristic feature of the system for 6 = 0 
is that for a certain setting there is no self-oscillation in the system, This self-oscillation is unavoidable in extre- 
mal systems with constant steps. 


* The application of impulse correction methods in extremal pulse systems was proposed by O, M. Kryzhanovskii in 
{13]. One of the variants of the practical realization of such systems was also introduced in this work. 
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A further improvement in the dynamic properties of extremal pulse systems (increasing the rapidity of action 
of the system and decreasing the overcontrol) is obtained by the simultaneous application of both pulse and nonlinear 
corrections, To confirm this fact, we will analyze an extremal sy em with the control given by the equations (17) 


and (19a). 
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Fig. 5. Phase diagram of the extremal pulse 
system (6) with nonlinear correction, in the case 
when the control is described by the equations 
(20) and (20a); a, For r=1,A=1,a=1,6 =1, 
b= 0,05, c= 4, Xp = 20, yo = 0; b, the same, 
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Fig. 6, Phase diagram of the extremal 
pulse system (6) with nonlinear and 
pulse correction, in the case when the 
control is described by the equations 
(19) and (21a): a, For r=1, A=1, 





but with b= 0 and a@ = 1,58 a=1, 6 =1,c= 4,b= 0,05, 6 = -0,25, 
Xq = 20, yo = 0; b, the same, but with 
8B = -0.06 and c = (without any limit- 


ations on the absolute value). 


Just as for the systems considered above, the continuous part of such a system is described by the difference 
equation (6), where u is obtained from (18) and (20a), 


In Fig. 6 we show the phase diagrams of the system with the control described by the equations (18) and (20a), 
These diagrams were drawn from results obtained from (7) for two different adjustments of the system, It can be 
seen by comparing Fig. 5 with Fig. 6 that the simultaneous use of both pulse and nonlinear correction permis a sig- 
nificant improvement in the quality of operation of a pulse control system, The introduction of pulse correction 
with larger coefficients (see curve b in Fig. 6, for which 8= 0.06) makes it unnecessary to impose restrictions on the 
absolute value of uy. 


The example given above of the investigation of systems with pulse and nonlinear correction confirm the ef- 
fectiveness of introducing such correction methods in pulse extremal systems, and they also show that the study of 
such systems in the difference phase plane is not basically complex. 


4. Nonlinear Pulse Systems of the Third Order in the Difference Phase Plane 





We consider in conclusion the possibility of extending the application of the difference-phase-plane method 
of investigating nonlinear pulse systems to pulse systems of the third order, The analysis of the motion of pulse sys- 
tems of the third order, generally speaking, requires the consideration of the three-dimensional phase space of ax, 
4x and x. The investigation of such systems can, however, be reduced to a two-dimensional analysis. To do this, 
a method can be used that is similar to that used for the same purpose in the study of continuous systems, We con- 


a 


*It should be noted that it is not implied that there is no necessity of imposing limitations on the absolute value of 
Un in the presence of random perturbations in the control input. 
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sider the system shown in Fig. 4, a, Let its continuous part have the transfer function 





_— a 
YS) ~ Sapa eee 





(21) 
Then we can write 
X* (q) = 
__ te [(ba—br) (ee) (0-2) bs (eA) (00) +-by (€*—1) (€* 9] 7 a) (3% 
babs (ec — 1) (e%@—- e) (et — er) a 
The difference equation corresponding to (22) is 
A®z,, -}- A*x,C -- Az, D = uni, A + Uni, B (23) 
where 
C = 2 me e's penn et, D= eT ae ee: + e—rbs 
A = by— by + byes — byes. BB = e-(by — by) — bye + bye. 
We introduce new variables satisfying the equations 
z=A*z. v=Ar. w= A*zr + AzrC + 2D. (24) 
Then the. equation (23) can be replaced by the equivalent equations 
Atv, +- Avg€ + UnD = UnygA-+ Uay,Bo (eg = Av) (25a) 
Awn = UnzegA + Un4B. (25b) 


The investigation of the system is thus reduced to that of the difference equation of the second order (25a), the 
solution of which can be obtained by known methods, The variation of the coordinate w with time is easily obtained 
from a recurring formula obtained from (25b): 


Wnt = Wn UnzeA + Un4 iB. (26) 


SUMMARY 


It is shown in the present article that the difference-phase-plane method can be used in the analysis of a wider 
class of systems than was originally proposed in [6]. A further development of this method has made it possible to 
investigate transitional processes ir nonlinear pulse systems with delay in the simplest extremal systems, and also in 
systems using nonlinear and pulse correction, A method is indicated of determining the “switching region,” which is 
analagous to the switching lines for continuous relay systems, A generalization of this method for investigating non- 
linear pulse systems of the third order is described. 
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RELAY SYSTEM WITH SELF-OSCILLATING MODE OF OPERATION 
DISTURBED BY RANDOM SIGNALS 


V. Ya. Katkovnik and A, A. Pervozvanskii 


( Leningrad) 

Translated from Avtometika i Telemekhanika, Vol, 22, No, 5, 
pp. 599-604, May, 1961 

Original article submitted September 19, 1960, 


A basic problem in the theory of automatic systems is precise performance in the presence of random effects. 
Accuracy evaluation is based on the analysis of system motion under the effect of forces and random disturbances. 


The object of the present work is to investigate the dynamics of a self- oscillating system having relay-type 
nonlinearities, 


Fig. 1 shows the block diagram of such a system, 


The problem is solved by the method of adjustment [1]. The effect of small random forces is investigated for 
a system with an oscillating mode of operation. Discussing the problem of small fluctuations of a system, it is natu- 
ral to assume the undisturbed mode to be that periodicprocess which the method of adjustment describes qualitatively. 


The equations of motion of the system have the form 


r=—K(p)y, y=f(6,6), ¢=274¢ 


es aN 
(K@= ogy =a) 
(1) 
where €(t) is a stationary random process (with M{t} = 0 and M{t?} = oF; f(o, 0) is a relay characteristic with 
hysteresis loop and current brace (Fig. 2). 


f(4,d) 








= x7A 
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Let us introduce the new variables 


vy ee a 
p 20 p—h, (pom 4.2, . . om) = 
where A,(p = 1,2,...,M) are the roots of polynomial (Xp), which for ease of computation are assumed to be dis- 
similar (Re Xp < 0), Let us write equation (1) in the canonical form [1]: 


vA 
~ 





0 


Zo = Y- 29 — Ap Zp = y. <= >) BeZp. y = (6, 6). s=2-+C. (3) 
e=0 


Here, 


i hig etae Bi Pmt 
PC ae nS 


We examine the initially undisturbed mode (¢ = 0), Integration of system (2), at y = fp, yields 


Zo (t) = fot + Cy. 2p (t)= Cer! to | 
e 
Let us denote the coordinates of the system which describe a motion under the effect of signal fy, and a motion 
under the effect of signal -fy, through subscripts ; and 2, respectively, T, and T, are the corresponding time intervals, 


Using the continuity and periodicity conditions of the process in the system (the time count starts from the mo- 
ment of the last switch of the relay), we write 


Ze,1(0) =2.2(72), 2%.3(71) = 2.2 (0) (p=0.1..---2). (4) 
Then follows the systems of relations; 


Coa = —foT2 +Co,2 foT 1 + Coy = Co,» 


(5) 
Coa — ie = C,,,¢°°7* + ro 
) 


0 C eel feo __¢ , a fo 
" 6,1 i, e,2 i, . 
By adding the relay performance conditions 

%(0)=x+A. 2,(0)=x—A 


to relations (5) we get a system of equations with the unknowns T,, T2, Cp,1, Cp,2, Cp 1, C p,- From the first pair of 
conditions(5) follows 


T,=T,=T. 


The self-oscillation period is determined from the transcendental relation 





Bor. Sy Be det A (6) 
2 +25, th-- =— =. 
For the constants of integration we get 
1 1 , 
Co1= 3p, (2* — foBoT’) Co.2= 3g, (2% + foBoT’) (7) 
C,, ae s - —— ry C,, dine ae 


he 4 oo eet he 4 oo eet 
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The magnitude of period T, is independent of the x-shift in the relay characteristics, because % affects 
only the values of constants Co,; and Co,», It is necessary to verify the value of T, obtained from equation (6) with 
the aid of the condition for an adequate direction of switching 


n 


x (0) =- Bofo -i- >) Bp Ae Ce. 1 >> 0. 


el 


Let us now analyze the disturbed mode of operation. 


Denote through T,,(k = 1,2,...) the consecutive time cells ("half-periods") during each of which the relay 
output signal, y, is constant (takes the form fy or -fy). Let us assume for definiteness that the disturbance ¢ (t) is 
applied across the system beginning with the instant to (within the time cell T,) during which y = fy. Denote through 
$y the value for ¢ (t) at the end of interval T,,, i.e. , at the moment of switching. 


Let us further introduce the subscript k to denote the processes involved in the variation of all the variables 
over the interval T,. 


Then : . J 
29. k Cok |- (- i)*—1 fat. Zo.k = = C,. re rot - (— 1) ° 1.2 
A, 
(Time t is counted from the beginning of each interval.) 
The continuity conditions 
Zp n(T',) = Zo. k+1 (0) (p =0.4...., a} 


yield the recurrent relations 


Co, w+ (— 1) fo Tr = Co, eps, (8) 
C,. ,erotk -| (— 1)" fo = Ce. k+1 +- (— 1)*+1 #. (p =: 4. 2. — n). 
e 


Introduce the values which characterize a deviation from the undisturbed mode: 
AT, = (-- 1)* (7 —T,). AC,. x = Co oO, 


where T is determined from the period equation (6), while 


tin ye m4 eS te n), 
h, ite 
Ce at k = 2m —1 


Cok = Coo at k=2m. 


Furthermore, let us assume that the disturbances ¢,, are small compared to the amplitude values of the input 
signal,x(t), in the undistutbed mode of operation. Then, on account of the stability of the said mode, the values 
ATy, 4Cp i{p = 0,1... ), speaking generally, are also small compared to the corresponding values of T, Cp, 


Using relations (5), we linearize the discontinuity conditions (8) as follows; 
ACy, x = ACo. x14 ho ATy.. AC, y= e*et AC, k—1 “| C, ee? AT y_,. (9) 


Considering the formulas (9) to be first-order difference equations in relation to-quantities OC, 4{p = 0, 1.. + 
n), we readily obtain the expressions for the dependence of these quantities upon the fluctuations of the half- periods 
AT,,, as follows: 


k—1 k—-1 
ACy. x = fo bs AT,. AG. ¢= "ai ere kD AT. (10) 
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Let us now use the switching conditions 


ry (0) = (— 1) TA -- x—Ge 


(11) 
or 
Az, (0) = — Sx--- (12) 
At the same time 
Ax, (0) = 5 BpAzp, x (0) = > BAC, 4. (13) 


o=0 


Substitution of (10) into relation (13) yields the following infinite system of linear equations for the determina- 
tion of AT): 


k—1 k—} 
Bofo >} AT; +- s B,C, >) ee Te OAT, —- —t_, (k =2.3....). (14) 
i=1 e=0 t=] 


By subtracting the (k- 1) equation from each k-equation, the system is brought to the form: 


k 
> a ATi = — Oe + Ce (15) 


i=1 
where 


n 
ay = foBo-t- J BodgCge?*™ 


e=1 


n 
a = >) BAC, (e'e? —tpe OT © 2.3..... k). 
e=1 
It is obvious that the coefficients a; decrease rapdily with increasing subscript i. Their physical meaning is 

quite simple, The quantity a; constitutes the difference between the values at the moments iT and (i - 1)T of the 
response of the linear group of the system to pulse disturbance; however, not at zero initial conditions, where it 
would be the impulse function that corresponds to the operator — D > K(p), but for conditions which provide for the con- 
tinuity and periodicity of the processes in the undisturbed system. 


The solution of the system (15) is not too complicated because the matrix of the a-coefficient is triangular, 
and the quantities AT; can be expressed in sequence through all the ¢ j(i =k). For an arbitrary k we have 


k 
AT, = — >) Adn—its (16) 


ae 
where the coefficients A, are defined from the recurrent relation 
k~1 


A, = —+> Ajx—j44 (k > 3) (17) 


j=1 
it being evident that 


agel = . mem (29) 


— 


(the validity of relation (17) is readily proved from induction), 


Speaking generally, equation (16) fully determines the disturbed mode of operation. Substituting (16) into (10) 
and changing the order of summation, we get ae 


AC. = — fo > ti D Aj 
“ee as 
AC,,. r=—C >) O53) A er tt, T 


j=1 i=1 








The fluctuations of signal x(t) are determined from formula 


Arg (t) = BoACy, x + >) Boe’*'AC,, x 


p= 


or, if we allow for (12), from 
Az, (t) = Gy. — bs 8 (1 — e’e') AC, x. 
e=1 


The first item in the last expression characterizes the direct effect of disturbance in k-switching, while the 
remaining items characterize the effect of the fluctuations of the preceding half-periods, weakened due to system 
attenuation. Finally we write 


k—1 
Ary (t) = 5) Ce—iDj (0) 


j=1 (19) 
where 
1— A, }) B.C,.(1 — ee!) ere! numa J = 1 
D;(t) = Mee on 
? n 
<2 A; >) B.C, (1 — ee") * sities 2 nan jf = 2.3.... 
i=l p=1 


Let €(t) be a stationary random process with an assembly average equal to zero, and a correlation function 
Ry (t) = op, (1). 


It is obvious that the assembly average of the period fluctuations and the signal x(t) is also equal to zero, 
The dispersion of the period is given by formula 


oT, = M {ATi} = oe >) > AjApp; (i — J) T) (20) 


i=1 j=1 


and the dispersion of signal x(t) is given by formula 


k—1 k—-1 


: 21 
oz, (t) = of >) Di(t) Dj (tec —1) 7). (21) 
i=1 j=1 
If the random disturbance is a wide-band disturbance and 
Pe ((i—J)T)<1 upn i] 
the simpler relations 
k 
or, = 325) Af (22) 
i=1 
k—-1 
Oi, == & ye Di (23) 
i=} 


can be obtained, 
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In a similar manner it is possible to find the correlation coefficient between the half-periods 


koem™ 
M {AT AT m) = 0% 5) 3} AiAjpe (i — j -- m — ky}. (24) 
i 


The relation for the stable value Or is found from (20) or from (22), provided that the upper summation limit 
is extended to infinity . 


The fluctuations of signal x(t) are of a nonstationary nature, However, the change Ax(t) over the period T,, is 
fully defined by the fluctuations of the initial values AC, K p=0,1,...,), which constitute elements of random 
sequences, each of which has a stable character in the case of sufficiently large k. In these terms, it is possible to 
speak of the stable fluctuations, X}{t), as of fluctuations over a period that is sufficiently distant from the initial mo- 
ment at which the random disturbance ¢(t) was applied. For the dispersion of the stable fluctuations we have the 


following expression 


o-(t) = o& >) >} Dilt) Dj (t)e i — 7) 71. (25) 


i=1 j=1 
and at small correlation, the expression 


s(t) = at 3} Di(t). sae: 


i=1 


In spite of their apparent bulkiness, the given formulas are not too tedious in computation since the -oefficients 
Aj and D; diminish rapidly with increasing number as a direct consequence of mode stability. 


As an example, let us examine the following application of the technique described, Let operator K(p) have 
the following form: 


1 
(tip + 1) (tap + 1) 


and let the relay have no hysteresis loop (4 = 0), A half-period of the oscillations is defined by the expression 


K (p) = 





2 

BT btn ~ ER 

>= +> + Pike mene 
p=1 


At T,= 1 sec, % = 2 sec, \y=-1 sec”, A, =-0.5 sec”’, By = -1, 8, = -1, By = 2, we have T = 4.6 sec, 
The computation of a), may be limited to the three terms: 
a, = — 0,656 fo, az = — 0,326 fo, ag = — 0,0328 fo. 


Let us use the values obtained for a;, a2, a3 to calculate several values for the A,,~coefficients: 


Ar = 51.53. den — 2.330. dom 2.108. Age — 0.44. Am -1.0.153. 
0 


ho 0 ho to 
When the noise correJation function is approximated by the exponent 


R, (t) = oe*h"I 
from (13), we get 


Sp = 0,84 Fo aaa a = 0.01. Op = 1.67 0, nan a=0.1. 








In the case of white noise (Qo) = 5(T)), we have 
1 
Sr = 3,0-— «,. 
. fo * 


The wider the band of white Gaussian noise, the larger the spread in values of the oscillation parameters; at a 
very- low-frequency random effect, the standard deviation of the half-period diminishes. 


The discussed procedure, practically unchanged, may be used to investigate the disturbances in other periodic 
modes of relay systems. In principle, it may be extended to cover piecewise-linear systems [2]; it seems, however, 
that the impediment of the excessive computational labor involved can be circumvented in the practice only for 
systems of low order. 
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COMPENSATION OF DELAY IN A RELAY SYSTEM 


N. A. Korolev 


(Moscow) 
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pp. 695-612, May, 1961 
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A method of stabilizing a lagging automatic control system employing relays is explained, using a 
second-order system as an example, The method is based on compensation of the lagging by intro- 
duction of an anticipatory signal by a feedback to the system input from the output of the element 
that causes the delay. The transient processes in such a system are calculated and plotted. The con- 
trol function is computed, taking into account the derivatives of the quantity to be controlled, Com- 
pensation of the hysteresis loop in a static relay characteristic is examined, 


1, Fundamental Equation of the System 


Control systems with relays, as a rule, have a time delay, originating in the relay components, which notice- 
ably lowers the quality of the operation, increases the time during which transients are prevailing, and sometimes 
lead to instability of the control system. 
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R, W. Bass [1,2] has offered a method for compensation of the lagging in relay operated systems by the addition 
of a special feedback. In the present paper, an attempt is made to evaluate this method by means of an example of 
its application to an actual relay operated tracking system; the circuit diagram of this system is shown in Fig. 1, a, 
and a function block diagram of it is in Fig. 1,.b. 


The output shaft (ultimate power application) must reproduce the movements of the control or pilot shaft. The 


amplifier y provides a gain k,; the polarized relay Ip has a zone of insensitivity x, the repeater relay—a lag T. The 
transfer function of the motor is taken to be: 


K = tiles SS 
«(p) p(Mp +1) (1) 


and that of the repeater relay: 
Ky (p) = kpe-?*. (2) 


Hence, the transfer function of the linear part of the system is: 
k 


, . k 
K (p) = Ky(p) Ky (p) = sip rh i (3) 


The characteristic of the polarized relay is presented in Fig. 2 and is expressed mathematically as follows: 
1 at kyx > x. 


y = D(kyx) = 0 at % > kyr > —%. (4) 
—i at — % > kyz. 


Now we can write down the equation of the complete system in the open loop state: setting kiykp = k, this 
equation is; 


My (t) +  (t) = k® [ky x(t — T)}. (5) 
The equation for closing the loop is: 
Hr (t)— H(t) = 2 )- (6) 
In general ¥,{t) (the reference input setting) can be any function of time, but we shall consider here only the 
particular case when this function represerits a stepwise disturbance, that is, a sudden instantaneous change of the 
control angle setting by a constant value. For the sake of simplicity, we assume to have at t= 0; ¥(0) = %» and 


¥,= 0, Then substituting (6) into (5), and introducing © as the independent variable, we obtain the equation of the 
system in the closed loop condition: 


$ (0) + (0) = — Dx 1(9)}, n 
where 
i t T 
9) = pr Vi)» 9=ay tH 
—®, [9 )]=0.[— 9 0)1 = PI—9 0—)], 
1 at @> e. 
O(g)={ 0 at e>o>—«e. (8) 


—1 at —et>®. 

















When ky? = * we have g = €, Hence; 


kk *° (9) 


For any of the possible values of the function 6(¢), 
the solution of equation (7) will have the following form: 
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2,.Compensation of Lagging 





Using known techniques [3], it is possible to determine and plot the transient processes in the closed loop sys- 
tem described by equation(7). In Fig. 3 are shown two curves representing the transient process, one in a system 
with time lag (I), the other in a system without time lag (II). Both curves were computed for a system having the 
following parameters; M = 0,02 sec, k, = 30 ma/deg, ky = 4 deg/v*sec, kp = 25v, X= 3 ma, ~(0)= -1 deg. It is 
evident that lagging makes the effect of the transient process very much worse, 
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Fig, 3 Fig. 4 


The compensation method described in [1,2] is based on a simple idea; in a relay operated system, the reac- 
tion of the linear part due to the operation of the relay component is well known, Therefore, it is possible to create 
at the systems input, after time T has elapsed, a controlling command action 0(@ -T), which will completely com- 
pensate for the time lag, or, in other words, to introduce in the original command function (reference input) ¢(®) a 
forward shift of the time count reference point by T. We thus have: 


6(8—t)=@(6) wan 5(6) = (0+ 1) (11) 


where 9(6) is the command function of the system without time lag, Obviously, to make this possible it is necessary 
that the time interval between consecutive switching operations of the relay be greater thanT, 








The equation of the system with lag compensation is obtained from (7) by substituting for the control function 
(8) the function o@-T). We then have: 


@ (0) + 9 (0) = —®, (o (8)) 


(12) 
where 
1 at o>e 
(sc) = 0 at e>s>—e 
—i at —e>s. 


Let us find the control action o(8) by combining all the signals present in the system in a linear fashion. We 
then have; 


3 (8) = Be (0) + y@ (6) + aD [— 6 (0 — 1)}. (13) 


Taking into consideration (11) we rewrite (13) for the shifted time base as follows: 


6 (0 — t) = Bp (0 — t) + vp (8 — +) + aD [— o (0 — 2)] = = (8). (14) 


The functional scheme of such a system is shown in the block diagram of Fig. 4. 


Let us determine the coefficients a, 68 andy by solving equation (12) for the time interval between two con- 
secutive switchovers of the relay, Setting 


we obtain from (10): 


9(9)=G—1+60+e°. (15) 
Shifting our variable (argument) 6 byT , we have: 
g(6—t)=q@—1+6—t+ ete. (16) 
Substituting solutions (15) and (16) into(14), and making use of the identity 6(0-T ) = ~@) we obtain: 
B= 1. ¥=1i1-—e. a=t—i+e*. (17) 


Sometimes a system has a bad transients characteristic even after compensation of the time lag, It may be re- 
quired that the control signal not only be compensated with respect to the time lag, but also that some correction be 
introduced in it depending on the first and second derivative. For instance, let it be required that the control signal 
have the following form: 


A (8— t) = (8) + cy (8) + cx (8) 
where c; and c, are predetermined coefficients. 
Leaving the functional scheme of Fig. 4 unchanged, and taking into consideration that feedback produces a 
signal which is proportional to the second derivative of y, we can rewrite (14) as follows: 
A (8 — t) = Bp (8 — t) + 1p (9— 1) + a [— 2.(0 — 2t)] = 
= (8) + crp (8) + cop (8). 
The new values of coefficient: «, 8 and y are determined in a similar manner as before. We use the solutions 


for equation (12) within the time interval Ty between switchovers of the relay (provided Ty > 2r), that is (15) with 
the original variable and (16) with the shifted by T variable, and obtain from (19) the following equation: 


(18) 


(19) 











A (9— t) =B[@o—1 -} 09—t-+- ete) +- x4 —e'e*) 1a = 
(8) + exp (0) -i- cop (8) = Po — 1 4-0 + e 4 ey (1 — e-%) ego. (20) 


This equation must be satisfied for any value of @. Let us substitute in it for @ the variables 0 + ©, and @ = 0, 
and in addition equate the members free of variables on both sides. We then obtain: 


B = 1, 
Y= 1—e-*+4 (c, —e,)e-*, (21) 
& = Cee—* +- c, (1 — e-*) + e-* 4+ tr — 1. 


Substitution of these values for &, 6 andy into(20) turns it into an identity, 


When c; = cz = 0, the coefficients a and y assume the values given above in (17). 


3. Plotting of the Transient Process in a Compensated System 





For the sake of simplicity, let us plot the process g(@) at the condition c,; = c,= 0, We use equation (14), and 
the values for «, 6 and y derived from it and given in(17), To determine the function o(6-T), it is necessary to 
find its components 0, 02, and og (see Fig, 4), It can be deduced directly from this diagram that for 0< rT, there 
always is 0; = 06,= 63= 0, For6 2T, there will be: 


Gy =a 

Gg = 7 (1 — e*e—*) = 7 (1 —e-) (22) 
G3 = B (8) — 1 + e-) 

where 69= 6 -tT 20, Relationships (22) are valid until the effect of the next switchover of the relay comes through. 


Suppose that o reaches the threshold value for throwing over the relay o = - €9 at the moment of time @* . Hence, 
expressions (22) are valid when © varies within the limits of 0=6,9 =6},+T. 


In figure 5,a, b, and c the components 0,, 6, and 

be Gs are plotted, and in Fig, 5d, the sum obtained as a re- 
” sult of their addition (graphically or by tabulation and 
arithmetic). The point of intersection of the curve o 
with the horizontal straight line - €9 corresponds to the 
~025}- x moment of release of the relay 0}, Beginning with the 
: moment of time 63, +7, there will be other equations 
operative, derived in a manner similar to the derivation 
of (22), These are: 























4 =0 
0 Oo, = (_tu+* afl e~? 
, ’ (23) 
~Q25 Gs a Bo: ai (2% —1) e~*. 
~ 050} 
Q50 
Equations (23) will be valid within the limit of 
Q25+ fre: Obs + T 36, SO 55+ Tt where 801 is the moment at which 
? NF ot | 4fis : 78 — 0() crosses the straight line 6 = + &g the second thres- 
7h : fA hold value, corresponding to the reverse switching opera- 
-a25+ Ww tion of the relay. These quantities are determined graph- 
d ically. Further, equations for the next operational sections 
-050 are made up, similar to (22) and (23), The values obtained 
for 04, Og, Og and o are shown in Fig, 5. 
Fig. 5 
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We are interested in the output parameter og) of the system, which represents the degree of completion of 
the course run by functiong@). Comparing Fig. 5.c with curve II of Fig. 3, we see that at the beginning a complete v 
compensation of the lagging is taking place (no shift by T present). After that, however, some difference is observed 
between these two curves, and from there on, the control action o(®) assumes the character of undamped oscillations, 


This is due to the fact that one of the conditions for achievement of complete compensation is not satisfied, namely, W 
the requirement that the time interval between consecutive switching operations of the relay must be greater than the ¢ 
lag time T: already the time between the first and second switchover of the relay is smaller thant. During the tran- i 
sient process this incongruity is repeated many times, and finally exerts a predominant influence. During this time F 
the output function og fluctuates with a small amplitude and at a relatively high frequency. The latter fact is in 
contradiction with the statements of the author of papers [1] and [2): the suggested method of stabilization does not a 
eliminate the problem of wear of the relay contacts. 

An attempt to investigate compensation of time lag by the phase plane method in a manner, similar to the in- t 
vestigation of the lag itself, as for instance in [4], does not lead to desired results; the basic virtue of the phase plane P 
method easier determination of operational complexes of a system by means of graphical construction and plotting ¢ 


of switchover lines and point loci transformations—can not be made use of in the present case, because the control 
operation if governed to a great extent by whichever derivative of the system happens to be the “dominant” deriva- 
tive, and it is not knowna priori which derivative that is; this can be determined only by studying a whole family of 
switchover curves, and selecting the one having the proper crossover with the curve of the point loci, which in turn 
is dependent upon the initial conditions, The only thing left to do is to use the method of sectional curve fitting, 
that is, to solve the equations of motions in several consecutive sections for which the initial conditions are given 




















Fig. 6 


(Fig. 6), Therefore, representations of this operation on the phase plane does not offer any advantages in comparison 
with the representation by an explicit time function (Fig. 5), either in simplicity or in clarity. 


4. Compensation for Hysteresis Loop of Relay Characteristics 


Until now only single valued relay characteristics were considered (see Fig. 2), However, characteristics of ac- 
tual relays have sections consisting of hysteresis loops (Fig. 7, a and b), In order to make actual use of the described 
stabilization against timelag effects, the hysteresis loop sections of the relay characteristics must also be compensa- 





ted for, 
Let us examine first the characteristics of a relay shown in Fig. 7, a these characteristics are described by the 
following equation: ‘ o>. 
b>so>—b yw=i 


y=9(s) = 





s>o>—4, You—ti 
—1{ (24) 


’><, 


where yo is the value of y just prior to one of the switching over operations of the relay. Let us introduce a new 
variable x acting at the input to the relay and subject to the following condition: 


‘ % = o—bdsign  (X%) = o — 8M (X). sss 

. 

f which expresses the fact that there is a rigid feedback provided from the output (Fig. 8a). The operational threshold 
e of the relay is + 5 (Fig, 7a): that is, the value of x must fall to at least x = -5, or rise to at least x = +5 for the re- 
- lay to carry out a switching operation. In either case, until this threshold is reached, o will equal zero (Fig. 7, b). 


However, such a compensation of the hysteresis loop by a simple rigid feedback across the relay itself results in faulty 
operation of the latter, because, as is well known, at the moment of the switchover the input signal also is subject to 
a jump, which moreover is equal to +26, 


In order to avoid contact chattering during the makes and breaks of the relay, it is sufficient to introduce in 
, the feedback circuit a small delay action, Inasmuch as there is already a time delay T present in the system, it ap- 
. pears advantageous to make use of it for this purpose, To do this, if the component causing the time delay is in the 
circuit after the relay, part of the signal is fed back to the relay input from the output of the time delay component 
(Fig. 8,b),if the time delay element precedes the relay in the 
circuit, the feedback goes from the output of the relay to the 
input of the time delay component (Fig. 8,c), It must be noted, 
Yo) PX) that complete compensation for the hysteresis loop will take 
t . tg ad place only in case the interval between successive switching 
° / operations of the relay is not less thant. When feedback not 
a | across the relay itself is used, but across the relay plus the 
other delay causing component, equation (25) takes the fol- 
low ing form: 








0 
— 
5 

















2 a % = 6 — 8, (X), 


re 
| 1 
A... 
0 g | 0 a é 
eav bh nade Substituting into (26) the value of o from (13) and tak- 
: ing into account (17) and also the fact that to the input of the 
relay now the signal -x instead of -o is applied, we obtain 
the following: 
X= (0) + 19(0) + aD, (—%) +80, (—%) = 
= 9 (0) + 79 (8) —(a + 8) ®, (X). 7" 


We see that the compensation for the relay hysteresis loop, mathematically, is due to an increase of the coeffi- 
cient & of the additional feedback by the relay hysteresis 5 (see Fig. 7, a). 


~~ 

s 
~ 

- 


(26) 
d H where the subscript T , at ® symbolically designates that now 
/ 








the delayed byt value is used as the independent variable in 
function @, that is 4[x(9-T)] =, (x). 























Fig. 7 
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Fig. 8 


The compensation of a system having a relay characteristic of a more general type (Fig. 7,c) is carried out in 
a similar manner. The additional feedback around the relay and the delay components is expressed by an equation 
similar to (26), namely: 


% = 6— vO, (X). (28) 


Substituting into (28) the proper value for o, we obtaim 











% = 9 (8) + rp (0) — (a + v) D(x). 
(29) 


We see that the compensation for relay hysteresis amounts to an increase of the feedback coefficient « by the 
value of the relay hysteresis v. The effective dead area of the relay (when the input to it is below the threshold) 
with this hysteresis compensation is equal to € + v, It should be noted that in [1] it is erroneously stated that the 
feedback coeffitient is decreased by the amount of hysteresis, and that the dead area is increased not by v, but by 
v/2, 


SUMMARY 


An effective system for stabilization of relay systems subject to time delays is described which was suggested 
by R, W. Bass [1,2]. The time lag of the system and also the effects of the hysteresis loop sections in the relay char- 
acteristic are compensated by additional control action of a kind specifically developed for this purpose, 


The method is applicable to a wide variety of relay systems, and particularly in cases where the delay is an in- 
herent property of the relay itself. The wearing out of the relay contacts, however, may be increased by this method, 
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INFLUENCE OF HYSTERESIS ON THE MODE OF PERIODIC 
PROCESSES IN PULSED RELAY SYSTEMS 
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The approximate method of defining periodic modes of operation in automatic pulsed relay systems 
described in [1], is extended to include a pulse-unit relay characteristic having hysteresis, 


Some of the problems of the dynamics of pulsedrelay systems having an ideal relay characteristic have been 
examined in[1], The presence of hysteresis in the relay characteristic of a pulse unit may have a marked effect 
upon the mode of the processes occurring in such systems, 


The dynamics of pulsedrelay systems having both hysteresis and a second-order linear group have been investi- 
gated in (2) by the phase plane method, 


In the present work, a method is outlined in which the effect of hysteresis upon the periodic- mode parameters 
of pulsed relay systems is taken into account irrespective of the complexity of the phase plane, Fig. 1 shows the 
schematic diagram of the pulsed relay system investigated. The pulse unit of such a system can be represented in 
the form of a two-position relay with hysteresis, series connected to an amplitude- modulated pulse unit. 

















Let the system have simple self- oscillations with a relative half-period N, where N is an integer, Then, in 
analogy with continuous relay systems [3], the conditions for the existence of periodic solutions in a pulsed relay sy- 
stem may be written in the following form: 


Z(t) =%) at (N—1)<i,<N 
z(,)>0 a (N—1)<i<N 1) 
Zs (t) < %, and =%4(t)=—1 at O<St<N. 


Here, %» is the half-width of the hysteresis ZOne; t= t/T,, is the relative time; to is the moment of switching, 
that is the moment at which occurs the intersection of Xt) with the line %»; N is equal to /Tps T is a half-period 
of the oscillations; and T), is the repetition rate of the pulse unit. 


The relations (1) constitute the conditions for the proper moment and direction of switching as well as the con- 
ditions for the absence of switching throughout a half-period, It is assumed that the linear group of the system has 
the properties of a low- frequency filter and, hence, the stable response ¢ of the linear group to the first harmonic of 
the periodic pulse train, XA t), may in approximation be substitued for x Xx). Here, the conditions for the existence 
of periodic solutions may be expressed through the values of X<0 and Xx t) in the discrete points N - 1 and N, in the 
form of the follow ing two pairs of inequalities: 


is(N—1)<x.  %(N—1)>0 


~ x (2) 
t3(N) > x,. t(N) >0. 


At %»= 0, we arrive at the conditions cited in [1]. 


Fig. 2 illustrates the two end positions of the ‘X¢t) 















































Ne a Rather ig Bah 7 curve relative to the pulse-period moments N-1 and N, 
z,(i)\ ff £,(8) w3I ! Linear z,(t) which provide for the existence of a periodic mode with a 
\ i group half-period N=1, It may be seen from the figure that the 
ee a ae first pair of inequalities (2) is indicative of the existence of 
a periodic solution for the case where the X(t) curves are 








located on the left side of the dotted curve; the second pair 























_ of inequalities defines a given periodic mode for the case 
where the xg(t) curves are on the right side of the dotted 
7,3) durve, The existence of periodic solutions with a half-peri- 
th Pa od N =2 is defined by the complete set of inequalities (2), 
a N-0 As has been shown in [1], the output of the system and its 
N-I}+0 7 xy first derivative pe tne form: 
(N-1)-0 ,/ “x, P ses 
ran Kee z(t) = *% MK [(2m —1) w| sin {(2m —1) +t 
}-* - 4k , 
W-1 N +940 [(2m—1) al »=7 3) (2m—1) MK 
Fig. 
- n 
{(2m—1) x | cos {(2m—1) 5 t+-- o[(2m—1) 9} 
where 
sin (2m — 1) = ¥ 
Mins + oN and p= (2m—1) oy (1 —7)- 





7 sin (2m — 1) or 
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Hence, their first harmonics are: 


z3 (2) = tte MENT (8) gin 2 9 + awit —1) +0(+)| 
sin oN (3) 
z, (i) = the SIN (8 oul € + sy (t—1) + 0(7)]. 
sn 5 (4) 


where k, is the simplification factor of the pulse unit; y is the relative duration of the pulses acting at the output 
of the pulse unit; K(/N) is the frequency response modulus of the linear group of the system, and 6( 1 / N) is the 
linear- group phase, 


Substituting the time values t = N-1 and t = N into expressions (3) and (4), we get 
pare: | pa xn i+y7 cid 
; Z,(N — 1) A sin|-F- —ft_ | 


ATi) #,(N) = Asin[—y “>t —0( 














%(N—1)=—Fcos[ F441 g(a] © 
(WN) = — +-cos [ — Fr +zt—e +): 
4k sin 
where A=—*——~_ K (+). 

si oN 


Fig. 3, 


By subordinating the equations (5) to the inequalities (2), we obtain the conditions for the existence of a 
periodic solution for an arbitrary N, as follows: 


ll baer) OG +)< 


6 
—|x (2k -+1)—n—* x aa (k= 0.1.2. .... 00). ” 
Here, n = arc sine'*/A, 


From this inequality it follows that if the point on the frequency-response curve of the linear group that corres- 
ponds to a certain value of the relative frequency ™/ Nj has a phase 6 (1m /Nj) that satisfies inequality (6), then the 
system will have a periodic solution with a relative half-period N = Nj. Geometrically, condition (6) defines the 
vectors on the complex plane which have their vertices at the origin of the coordinates, for any value of the relative 
half-period N ad the relative frequency"), The angle of sector 8 is numerically equal to the value of the rela- 
tive frequency ™ for which it is plotted, It can be seen from expression (6) that the effect of hysteresis becomes 
apparent in a deflection of each sector through an angle nj in the direction of the lower relative frequences (larger 
N). 


Since at %= const, the angle n is a function of the amplitude of the first harmonic Xt), which for given y 
and ky is dependent on frequency - alone, each sector will be turned anti-clockwise through its angle n; . 


Fig. 3 shows a method of determining 1; graphically directly from the frequency response of the system's 
linear group expressed in values of N. The frequency response modulus is multiplied by the coefficient kya, where 
a= 4/ [(sin 1/ 2N)y / sin / 2N), Values for a calculated for variousy and N, are given in [1]. 


Since the position of the sectors depends upon the amplitude of the periodic mode of the system, a variation 
in circuit gain—as distinct from a system employing an ideal relay—will cause a change in frequency of the perio- 
dic process, this change being of a discrete nature. An angular position of the sector that determines the existence 











in the system of a periodic mode with a half-period N =Nj can vary in the positive direction through an angle a;, 
and in the negative direction through an angle a, without affecting the frequency of the periodic mode involved 
(Fig. 4). Thus, a change in system gain over a certain range will not affect the frequency of the periodic process. 
At larger gain variations in one or the other direction, the oscillation frequency of the system will vary in a jump- 
like manner, Let us denote by nj = 14 + a; the limit value of angle n at which the system transfers from the exist- 
ing periodic mode to another, lower-frequency mode, Here,it is necessary always to observe the condition ny = 
= , which expresses the fact that the amplitude of the periodic mode may not be less than %», Bearing in mind that 


as . Xo \ Xo 
= arc Sin = ame 
" 1 a arc sin A; + &% 
where Af is the amplitude value of the oscillations having a relative half-period N = N,, at which there occurs a 
change in frequency, we find the corresponding change in system gain which provides for the occurrence of a pericd- 
ic mode with a relative half-period larger than N;- This change is equal to 





4B = 
A %o+ A, sin a ° 


From condition nj = 1/2 it follows that the magnitude of angle a, must obey the inequality 


a1 < are sin (4 — “t) + 


The boundary condition 1, at which the frequency of the periodic mode increases, we will denote through Te 
= Nj - %g, with the restriction that an increase in gain does not render the angle n smaller than zero, i.e,, m4= 9, 
Since 


aide ice. os . 
We wosn~- ware sin A, 
the relative increase in gain which leads to an increase in periodic-mode frequency is equal to 


es Xo 
~~ Xq— Aysin a © 





>| > 


The feasibility of such change is limited by the condition 


Oe < 1,- 


In this way, the presence of hysteresis in the pulse-unit relay characteristic, in the case of a limited value of 
open-system gain, leads to a drop in frequency of the periodic mode, and makes it dependent of the gain value 
selected, 
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Three variations of an automatic search method using sufficient conditions for the extremum are 
described, Sufficient conditions for the extremum are found which are suitable for use in computers. 


In solving problems about the automatic search for the extremum of a function of several variables, methods 
which involve the necessary conditions for the extremum are ordinarily used (see, for example [1-3]), The search 
ceases at point My if each variable individually has an extremum at this point, This condition is necessary, but not 
sufficient, Such methods give the desired results for a large class of functions, but they are, strictly speaking, appli- 
cable only to the case n= 1 (n is the number of variables), Even when n= 2, convergence to a saddle or crest of a 
function Q, for which the extremum is sought, is possible. This happens because certain search methods, which theo- 
retically converge strictly to the true extremum point if the function is continuously differentiable, in practice do 
not converge to the extremum since they use a finite step. This proves that if the step is finite, it is impossible to 
consider the convergence of the method used only for the case of continuously differentiable functions, Practice 
shows that it is better to consider the function Q continuous, but differentiable only almost everywhere, and to use 
only those analysis results which are justified for such functions, For example, Fig, 1, 4 shows the movement of an 
automatic search system for two methods: the Gauss- Zeidel method (heavy solid line), and the gradient method 
(heavy dashed line), The light lines designate equal levels of the function Q. The search converges to points K; 
and K,, respectively. There is an extremum for each of the variables individually at each of these points, The 
points are not extremum, however, since the sufficient conditions for the extreme are not met for the aggregate of 
the variables, 


Methods are described in the literature which permit one to overcome this difficulty to a greater or lesser 
degree [4]. 


Three variations of one search method are described below, which ensure convergence to the point Mo, at which 
the sufficient conditions for the extremum for the aggregate of variables are met, and which hold for a broad class of 
functions, 


The essence of the simplest variation will be explained by an example, Let it be necessary to find the extre- 
mum of a function Q = Q(x;, xg, ..., Xp) Of n variables, We will proceed in the following manner, 


1, We will seek the extremum of Q in the usual manner. However, an extremum may or may not exist at the 
convergence point of the search. 


2, We will fix the variable x;. 


3, We will seek, by any known method, the extremum of Q for the remaining variables, not concerning our- 
selves about meeting the sufficient conditions, 


4. We will remember the value of Q at the point found, 
5. We will change x, by a small amount 4x,, leaving all the other coordinates unchanged, 
6, We will repeat the search for the extremum of Q with this new value of x;. 




















7. We will compare the former and newly found values of Q, On the basis of this comparison, we will make 
the next change in x, in such a way that Q changes in the desired direction, The path of movement for the case of 
two variables is given in Fig, 1, b. 


This method can be called "the method of fixed coordinates.” 
It can be shown that this method converges to the point My, which is 
extremal for a certain class of functions Q in which the convergence 
points of the search in n-1 variables form curves in the coordinate 
space, for different values of x,, regardless of whether Mo lies in the 
b domain of the given function or on the boundary. 


The block diagram of a system which operates according to 
the method described is given in Fig. 2, The output of object O is 
fed to two automatic optimizers: a multichannel optimizer AO,, and 
a single-channel optimizer AO;, The variable x, is set by optimizer 
AO, after optimizer AO, finishes its search, The remaining variables 
Xe, Xge++ +» Xp are set by optimizer AO,, AO, is a step optimizer; it 
takes a step on command from AO,, which acts after the completion 
of a search, If changes in x; in one direction cause undesirable changes 
in Q, optimizer AO, changes x, in the other direction, AO, can 
Fig. 1 be any single-channel optimizer, The simplest circuit of such an 
optimizer is shown in Fig, 3, 








To explain the operation of the circuit, let us first assume that capacitor C, and the circuit with Ps are not 
present, and that R, is very small, P, is a two-winding neutral polarized relay, connected so that the passage of cur- 
rent through coil A causes contact 1P, to close, operates P,, and transfers the moving contact 1P, to position B, Cur- 
rent now flows through winding B of relay Py. This causes contact 1P, to open, etc, A pulse pair is obtained as a 
result, If R; is increased, the pulse-pair oscillations will cease when a certain value of R, is reached, If a short posi- 
tive pulse is fed through C,, one pulse pair will occur, and the process will once more stop, The contact 2P, changes 
its position in this manner, The first positive pulse switches contact 1P3, for example to +E; and the second to -E, 
Contacts 1P,, and 2P,, are closed by optimizer AO, when a search is finished, Capacitor C, is charged to a voltage 
equal to Q, and capacitor Cs to a voltage E, The polarized relay Ps is connected so that the charging current does 
not cause it to operate, After contacts 1P,, and 2P, drop out, there is no discharge current through C, (the capaci- 
tor “remembers” the value of Q), but the Secheree’ current = Cs flows through relay Ps and causes contact 1P, to 
close briefly; x, is increased or decreased in this case, depending on the position of 2Ps, The search for the extre- 
mum of Q is conducted for the variables x2, Xg,..+» Xn» 
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After the search is finished, relay Po built into AO,, operates again. If Q has decreased, no pulse pair is 
triggered, and contact 1Ps causes x, to change in the former direction. If, however, Q has increased, the pulse pair 
is wiggered, and x, will change in a new direction after contact 2A, drops out, The optimizer described was tested 
and “felt” a difference of AQ = 0,2 volt. If it were necessary to search for the maximum of Q, it would be neces- 
sary to change the polarity of E and to reconnect the windings of the polarized relay. 


The suggested method and the conditions of its applicability have a mathematical basis, which is given in 
the Appendix, 

















The proposed method is very simple. However, if an arbitrary variable is fixed, it does not work for all func- 
tions, It gives good results for a number of functions if the variable tc be fixed is chosen correctly. In this case, the 
optimizer can automatically switch the output of AO, to the different variables x;,x2,..., X, in turn, by means of a 
special commutator, Functions may be encountered, however, for which even this complication of the search(second 
variation) does not ensure convergence to the true extremum. 


In this case, a method (third variation) can be recommended, whose essence will be explained by the method 
of induction, 


For the Case of Two Variables 





1, Using any known method, we will find the extremum of a function of two variables, not concerning our 
selves about meeting the sufficient conditions, 


2. We will fix the variable x. 
. We will find the extremum of a function of one variable, x, (sufficient conditions are met). 


. We will remember the value of the function Q at the point found, 


. We then repeat the search for the extremum of Q with x, as the variable, 





3 
4 
5. We will change x, by a small amount 4X;,, leaving all other coordinates unchanged. 
6 
7 


. We will compare the former and newly found values, and make the next change in x, in the proper direc- 
tion on the basis of this comparison, 


For the Case of Three Variables 


The variable xg must be fixed, and a point sought by the above-mentioned method at which the sufficient con- 
ditions for the extremum of a function of two variables (x,;, x.) are met. Then xs is changed, and the search is re- 
peated with x; and x, as variables, and finally the effectiveness of the change in xg is evaluated, 





For the Case of n Variables 


The variable xp must be fixed, and a point found at which the sufficient conditions of the extremum are met 
for the (n-1)st variable, and the value of Q found at that point must be remembered, Then xp can be changed, the 
search repeated in terms of the (n-1)st variable, and the effectiveness of the change in xp evaluated, The mathema- 
tical basis and conditions of applicability of the suggested search method can be found in the first half of the theorem 
in the Appendix (for the case n =2), if one considers x, to be a vector with coordinates (x2, X3,..., Xp)- In this case, 
a line of extremums will exist for a very broad class of functions. 





We note in conclusion that the use of these methods causes practically no increase in difficulty or search time 
(even the last, most complex method), The fact is that the variations of the method suggested do not supplant other 
search methods, but only supplement them and one another, If a simple method gives the proper results, there is no 
sense in using a more complex one, and if the simple method does not give correct results, it cannot be compared to 
a more complex one which does, Thus, for example, a search for the minimum of the function z= x’ + y’ is easily 
accomplished by any means, The tuning of thyrite square-law devices was greatly hastened when the first of the 
suggested variations was used, The search for the minimum of the function y = max{ x4, Xp, oss Xn} » Which easily 
converges when the third search variation is used, does not converge to the extremum when the gradient method, 
steepest slope method, Gauss-Zeidel method, or the two simpler variations of the suggested method are used, The 
information presented may prove useful even in the case of manual search for the extremum, 


APPENDIX 


Let us find the sufficient conditions for the maximum® of a function of n variables, suitable for use in compu- 
ters, The meeting of sufficient conditions for the maximum of a function of one variable is easily confirmed by an 
automatic machine, It is desirable to use this fact in solving the stated problem, 


* The result is found in an analogous manner for the minimum, 
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In practice, it is important to find not the strict maximum, but the largest value in a certain, limited, closed 
region. If the function is continuous, this value may either be a maximum, or may lie on the boundary, It is desir- 
able to combine these two cases by the term “limiting value,* 


Definition I, The function Q = Q(x;, Xg,..+» Xn), given in the region G(x;, X2,..., Xp), is said to have a 
limiting value at the point Melxy, X20, ++. Xnot> if 1) in any arbitrarily small area around the point Mo, there are 
points belonging to G(x;, Xg,...+, Xp); 2) am area around Mp exists in which the value of Q(M)<Q(Ms) at all points M be- 
longing to G(x), Xg,.++, Xp)- 


Let us now examine the function Q(x;, X2,..., Xp) at the point M{x%4, Kees Xeee sees Xnel - Here x,= Xte and 
is arbitrary and fixed, while Xge, Xge,..., Xne are certain numbers, 


Definition Il. The point M{x}, Xse, Xges-++» Xne} will be called a point of partial maximums for x, = Xts if 
Xees Xgee+++» Xne are such that each of the functions: 


QAX2) = Ax}, Xee Xge,-++ + Xne) Of one variable x», 
Qg Xg) = Axf, Xge, X34-++e "Xne) Of one variable xz, 


QnlXp) = Axf , Xe, Xge,-++» Xp) Of one variable x, has a limiting value at the 
point M, 
It is evident that 
70> Pa (21). 
Z3e> Ps (71), 


(A) 


Tn es Pn (21), 


that is, the position of the point of partial maximums depends on x, and on the form of the function Q, At least one 
point of partial maximums corresponds to each x;, 


Definition II], The curve given by relations (A) in the region G(x,, X,...,» Xp) is called the curve of partial 
maximums, 





A curve of partial maximums sometimes exists for continuous functions in a bounded, closed region, Another 
variable can be used instead of x; to construct the curve of partial maximums, Our proof does not depend on which 
particular variable is used for this construction. If a partial maximum curve does not exist for a given function (but 
a surface does), the theorem does not hold, We note, however, that the partial maximum curve exists for the case 
n= 2 for a very broad class of functions, This fact can easily be used as the basis for the third variation of the sug- 
gested method, 


Theorem, If: 1) Q(x, Xg,...,» Xp) is a continuous function of n variables having a continuous partial maxi- 
mum line Ly,; 2) relations (A) are the equation of a curve of partial maximums; 3) Mol x10, Xg0,--+» Xingl is @ 
point lying on this curve; and 4) there is a limiting value of the function Q¢x,) = Q[x;, %g. (X1), 3 (Xp... . Ok XD) 
of one variable x, at point Mo, then a function Q of n variables has a limiting value at Mo. 


We will prove the theorem for n= 2, and generalize the result by the method of mathematical induction, 


The Case n= 2. According to the definition of Lx,, a number 5{x,)> 0 exists for each x, such that if the 
inequality |x, <7 {x,<5 {x;) holds, the inequality Q(x;, Xp¢)>Q(x;, Xq) will exist everywhere in G, 

Since 5 {x,)>0, it has a lower bound, that is, there is a lower bound 5,; 5, >0, since 54x,) 0 for every x;, 
and the region is closed, 


Since Q,(x;) = O[x;, ¢Xx,)], and has a limiting value at Mg, a 5 PO exists such that if the inequality | x, - 
X19| < 54, holds, the conditions Qj, (x1) <Qq(x19) = AX10, Xgo). 





Finally, we find that when | x; - xyl< 5, and | x - Xge| <5g, the inequality is justified, which proves the 
first part of our theorem; 


Q (21, 22) <Q (#10, 220). 





The case of Arbitrary n, Let us examine the function Q* = Q[x,, Xg, «+++ Xp-1e 9p X1)) instead of Q, This is 
a function of the (n-1)st variable, which meets the conditions of the theorem at Mp. By induction, there exist posi- 
tive 54, 52, ++, 5n-q such that the condition 





Qazi. 2. .... En _4> Pn (%1)] <Q (Mo) 


| 2; — 240 |< 1: | ts — 220 | < de: . «+2 1S __y — Fn—yy 9 | << Ona (1) 
hold simultaneously, | 
Since Q(x), Xg, «++» Xp) is continuous in x, everywhere in G for the case when 


| 1 — 2191 + | t2 — 220 | + --- + 12,4 —2n—1) 0! > 9 (2) 


according to the lemma about the conservation of sign, a 5 (Xj, Xg..++» Xp-4), can be found such that if relations (1), 
(2), and| Xp - Xn_l <5 pn, hold simultaneously, the inequality 


Q (2. C? Z,) < Q (219, Zag. +++: Zno) (3) 
also holds, 


Here §,709- the lower bound for 6 1 X;, Xg.+++» Xp-1) exists, since 57(X,, Xg,.++» Xp-1) > 0 and the region is 
closed, If the condition 


| 2 —210|+|2%3— 20/1 +... + |p, — %n—1) 9 | = 9 


occurs instead of (2), then a 5, can be found according to the definition of Lx; such that (3) always holds when 
Xn - Xnol <4) ; 


Thus, the existence of an area around point My in which (3) holds everywhere is proven, The bounds of this 
area are defined by relations (1), and |xp - Xpel<min {5p, 6p}. The theorem has been completely proven. 
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A servomechanism with an asynchronous motor operating on alternating current is examined, in 
which the control winding is powered by a current generator. An expression is found for the output 
magnitude of the open-loop system, taking initial conditions into account, and also for the transfer 
function, The stability of the system is studied with certain simplifying assumptions, 


In AC servomechanisms, the actuating element is often a two-phase motor with a squirrel-cage rotor, One of 
the stator windings— the excitation winding— is connected to a voltage of constant amplitude; the other winding— the 
control winding —is part of the control circuit of the servo system, It is ordinarily assumed that the excitation wind- 
ing resistance is small in comparison with the inductive reactance, 


The control winding can be powered by a controllable power amplifier having a sufficiently small resistance 
(from a voltage generator), It is then assumed that the control winding circuit can be neglected, and only the induc- 
tive reactance considered, Servomechanisms with asynchronous motors operating under these conditions have been 
examined more than once [1-3], 


The control winding can also be powered by a power amplifier with a high resistance, so high that the induc- 
tive reactance can be neglected in the control winding circuit, This kind of power amplifier is a current generator. 
A servomechanism with an asynchronous motor powered by a current generator will be discussed in this article, 


1, Motion Equation and Transfer Function of the Servo System 





Let us examine the servo system whose diagram is shown in the figure. In this diagram, I is a measuring ele- 
ment, M is a modulator, 1, 2, 3, are linear passive four-pole networks, A is an asynchronous motor, TG is an asyn- 
chronous tachometer generator, and P is a gear box, Using these designations, we can write an equation for the 
measuring element 


x(t) = ©, (t) — 93 (¢) 


(1) 
an equation for the modulator 

Uz (t) = Umx(t) cos (@ot + 9) (2) 

and one for the voltage across the tach generator output 

k ; 
ug (t) = £0, msin (wot + a) S. (3) 
Neglecting the resistance of the excitation winding circuit, we write the equation of this winding 
ap 

Us (t) = U om 8in (Wot -}- p) = a> - (4) 


Here ¥, is the flux of the excitation winding, 
Let the transfer functions of four-poles 1, 2, and 3 be Yx(p), Yfp) and Yg{p) respectively, 











Using the diagram, we can write 





le (p) = 2 1¥,) Ux () —¥20) & i. (5) 


Here R is a constant having the sense and dimensions of resistance, I,(p) is a representation of the current in 
the control winding. Equation (5) can be considered as the control winding circuit equation, 
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Let us write an expression for the rotor flux referred to the axis of the control winding. 
Ya = Lig + My .. (6) 


Here L is the inductance of the rotor, My is the mutual inductance between the rotor and control winding, i, 
is the rotor current, referred to the control winding axis, 


The rotor equations can be represented in the form 


ap, d® 
a + 1% a+ ter = (7) 
ms 


— MPa 5 2 + ir =0 

(7a) 
where r is the rotor resistance, and n is the transfer coefficient (the ratio of the angular rotation of the motor shaft 
to that of the drive shaft, ©,). Equations (4, 6, 7 and 7a) are written under the assumption that one can neglect the 
rotor and stator leakage flux. 


Let us write the equation for the rotation of the drive shaft 


de d8. 
ee a a —Ya ]--4 a 








Jn? SS +. [n (ape + a) ci 

Here J is the moment of inertia of the rotating mass, referred to the motor shaft, and k d©, /dt is the external 
torque applied to the drive shaft, If the term vs (d6,/ dt)in the left half of equation (8) is neglected, one can say 
that for small signals, equations (1-6, 7 and 8) form a set of linear differential equations with periodic coefficient, 
from which equation (5) can be written in a form suitable for Laplace transformation, This set can be solved by 
means of the method given in [1,3]. Writing (1-6, 7 and 8) in a form which takes initial conditions into account, 
and solving for Op), discarding terms involving X and ©, as functions of the operator P with arguments p+ jw» on 
the basis of the reasoning given in [1,3], we get the following expression for ©,(p): 
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Here the designation [ ]* represents the complex conjugate of the expression included within the brackets [ } 
appearing before it, 


In equation (9), the operator coefficient for X(p) is the transfer function, and the remaining terms express the 
effect of the initial conditions, Introducing the operator s = pAvg, and the constants 
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we write the transfer function of the open-loop system in the following manner: 
O2(s) _ 
X (8) 
1 { — 2 — % ; 1 
' 2 alccesies re—Ave—D"|+| |} (11) 
ge j (s — j2) | ; \_Ts (72 j2) 
Ts <T -— wr Y 
{ wHt+s (tt (s—))+1 Yi (s—) Ya(s—y)e | ) =| ( Tis— Srl I" } 
Using equation (1), we can find the transfer function of the closed-loop system, 
If the voltage 
d 
ty (t) = = USm COS (Wot + &) >= ot 
(3a) 
is put into the feedback loop inst-ad of (3), we get the expression 
s—j2 - ; . 
62 (s) ey » Yi zilrw= per Ti(s—j+i esi ea. I*} 
¥ (s) (11a) 











T (s—j2) 
1 y — pM! am ° . 
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for the transfer function of the open-loop system, 
It is evident that the transfer function of the closed loop system can be found from (11a) by using equation(1), 


2, Stabilization of the Servo System 





Je< us examine the servo without feedback, in which case the only member with an operator coefficient is the 
two-phase asynchronous motor, 


Substituting K, = 0, Ys) = Ky, and Y{s) = Kg into. equation (11) or (11a), where K, and Kg are amplification 
coefficients, and assuming for simplicity that ¥ = ¢, we get the following expression for the transfer function of the 
closed-loop system 
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K,K 
K, = See. 





Let us examine the case of a purely inertial Joad(k= 0), Then T= 1, as follows from(10), The inequalities 
T*» 1 and TT, » 1 ordinarily hold well for an asynchronous motor, Keeping this in mind, and considering the ap- 


proximate transfer function [2] for|s*| < 1, we find, from (12), 


rs) _ Ko (s + 2T) 
G(s) TT ys® + TT s* + (Ko— T*) 8 + 27K, 





(12a) 


Using the Routh- Hurwitz criterion, we find that the system described by equation (12a) is unstable, 
Let us consider the system with feedback, under the condition that the voltage of expression (3) is introduced 


into the feedback loop, 
Let four-poles 1, 2 and 3 be inertialess with amplification coefficients Y¢s) = K;.,¥4s) = K,, and Y(s) = Kg, 
Setting 5= y= 0 for simplicity, we find, by using (11), an expression for the transfer function of the closed-loop 


system 





O2(s) _ Ko [T (s* + 2) + 8] 
(s) 8 (T+ 1) (7? (s*+ 1) + 27's + 1) + p (Ts + 2) + (Ko— 77) (T (8? + 2) +] * 
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Here 
p = Ky KK). (13) 


If again we set T=1, T?> 1, TT, > 1, and examine the approximate transfer function for|s*| << 1, we get 


(s) Ko(s + 2T) 
Oi(s) AP Tis + (777, + pT) 8? + (Ko + 2p — T) 8 + 27K, * (13a) 
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The servo will be stable if the following conditions hold; 
Ky +- 2p >T* (14) 
p (Ay + 2p + 277, — 1°) > TT, 43K, + 7%). i 


It is evident that inequalities (14) and (15) will hold if is sufficiently large. 


Let the voltage of (3a) be introduced into the feedback loop, Then under the previous conditions (5 = 
YAs) = Ky Ys) = Kg; Y9(s) = Kg), by using equation (11a), we get an expression for the transfer function of onsen 


loop system 


y = 9; 


6x() Ko [T (s* + 2) +s] 
Sie) 8 (T8 F 1) (T*( $1) + 27s $1) + 1K +e — TTT +2) +8) (16) 





If we again set Ty = 1, T* > 1, TT, > 1, and examine the approximate transfer function for|s*| « 1, we get 








Ole) _ Ky (s + 27) 
6:0) PT F(T Fw) + (Kot opt — 1) 0 + aT Ky ts 
The servo will be stable if the following conditions hold 
Ky + 2pT >T? an 
K 
of $+ 28 +7 (277, —19)]>77, (3K, +7. je 











It is evident that inequalities (17) and (18) will hold if u is sufficiently large, 
A comparison of conditions (14) and (17) shows that (17) will hold for smaller values of , 
Let us now compare conditions (15) and (18), Since the right halves of the inequalities are identical, we will 
examine the inequality 
K ‘ 
- +2447 (2TT,—T*)> K, + 2p + 277, —T? 
from which we find 
K, 
27, —T> 7 - (19) 


If (19) is to hold, it is better to introduce the voltage of (3a) into the feedback loop in order to stabilize the 
system, 


The quality of the servo can be enhanced by putting a differentiating four-pole in the loop, Since. the system 
is unstable without feedback, it is desirable to use feedback to gain stability, and to introduce a differentiating four- 
pole into the loop, thus enhancing the system quality, Let such a four-pole be the four-pole 3 (see figure) with 
transfer function 





mm s? +-kys +1 
Y,(s) = Ky Stink? 
Let us examine a system in which the voltage of (3) is introduced into the feedback loop, We will seté = y 
= 0 for simplicity, We then get 








Oi(s) Ko (2s + hy) (28 +- hea) (8 4+- 27°) — (hp — hy) (T's + 2) 5°] (20) 
Yis 8 Ts +1) QT, + T? +1) +p (Ts + 2) —TT3 (27 + 8)) (2s + ha) 
for the approximate transfer function of the open-loop system for |s*| « 1 
If the voltage of expression (3a) is used in the feedback loop, we get 
Gals) Ko (28 hy) (28 + ha) (8 4+ 27) — (hp — bs) (T's + 2) 8°] 
XG) Ts +1) (2Ts—T* + 1) + 2T + 8) (@—TTs)) Cs + hay (21) 


under the same simplifications assumed for the transfer function of (20), 


If the diminuend in the numerator within the square brackets of expressions (20) and (21) can be neglected, 
the simplified transfer functions obtained can be used for a preliminary determination of the coefficients k, and k, 
of the differentiating four-pole, and of the coefficient which characterizes the feedback, 


SUMMARY 


An expression for the output magnitude was found for a servomechanism operating with a two-phase AC 
asynchronous motor whose control winding is powered by a current generator, initial conditions being taken into 
account, An expression was also found for the transfer function of the system, W ith certain simplifying assumptions, 
it was shown that a system having the asynchronous motor as the only inertial element would be unstable, Stability 
conditions for the system were found as functions of the phase of the voltage acting in the feedback loop, 
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The paper presents an analytical and experimental investigation of the steady-state and transient 
responses involved in the conversion of discrete electrical pulses into mechanical angles of rotation 
in a pulsed electric drive with a three-stator step motor, Formulas are derived for choosing the 
optimal parameters of the system, and an analysis is made of certain of its operating modes, 


In automatic control systems a converter that converts discrete electrical pulses into mechanical angles of ro- 


tation is often used as an actuating drive [1-3] (Fig. 1). The converter consists of dc source, a three-stator motor of 
the step type, and a network which amplifies the power in the pulses, 
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_ Fig. 1. Circuit diagram for the amplification and conversion of discrete 
pulses to angles of rotation, 1, 2,3 are the windings and poles of the 
step motor; fr, 2, fy are the ballast resistors, 


The converter operates as follows, As the control pulses are successively applied from the input unit to the 
grids of the rectifiers (at discrete instants of time) the current is switched from one winding to another; as a conse- 
quence of this the rotor rotates through a definite angle, 


Each control pulse corresponds to an angle of rotation equal to '/; of a pole division, This angle is henceforth 
called the “step,” 
In the literature, we are acquainted with [4,5] which investigate the dynamics of step motors, However, these 
investigations were performed from the point of view of the mechanical behavior of a rotor in a motor without taking 
into accountthe processes which occur in the pulse amplification network and the pulse conversion device, 


In this paper we investigate the entire converter as a whole, taking into account both the mechanical and 
electromagnetic processes, 

















The investigation of the processes in the converter is a complex nonlinear problem in view of the presence of 
the thyratron rectifiers and the variable inductances of the motor windings, These inductances are in the general 
case a function of the current and the rotor angle of rotation, Each stator of the motor must be treated as an electro- 


magnet operating on periodic pulses with a variable frequency, 


We shall study the energy relationships in the stator for this mode of operation, We shall introduce the follow- 
ing quantities: w is the angular frequency for the periodic pulses of arbitrary shape which are applied to the winding; 
6 is the time in angular units measured along the curve for the periodic pulses © = wt; x is the instantaneous value 
of the winding reactance corresponding to harmonic oscillations at a frequency equal to the frequency of the periodic 
pulses: 


z=oL 


where L is the instantaneous value of the motor winding inductance, 


We shall assume that the reactance of the winding is a function of the current i and the angle of rotation ¢ of 
the rotor: 


x= x(i, @). 


We shall make use of the known relationships for the voltage balance and the energy storage in the winding, 
Taking into account the definitions introduced above, these equations can be written as 


Note that in the latter equation the upper limits of the integrals, as well as the expressions in the integrand, 
depend on the current and the angle of rotation of the rotors. 


Taking into account the fact that 
d (zi) 9 (zi) di 0 (zi) dp 


a ~ oF d+ dp dQ 


we find the power drawn by the winding 





sui = tr + i 2 a (ei) e j a (zt) fet 


(1) 


Differentiating the expression for the stored energy, we obtain 


dw Ow di 
-S-Gaete a: 
Since 
dw 1 @ ( i (zi 
Tw wH [i(2i)— | ziai] = < oe) 
it follows that . 
8 d dw 
q=i PR toRs. (2) 


By analogy with circuit theory we shall call the quantity q the instantaneous reactive power for a harmonic 
input. Such a definition is justified on two counts, First, Expression (2) contains the winding reactance correspond- 
ing to a definite frequency of the harmonic current. Second, this expression corresponds only to the magnetizing 
power producing the magnetic flux, 











When there is no saturation 


io (;d@2 , _ di\dg _ fi dz di\d 
roM@tiiSt+sa le  omigHteR)H- (3) 


Studying expressions (3) simultaneously, we can write the following equation for the resistive power by analogy 


with circuit theory for a harmonic input: 
’ it dz d 
Pa= ir + 5a ae 


? (4) 
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Fig. 2, Circuit for amplifying and converting discrete pulses into 
angles of rotation, Lg is the inductance of the cathode choke; 

1, 2, 3 are the windings and poles of the step motor; A, B, C are 
uncontrolled diode rectifiers, 


After dividing (3) and (4) by i, we obtain the corresponding expressions for the total (apparent) voltage con- 
sisting of its resistive and reactive components: 


Wire tie (5) 

Ug=ir+ oo (6) 
di i 

W=2e+sa: (7) 


The resistive component of the voltage balances the voltage drop across the resistor and one half of the back 
emf, The reactive component of the total voltage balances the reactive quiescent voltage drop and the second half 
of the back emf, 


Taking the analysis above into account, it is easy to arrive at analogous conclusions when saturation is present 
if we take into account the magnetic characteristic of the motor, The resulting expressions provide a comparatively 
easy means of performing an analysis of the transient responses and energy relationships in the step converter network, 


The processes in the converter are studied in accordance with the circuit shown in Fig. 2, In contrast to known 
circuits, this circuit has an increased efficiency and better frequency characteristics; this can be explained by the 
absence of ballast resistors and the special connection of the capacitors and the diode rectifiers. 


The results of the analysis performed below can easily be extended to apply to the circuit in Fig. 1 and its 
modification which is the particular case of the variant we are studying; it can also be extended to apply to analo- 
gous circuits with silicon transistors, 











The reactance of the motor when the rotor turns can be expressed by the relationship (cf. Fig. 3) 


A 
%=Znint > (1 — cos a) (8) 


where & = 21g, 1 is the number of pairs of poles, and Ax = Xmax - Xmin: 


When (8) is taken into account 





- , Azl@i? . 
Pa = Ur + —5-— sina 
r (9) 
= zi < Asidi in @ 
q= um Ze 


(10) 
is valid where % = dy/dt is the angular velocity. 
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Fig. 3. Variation of the motor reactance 
when the rotor turns, t, is a pole division, 


Xmax 4d Xmin are the maximum and 
minimum values of the reactance, 












































Fig. 4. 4, is the current in the motor winding, i, is the 
plate current of the thyratron; u,, is the voltage across 
the thyratron, i, is the rectified current, Uy is the recti- 
fied voltage at the output of the choke; a) Ly = 7h; b) Ly 
= 0,5 h; c) Lg = 0, rg = 50 ohm, 


We shall study the steady-state mode. We shall assume that the instantaneous value of the current in the wind- 
ing does not depend on the angle of rotation of the rotor; the duration of this value is © = %4 1, (i.e,, the current has 
a rectangular or almost rectangular shape). 


For this case the equation 
6; => 6; 


is valid where 6, is the time required to process a step, 


From expression (9) we find that the instantaneous value of the electromagnetic power is equal to 





Aziéi? . 
Pe = —3,— Sin a 


Integrating this expression within the limit of a step, we obtain 


ama} 
3 Azld@i? 
Pa 2 





, 3V3Acdi , . 7‘ 
sinada = —_— i sin(a, - 3) 


a=a, 


where a, is the angle « corresponding to the initial instant at which a current pulse is applied to the winding of the 
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motor; I is the average value of the rectified current, The average value of the reactive power over a period is, as 
we know, equal to zero, 


For the electromagnetic torque we have the expression 





3V 3Ac/2 , 
M.= V —£ sin (a, + 3). 


4n@ 


(11) 
The resistive power drawn from the rectifier is equal to 
Pa=I1q Ug = Aulg + rli+ Pe + ral 








(12) 
where Au is the thyratron arc drop, Ug is the choke resistance andrq is the resistance in the winding circuit, 
From expression (12) we find 
F (U4 — Au) 4x0 
qu ; 

3V 3 lAzd sin(o1 + 3) + 4no (r + r4) (13) 

Substituting the value of Iq into Eq, (11), we finally find 

12V3 xnwAcl (0, — Au)*sin (a1 ao 3) 

es (14) 


» il [sVFeAct sin (a1 + 5) + 4x0 (r + ra) | 


In this eupeation the angle o, is a variable quantity. For synchronous operation of the motor this angle varies 
over the range Ansa,= "yr, depending on the load which is applied to the motor shaft; here the left inequality 
corresponds to maximum load, and the right inequality corresponds to no-load, At no-load it is possible to have two 
values of &, which satisfy the zero value of Me: %,2 /, anda,=- yn. The position for «,=- Ym is unstable, 


Substituting a, = hn into (14), we find the expression for the maximum (stalling) torque 


= 12V3 xnwAzl(U ,— Au)? 
max [3V3 @Acl + 4x (r + rq)" (15) 





For further loading the motor goes out of synchronism. 


Expression (14) provides a simple means for finding the dependence of the angle , on the load applied to the 
motor shaft for specified circuit reactances, motor parameters, and rotor speeds, 


The derived expression (14) is valid for comparatively low frequencies of the current pulses in the motor wind- 
ings. When the frequency of the pulses is increased the duration of the current in the windings increases, and the 
overlap of the currents in the motor windings begins to-have an appreciable effect on the magnitude of the electro- 
magnetic torque, This effect can be explained by the fact that at the instant of overlap currents exist simultaneously 


in two windings, These currents produce two torques which act in opposite directions, and this leads to lowering 
of the resultant torque on the shaft of the motor, 


The equations for the converter current can be formulated in the conventional manner used in the analysis of 
converter circuits with rectifiers (i.e., by the method of “aligning” transient responses by segments), However, the 
equations for each instantaneous electromagnetic circuit which are obtained under these conditions are very complex, 
The problem can be simplified appreciably if we determine the degree to which the choke in the converter circuit 
affects the processes in the circuit and the limits for the variation of the overlap angle between the rectifier plate 
current, 


The wave shape of the dc current as such is evidently not of interest to us, In the final analysis there is only 
one purpose for the choke in our case: to absorb the voltage drop in the converter at the instant when the commuta- 
tion of the current occurs, In the contrary case the operation of the converter would be impossible, since there would 
be no guarantee that the voltage would be removed from the rectifier during the period required for deionization, 
This time must be less than the time corresponding to the lead angle, In the general case, the lead angle 65 de- 
pends on the parameters of the entire circuit, For a specific motor and specific rectifiers the lead angle will be 


















determined by the magnitudes of the cathode inductance and capacitance, Here we find that in obtaining a specific 
lead angle an increase in the choke inductance causes a corresponding decrease in the required magnitude of the 
capacitors; a decrease in the capacitance causes a decrease in the overlap angles for the currents in the windings and 
increases the motor torque, However, the relationship here is not linear, and at a certain value the variation in the 
magnitude of the choke inductance ceases to affect the lead angle, 


Figure 4 shows oscillograms illustrating the effect of the magnitude of the choke inductance on a) the lead 
angle when the remaining converter parameters remain constant, and b) the pulsations of the rectified current from 
a Larinov rectifier circuit, From the oscillogram it is evident that for Ly = L it is possible to assume 1, = const. with- 
in the limits of a step, 


We shall demonstrate that the overlap between plate currents can also be neglected, For the commutation in- 
terval the equations for the currents and voltages can be written as follows in operator form 


Iq > ar ry = 
ra + tc, + tg + te, 

u = ir + pay, t, — tpi, (0) 
ic, = pic,u — %c,Uc (0) 

u = ists + parse 


ic, = Pic. 


The voltage drop in the thyratron rectifiers can be taken into account by subtracting it from the supply recti- 
fier voltage. Taking into account the fact that the parameters of both phases are equal, we find the solution for the 
equations cited above relative to u while assuming that i¢0) = Ig, ug(@)= Igr: 


— I gr (Pacey, + Pec + 1) 
2p (pect, + Pec +1) © 





Making the transition back to the original, we find 


ate 


2 


The results shows that for commutation the rectified current instantaneously makes the transition to the next 
rectifier without any overlap; i.e,, the duration of the rectifier plate current does not exceed 120°, 


The remarks made above permit us to treat each parallel circuit consisting of a motor winding and a capaci- 
tor as a certain nonlinear element of the over-all circuit which is subjected to rectangular current pulses with a 
duty ration equal to three, 


In order to determine the lead angle ©§ we find the voltage across the rectifier 
Ugk = Uy — Ue 
where u, is the voltage across the energized winding, u, is the voltage across the deenergized winding, 


We shall study the general case while assuming conditionally that before the motor winding is energized the 
capacitor is subjected to the voltage uc{0) whose polarity is the reverse of the source polarity (cf, Fig. 2), We shall 
introduce the substitution 


2’ _ Azl?@ . 
al Sin &. 


F = di 
u,=ir+iz +2 - 





The system of equations for determining the voltage u, is written as 
u, = (r-+2')i + pai 
ic = prcu, + tcuc(0) 


la > % 
—=i+i 
> + ic 


whence 
se U g— P2cuc (0) [(r + 2’) + pz] 
t  plpeact ple + 2) ze+hl 





Assuming 1(0) = Ig and u,(0)= Ig(r+ z"), we find 


= pez((r+ 2’) +2o(r+2'Y—-2 
ug = I 
Perot pir+2)eo+i “4 





in analogous fashion, 


After transition to the original 


21 t 
tor = Ta(r +2’) + (gre — fem) 4 
u-(0) 


+ =, Bre" — Pre**) for ¢>4 


0 
tar = (r+ 2’) la + 4 e* sin (ve + gy) — 2 sin (vt + gs)" 
for ¢< 1. 
In the expressions derived above: t =6/€ is the dimensionless time, 


B=Vaw. o=tt*. po VE. tamct VOT 


2vG v 
v=Vi-—o*. tan &=sa—7- tn h=Z- 








Taking into account the fact that for T = Ts; we have u,), = 0, we arrive at the following transcendental equa- 
tions after certain transformations: 


Be\ Ars Bi\ Arts 
(B. + Bh) e (i+ gs)e = B.— Bi for °>1: (16) 
sin (vts + @2)— +sin(vis+q,)=ve"* for <i. al 


From this Ts is found using known methods. 
We shall find the overlap between the currents in the windings, The current in the winding of the rectifier 
which is cut in is equal to 
iy = La {4— > e* [sin (vt + g) — 208in ve] | for g¢<{4 
iy = Ta {1— 5g (81 — 20) e-* — (8, — 20) eH] | for s>4. 


For T = ty,{;= Ig. However, since beyond values i, the value i; slowly approaches Ig (especially foro > 1), 
it is expedient to assume 1 = Ty for i;= 0.91y, Then t is found from the transcendental equations: 


sin (vt, + @2) + 26 sin vt, = 0,1ve*** for o<1 (18) 








(By — 26) e~** — (8, — 20) e~** = 0,4 (Bj — By) for o>. ees 


The expressions for 1 and t, contain only one parameter o which takes into account all possible relation- 


ships between r, z’, L and C, 
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Fig. 5. The dependence of the angles 
of overlap and lead on the parame‘ers 
of the converter in relative units, The 
solid lines are for uc{0) = 0, and the 
broken lines are for u,(0) # 0, 
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Fig. 6, Time variation of the plate and 
phase currents in the converter, ig, ip, 
and i, are the plate currents of the thyra- 
trons; i,, i, and ig are the currents in the 
motor windings; g, is the overlap angle 
for the currents in the motor windings, 
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Fig. 7. The computed (the solid lines) and the experi- 
mental (the broken lines) characteristic of the drive, 7 
is the efficiency of the converter, ty is the time for 
which the overlap between the current in the motor 
windings occurs, f is the frequency of the control pulses, 


Figure 5 shows the functions T,, Ts = f( 0) for 
Uc (0) = 0 and ue (0) = Igr. For a specified value of 0 
the time Tg for u, (0) #0 is greater than the time Tg 
for U, (0) = 0, and the time Ty for u,(0)# 0 is less than 
the time Ty for u, (0) = 0, 


For stable operation of the converter the lead 
angle must be greater than the angle (time) for natural 
deionization of the rectifier as indicated by its ratings, 
The connection of uncontrolled diode rectifiers into the 
circuits of the windings for o < 1 (u-,(0) 0) makes it 
possible to assure the optimal lead angle from the point 
of view of stable commutation for smaller capacitances; 
therefore it makes it possible to operate with smaller 
current overlap angles in the windings and with increased 
torque, 

It is necessary to keep in mind the fact that all of 


the derived expressions for voltages and currents when 
uncontrolled diode rectifiers are present are valid only 


during the initial interval of the process and permit us merely to find the angles of lead and overlap, A determina- 
tion of the laws governing the variation of the current and the voltages over the entire interval requires, as we know, 
the solution of the problem by segments, For transition from dimensionless quantities to the angles of lead and over- 


lap it is necessary to make use of the relationships: 


0, = t,£. 
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Since T, andt,, are functions of o, it is evident that they are functions of the velocity and the angle of 
rotation ¢, In de ng the lead angle we shall base ourselves on the following, 


The defonization time t, for the rectifier is usually specified and the lead time ts must be found while taking 
the condition ts> ts ; this condition must be maintained for all the operating modes of the converter, The motor 
parameters are also specified, Since the overlap angle increases with the capacitance, it follows that the magnitude 
of the capacitance required to satisfy the indicated inequality must be minimal, 


We shall study the expressions for o and £: 

















r+ z’ V 2r6 [2@r “+ Aci@® sin a} 
atk alae 40Y 22.45 + Az (1 — cos a) 
ya i J @l2%yun t+ 2 (1 — cos a)) 


V 5 


When the angular velocity increases, o increases, On the other hand, as o increases the quantity Ts (and 
therefore 65) increases, 


Since the lead angle 9 must be minimal but greater than critical, the magnitude of the capacitance must be 
sought for the worst case: i.e,, when 9 = 0, « = 180° and 








Spang 4 / “c (20) 
§ = V zctpax- (21) 


In order to achieve an exact determination of the overlap for the transient modes it is necessary to find o = 0 
(9,¢) andé =€ (¢), After substituting these expressions into (18) or (19), we obtain Oy = f{(>, ~). Then, substitut- 
ing this expression into the dynamics equation, we find the transient response, Note that 6, is found after 95 is 
determined [i.e,, for the case where all of the converter parameters are known (chosen) J. 


However, in the majority of cases it is possible to limit ourselves to a simplified solution of the problem, First, 
the inequality Ax/2wr « 1 is usually valid; if we assume = 180°, then the maximum error does not exceed 10 to 
20%, and this is an adequate approximation in practical computations, Second, the transient responses are usually 
determined for dynamically stable modes of motor operation (without the loss of a step) for which the relationship 
® = w/1 is valid, Under these conditions the expressions for o and — can be written as 


eas * 2a ss Zc (mint max 
rms | Pa +z = 2 , 
max min 








We shall find the expressions for the currents and the electromagnetic torque while taking the overlap into ac- 
count, Considering the fact that in the commutation interval 


di, _— dig 


a9} 


it is possible to replace the curves for the decay and growth of the pulse currents with straight lines as shown in Fig. 
6. The area of the pulses is not altered by this, Then, taking into account the fact that a) at the commutation in- 
stant the magnetic fluxes of two adjacent poles act in opposite directions, and b) at the instant when the commuta- 
ted currents are equal to one another the resultant flux is equal to zero, it is possible to assume that the currents are 
also equal to zero, Thus the “active” area of each current pulse in the winding will be shaped as shown in Fig. 6 
(the shaded potion). Replacing this area by a rectangle of equal size that has a duration equal to ¥% the period of 
the instantaneous value for the "active" current and the electromagnetic torque when the overlap is taken into ac- 
count, we find 
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(22) 
| = la(1—20,). M, — SV 342! 12 sina, +4). 


If we take the basis torque to be that corresponding to 6,, = 0, then 


Y 
Me 4 =): 
w- (1%) 


Since is a function of the parameter 0, it follows that the resulting relationship can be plotted in the coor- 
dinates Mg/M, = f(c); this permits us to obtain the electromagnetic torque at the motor shaft directly for all possi- 
ble values of converter parameters, The resulting expression for ly is valid only for values 6, = 120°, For 4 > 190° 
the active area decreases more rapidly as 9, increases, The magnitude of the torque at the motor shaft is close to 
zero under these conditions, and an investigation of this region is not of interest. 


The angle 6 is proportional to frequency, and as the frequency increases both the active value of the current 
and the electromagnetic torque decrease, 


When (22) is taken into account, the expressions for the rectified current and the electromagnetic torque be- 
come 


ge 





Axl (8x — 36)? sin (a, + = + 49x%e (r + 74) (23) 
49n*wAsl [(Ug— Au)(8x — 39,)}? sin (0 + =) 


[Ar@! (8x — 36) sin (a, + =) 49n%w (r + r4)]* : (24) 
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The limits of variation for @,, the dependence of , on the load, and the stalling torque are found by an ana- 
logous method, It is not difficult to find the converter efficiency: 


Pe 8lx (U4 ~— Au)(8x — 36,) sin (a + 3) (25) 


jee ar re 
$ v,| o (8x ~-36,)?8in (ay + 3) +40n? 4 





From this expression, it is evident that the greater the constant SL/(r+ rg) the greater the efficiency of the 
converter when the remaining parameters stay constant, The correctness of our theoretical assumptions and the re- 
sulting relationships is verified by experiment. 


Figure 7 shows computed and experimental characteristics for the drive. In performing the experiments we 
used a three-stator step motor with a ribbon brake, The basic parameters for the motor and for the electric circuits 
in the drive were the following: AL= 0,24 h; r= 100 ohm; tg = 30 ohm; 1 = 30; 4Su= 15 v; Ug= 215 v. 


For the chosen circuit parameters and the assumed control pulse frequency, the motor was gradually loaded for 
various overlap times of the winding current, and the values of the stalling torque and the currents and voltages in 
the circuit were fixed, 


The author expresses his appreciation to A, A, Bulgakov and Ya, Z. Tsypkin for their help in performing the 
investigations described above. 
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THE CALCULATION OF THERMAL CONDITIONS IN TRANSISTORS 


A. N. Afanas'ev 


(Moscow) 

Translated from Aytomatika i Telemekhanika, Vol. 22, No. 5, 
pp. 641-647, May, 1961 

Original article submitted December 9, 1960 


A method of calculating the static thermal conditions of transistors is presented, The characteristics 
of a cooling plate are developed, Practical recommendations are made for cooling power transistors. 


For the investigation of problems connected with the cooling of power transistors, it is convenient to employ 
the electrical analog method, Essentially this method depends on the principle that the equations depicting thermal 


and electrical processes are identical, and a thermotechnical calculation reduces to a calculation for an electrical 
circuit, 


During the operation of a transistor, part of its power is wasted, being transformed into heat, which warms the 
transistor, If the power is constant in time, and is present for a long period, the temperature of the collector junction 
reaches a static value, In Fig, 1, ais shown the equivalent thermal circuit which reflects the process of power trans- 
mission to the surrounding medium, The thermal resistance ry, ;,,; representing in substance the resistance of the tran- 
sistor, determines the temperature difference 9 -8.45 between the collector junction and the casing. Itis a transistor 
parameter, and is called the internal resistance of the transistor. 


The external resistance of the transistor (casing resistance)Ry,.,, specifies the ability of the transistor to trans- 
mit to the surrounding medium the power wasted within it. 





























as@ 
4, ‘\ 
dr, 
' 94) |49, 
{0 Osur 
P dz 
) 
Fig, 1 Fig, 2 
The total thermal resistance of the transistor is 
Ry = Ty.int+ Ry. cas - (1) 
The maximum power which the transistor can dissipate is equal to 
_ 05 -max— 9sur 
Ss  — (2) 
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6 j.max is the maximum permissible temperature of the collector or junction for a specified collector voltage, and 
© sur is the temperature of the surrounding medium. 


The external thermal resistance of transistors is considerably higher than the internal thermal resistance, This 
prevents them from being operated at their full power, 


Transistor load capability can be increased with the help of a radiator, In Fig. 1, b, there is shown an equivalent 
circuit for a transistor cooled by means of a radiator. In this circuit, Ry, , isthe thermal resistance of the radiator, and 
I'y.¢ is the thermal resistance of the contact which results from the imperfect matching of the transistor and radiator 
surfaces, Usually it amounts to 0,3-0,5°C/w, The value of r'7,,canbe lowered to 0,15-0,2 *C/w if a pad of a soft 
metal— for example, load 0,1-0,15 mm in thickness~ is placed between the transistor and the radiator, 


The external thermal resistance of a transistor cooled by a radiator is equal to 


Ry ac (Ry, + "e.c) 


R;. ; 
i: Recas + Aa + a.c (3) 





As the thermal resistance of the radiator is reduced, the power which the transistor can dissipate is increased, 
For Rp, + tr. 0 it reaches the limiting value 


95. max — 9sur 
= . 4 
him e.toe (4) 
However, practically it is found unsuitable to make radiators with a thermal resistance less than the internal 
thermal resistance of the transistor, 





The maximum temperature of the junction is directly dependent on the collector voltage; this question has 
been considered in detail in [1,2]. 


The transfer of heat from a hot body to the surrounding medium takes place by convection and thermal radia- 
tion [3]. 


Ry rr °"C/w 
12 


Bic 





os © S&S & & 


0 Q4 10 


Fig. 3 Fig. 4. (1) 5 = 1 mm; (2) 6 = 2mm; (3) 6 = 4mm; 
(4) 6 = 2mm; (5) 6 = 4mm, 


The thermal transfer coefficient of convection for a uniformly heated vertical surface with natural cooling is 
determined by the approximate formula 


is 7 40 
% =18)599 (5) 


where A@, is the temperature difference between the hot surface and the surrounding medium in*C, and / is the 
height of the plate in m, 


Experimentally, it has been found that thermal transfer by convection from a horizontal surface turned upwards 
is 30% greater than from a vertical surface. A horizontal surface turned downwards dissipates 30-40% less heat than 
a vertical surface, 
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The thermal transfer coefficient of radiation for a uniformly heated surface is equal to 


a, = 4,60 (1 +38) (6) 


where € is the blackness coefficient of the surface, 9,y = 4(9,+ 9.y,;), and ©, is the surface temperature in °C, 


The simplest radiator for a transistor is a metal plate, Its thermal resistance depends on many factors, which 
are principally related to the geometrical dimensions, the thermal conductivity coefficient of the material, and the 
condition of the surface, 


In Fig. 2, a circular plate is shown to which power is transmitted over a circular region of radius rp. There ts 
also shown in this figure an equivalent thermal circuit consisting of a long line with continuously changing parameters, 


The temperature distribution along the radius of such a plate which is cooled naturally is given by the expres- 
sion (see Appendix): 
Ko (mz) I; (mRo) 4. Ky (mRo) Io (mz) 


= Ko (mro) 14 (mo) + Lo (mre) Ki (mo) 


Ad, + Osur 





(7) 
where Kg mx), Ip(mRg) , etc. are Bessel functions with the imaginary argument m= 742a75X; 2 is the thermal con- 
ductivity coefficient of the material of the plate in w/m* C; a= ay + a,; 6 is the plate thickness in m; r and Rg 
are in m, 
The thermal resistance of the plate is equal to 
|; 1 
rr © 2argdhmyp (8) 
where 


__ 14 (mRo) Ky (mr) — 11 (mre) Ki (mMRo) 
YP = Ki, (mR) To (mre) + Ti (Ro) Ky (mre) * 





The expression given is satisfactory for the calculation of rectangular plates where the ratio of the sides does 
not exceed 1,5-1.8, In this case, the radius of a circusar plate of equivalent area is used for the external radius. 


In order to indicate the inherent properties for a plate, we will develope an analysis of Eqs, (7) and (8) for the 
case when $ (0,;+@y,) = 45°C, 409 = 45°C, 


In Fig. 3, curves of the temperature distribution along the radius of the plate are presented, They indicate 
that with increasing distance from the center, the surface temperature of the plate decreases, The greatest amount 
of decrease occurs when the radius is small. In order to straighten the temperature curve, the plate must be thicker 
in the middle than at the edges. 


The gradual temperature reduction along the radius results in a continuous reduction of the amount of heat 
being transferred to the surrounding medium per unit area as the distance from the center is increased, It can be 
shown that there is an area above which a further increase does not produce anappreciable decrease in the thermal 
resistance, From the curves shown in Fig. 4 for Ry, » = £(S)s = const With bright (curves 1, 2 and 3) and with black ano 
dized (curves 4 and 5) aluminum plates it is seen that there is little value in making a cooling plate with a surface 
area (for both sides) of more than 1000-1300 cm” 


The surface condition of a plate has a marked effect on its thermal resistance. Thus, the thermal resistance of 
a black anodized aluminum plate is 25-40% less than for a bright one having the same dimensions. 


In practical work there are two problems to be solved: 


1) To obtain a given thermal resistance with the minimum plate area; 
2) To obtain a given thermal resistance with the minimum plate volume. 


The first problem is representative of the case when the plate must have minimum over-all dimensions. For a 
fixed area, the thermal resistance of a plate decreases with an increase in thickness, approaching the limiting value 


1 





r=] 





It can be shown that there is a thickness beyond which a further increase does not give an appreciable reduc- 
tion of the thermal resistance of a plate (Fig. 5). This occurs because for 5= 1.5aS/, the temperature distribution 
along the surface of the plate becomes practically uniform. A uniform temperature distribution occurs for aluminum 
plates at thicknesses from 1 to 4mm and a two-sided area of 300 cm’, This is indicated by the convergence of the 
characteristic curves in the first part of the graph shown in Fig. 4, 


The second problem arises when material costs are important, For a given quantity of material it can be shown 
(Fig. 6) that the thermal resistance of a plate is a minimum for mR, © 1, The two-sided surface of such a plate is 
equal to 


gs — 


6:1" (10) 


Let us compare two plates, one of which achieves the minimum over-all dimensions (area), and the other— the 
minimum material cost. The first plate uses the area efficiently, but is not economical with respect to material cost, 
The minimum quantity of material is used in making the second plate, but it utilizes the area inefficiently. Economy 
of material is achieved at the expense of area increase. 


°c/w 


Rp. “C/w 


V2) 
~~ & & Bw OH 





0 / 2 & mm 
Fig. 5 Fig. 6 
The relative rating of plates made of different metals is of some interest, This can be worked out, based on 
the curves in Fig. 7. 


The thermal resistance of a steel plate is considerably in excess of that of copper and aluminum plates of the 
same dimensions, It follows from this that steel is not a suitable material for making cooling plates. 


Copper plates have a somewhat smaller thermal resistance than aluminum, However, from the point of view of 
cost, it is preferable to use aluminum plates, 





R. ,c/ 1/ Reps wi'c 
Tr’ c/w 
J a9 
8 
7 
§ 
5 é&-2mm 
rhs mm 
4 
3 
2 
0 2 


0 cm 


Fig. 7. 1) Steel; 2) bright aluminum; 3) anodized Fig. 8. I) Anodized aluminum; II) bright aluminum 
aluminum; 4) copper. 
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Heat transfer from the surface of a plate depends both on the difference (convection) and on the absolute 
values (thermal radiation) of the temperatures of the surface of the plate and of the surrounding medium. Conse- 
quently, the thermal resistance of a plate is functionally connected to temperature. For an increase in the tempera- 
ture difference between the surface of the plate and the surrounding medium, the thermal resistance of the plate is 
lowered, On the average, it can be reckoned that for a change in the temperature difference of 1°C in the interval 
of AG, from 10 to 45°C, the thermal resistance changes by 0.4-0.% with anodized plates, and by 0,5-0.@% with 
bright ones. 


Equation (8) was given above for determining the thermal resistance of a plate. Because of the presence of 
Bessel functions in this expression, the calculation of cooling plates becomes cumbersome. Simplification of the 
calculation can be accomplished in two ways. The first method consists of constructing from Eq. (8) a curve for 
Rr.r = £(S)§ =const, Of 1/ RT. r = AS)5=const, and finding from ittheplate area needed to obtain a given thermal resist- 
ance. The standard curves for aluminum plates cooling P4 transistors (rp = 1.5 cm) are presented in Figs. 5 and 8, 


The second method of simplifying the plate calculation consists of approximating Eq. (8) with a simpler ex- 
pression, Good results for mRy = 1,1 are given by an approximation of the form 


10* “a 
Rar = a tkV (11) 


In this formula, a is in w/m* °C, S is incm’,A is inw/m°C, 6 is in mm, and k is a coefficient. 
In the calculation of plates for cooling P4 transistors, k is made equal ot 3.5. 


An experimental study of cooling plates showed that the calculated and experimental data agreed well with 
each other, The difference between calculation and experiment did not exceed 20%, which is completely accept- 
able for engineering purposes. 


SUMMARY 


The load capacity of a transistor under static conditions is characterized by the maximum permissible tem- 
perature of the collector junction, the internal and external thermal resistances, and the relation of the maximum 
permissible voltage to junction temperature. A considerable excess of external over internal thermal resistance 
prevents the use of a transistor at full power 


An increase in load capability is obtained for power transistors by means of auxiliary cooling with radiators, 
the simplest of which is a metal plate, From the point of view of cost and cooling ability, an anodized aluminum 
plate undoubtedly should receive preference. 


The thermal resistance of aluminum plates is practically unchanged for an increase in total surface area above 
1000-1300 cm? and a thickness above 3-4 mm. 


The choice of one or another radiator depends on the structural design of the apparatus as a whole, In some 
cases, it is suitable to utilize a plate or the surface of the apparatus chassis; in others this is not possible. 


I wish to thank B, S. Sot-skova for valuable comments on the paper. 


APPENDIX 


Derivation of the Formula for the Calculation of the Thermal Resistance of a 





Circular Metal Plate 





We will demonstrate the feasibility of using the electrical analog method to obtain the needed thermotechni- 
cal relations, 


For the development of the formula we will make the following assumptions, 
1) The plate is homogeneous and has a uniform thickness. 


2) The thickness of the plate is small compared to the outer radius; this permits us to neglect the power 
radiated from its end. 
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3) The thermal transfer coefficient is equal over the entire surface. 


Separated by a distance x from the center of the plate (Fig. 2), the element dx transfers power both along the 


line and to the surrounding medium, The change of the temperature difference between the surface of the plate and 
the surrounding medium, and the change of power along the line are equal to 


dA§ =—q dr, (1) 
dq = — Af dg,. (1) 


Expressions for the thermal resistance of element dry along the line, and its thermal conductivity to the sur- 
rounding medium dg- can be correspondingly presented in the form 





dx 
dr, = 2nx 6A (IT) 
dg,= 4xzxadz . (IV) 
From the joint solution of Eqs. (I)-(IV) we find 
d*A§ 1 dA®§ 
dx? a} dx —A Gm? = 0 (V) 
where m= V¥2a/ir, 


This is Bessel's equation, Its solution will be [4}: 


A§ = Cylo (ma) + C,Ky (mz). (VID 
Here C, and C, are constants of integration and Ip(mx) and Ky (mx) are Bessel functions with imaginary argu- 
ments, 


We find the constants C; and C, from the initial conditions; 


dA 
A6 = AGo for z= Tl, m0 for == Ao. 


After transformations we obtain 


ag = —Kalma) 1s (m Ra) + Ki (mR) To(mz)_ 
1“ Ko (mro) 1; (mRo) + 19 (mro) Ky (mRo) ” 








(VI) 
Substituting in Eq. (I) the value of d4@ found from the last expression, and in addition taking into account 
Eq. (Ill) we get 


Jo= 2xrdaAmpAGo, 
where 


_ 41 (m Ro) Ki (mro) — 11 (mre) Ki (mRo) 
— Zo (mro) Ky (mRo) + 1, (mRo) Ko (mrp) * 





(VII) 
The thermal resistance of the plate is 


ay test 
Tir go  2mrgdAmyp* (IX) 


Analogous relations were obtained in [3,5] by a method commonly used in heat engineering, 
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Original article submitted June 21, 1960 


The paper proposes a time relay network that includes a pulse generator and a pulse counter, 

The pulse generator consists of a contactless magnetic relay (CMR) and a thermister which is 

used as an inertial element. A number of variants of the pulse generator circuit are given for 
the case where a CMR is used with directly and indirectly heated thermistors, 


Time relays based on magnetic amplifiers or thermistors are usually designed on the principle of using the 
transient response time in one or another circuit of the unit for application of the input signal, The basic shortcom- 
ing of such units is the limitation imposed on the time interval by the inertial characteristics of the corresponding 
circuit. 





The generator relay whose network is shown in Fig. 1a is free from this shortcoming. This network includes 
a pulse generator consisting of a contactless magnetic relay (CMR) and a heated thermistor (HT), as well as a pulse 
counter (adder), The generator produces pulses with a specified repetition frequency, The time interval is produced 
by adjusting the counter for a corresponding number of pulses, The error in the time interval after which the relay 
operates is determined basically by the generator instability which can be reduced to allowable limits. 


The introduction of a CMR instead of a conventional electromagnetic relay [1] has made it possible to in- 
crease the service life and the operational reliability of the unit; this applies especially for conditions where vibra- 
tions and appreciable constant accelerations are present. 


1, The Circuit Diagram 





The circuit diagram of the generator time relay is shown in Fig, 1,b. The working body of a heated thermis- 
tor HT-1 is connected into the control winding of the pulse generator, and the heater is connected into the feedback 
winding of the CMR, 


When the switch is open By1,= 0) the maximum current (Fig. 1,c) flows through the heater which is the load 
of the CMR, and the working body of the HT is heated to the temperature 69+ 6 ,, where 69 is the temperature of 
the ambient medium and 9,, is the overheating temperature of the HT corresponding to the maximum current 
through the heater. 


When the switch is closed a current Ij, appears in the control winding; the value of this current is deter- 
mined from the expression 


*See English translation. 





where U is the voltage which supplies the control circuit, R(@9+ Ooy) is the HT resistance at the temperature 9)+ Gov, 


I’ 


U 





7 bas ’ 
Ra -HMow vo + | 


I, is the resistance of the control winding, r, is a multiplier resistance, 


The heating of the working body of the HT begins due to the current flowing through it, The magnitude of 
the HT resistance drops, the current Ic increases, aad for 1, = Ine] (point 2 in Fig. 2,4) the current in the heater drops 
abruptly; as a result the HT begins to cool, its resistance increases, and the current I, decreases to the value | 
(point 3 in Fig, 2,a). For 1, = Ie, the current through the heater again increases stepwise to a maximum and the HT 


is heated. 


The described variations of the current in the circuit of the working body of the HT are periodic in nature, 
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shown in Figs. 2, b and 2, c, 
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Fig, 1 


The graphs showing the variation of the control and load current with time, 1, = f(t) andl = ¢(t), are 


The repetition frequency for the pulses can be controlled by varying the resistance r, in the control circuit 
of the CMR or by varying the bias current in the CMR by means of the adjusting resistor tg, An electromechanical 
pulse counter (a step motor) is used as the added; this device is characterized by relatively small dimensions and a 
simple construction 


It is also possible to use other adders such as counters based on triggers, decatrons, ferrite elements, etc. 


The winding of the pulse counter is connected as the load of the CMR, The adjustment of the counter to a 
specified number of pulses is achieved by switching a special switch S, when the control winding is disconnected 
(Ic = 0), Under these conditions a number of pulses equal to the number of times the switching-on operation is per- 
formed is applied to the counter, 

































When the required time interval elapses, the group of contacts K is closed (Fig. 1, b); these contacts send a 
signa] to the control circuit and disconnect the input circuit of the pulse generator, 
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Fig. 2 


An essential feature of networks which use thermistors is the fact that the thermistor parameters depend on 
temperature, Therefore in order to assure the required operational stability of the pulse generator it is necessary to 
insure temperature stabilization of the unit, This is achieved by introducing an additional heating winding for the 
thermistor HT-1 which is shunted by a second thermistor HT-2, 


2. The Design of the Time Relay 





We shall dwell on the basic propositions which govern the design of the pulse generator. 


As we have already noted, the accuracy with which the time relay operates is determined chiefly by the 
stability of the pulse generator frequency, and therefore the generator should be designed on the basis of the require- 
ments for increased stability. The design of a CMR with increased stability was described in [2] and therefore is not 
described here. 


In a first approximation, the inertia of a CMR can be neglected, and therefore the computation of the dura- 
tion of the pulses reduces to computing the dynamics of the CMP control circuit which includes the thermistor. We 
should keep in mind the fact that the design of a circuit with a HT does not differ in principle from the design of a 
circuitwith a conventional thermistor; this consists of an approximate semi-analytical solution of the equation for 
thermal equilibrium in the circuit. The method for computing the transient response is cited in [3] and therefore is 
omitted here, 


In order to compute the dynamic characteristics of the generator it is necessary to make use of the follow- 
ing original data: the volt-ampere characteristics of the thermistor for 9» and 9+ 96,,., the heating characteristic of 
the thermistor, and the time constant of the thermistor. 


The method for designing the thermal stabilization circuit is cited in [5]. 


3. The Results of an Experimental Investigation 





For the purposes of performing an experimental investigation, we prepared a laboratory mock-up of a time 
relay designed according to the circuit shown in Fig, 1,b. We used mixed feedback in the CMR, and this permitted 
us to reduce the dimensions and weight of the unit. 


The CMR was constructed using toroidal ribbon magnetic circuits made of "79NM" permalloy 0,08 mm 
thick with an outside diameter of 42 mm, and inside diameter of 30mm, and a height of 3mm, 


The specifications for the windings were: w_ = 3500, We = 6000, W.i,, = 570, Wy, = 380. 


The experimentally determined characteristic of the CMR for a bias current of I,;,, = -20ma is shown in 
Fig. 1, c. 
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We used anHTof the type *MMT-4” which had a resistance of 100 kohm at a temperature of +20°C, Its 
temperature characteristic can be described by the expression R = 3,3e ®2 ohm, 





T 
The heater consisted of a coil made of nichrome wire 0,07 mm in diameter 
which was wound on the body of the HT; the wire was first covered with a tery- 
lene film, The resistance of the heater was 1000 ohm, The time constant of the 
heated thermistor was 66 seconds, 


The heating characteristic [4] is shown in Fig. 3. The ac circuit of the 
CMR was supplied from a source with a frequency of 1000 cps and a voltage of 
40 v; the bias circuit and the auxiliary heater were supplied from a 27 v dc 
source, The-rectifier bridge consisted of four diodes of the “D9E" type. 


Figure 4 shows the dependence of t,, tg and T = t, + t, on the control cir- 
cuit supply voltage for I,;,,= -20 ma, ambient temperature of +20°C, and 
a maximum heater current of 25 ma; this corresponds toan overheating of Oo, = 
40°C, 


For u © 80 vy, the curve T = f(u) has a minimum for which t = t,, and 
thus the generator instability which is caused by a variation in the control circuit 
supply vditage becomes minimal, 


This fact should be taken into account in choosing the supply voltage for 
the control circuit. 


As we have noted above, the duration of a period of oscillation depends on 
the bias current. 


As a result of the experiment, we established the fact that for u= 80¥y and 
69 = + 20°C a variation of l,j, from -20 to -25 ma led to a variation of up to 
30% in the duration of a period, This verifies the necessity of stabilizing the bias 
current, 


The effect of temperature is manifested chiefly in a variation of the volt- 
ampere characteristic of the thermistor, A variation in the parameters of the re- 
maining elements of the network due to fluctuations in the temperature of the 
medium causes an insignificant error, Thus, when the temperature varies from 
-40 to +50°C the temperature instability of a period of oscillation for the gene- 
rator did not exceed 1,4, 


We should keep in mind the fact that in these tests the thermistor was lo- 
cated inside a thermostat. 


Thus, when just the HT is thermostated we insure the satisfactory temperature stability of the generator, 


. As we have already noted above, this is achieved by ineans of an auxiliary heating winding which is shunted by 


another thermistor, 





it was approximately 2% ). 


When tests were made under normal temperature conditions the generator showed good stability. Thus, for 
u= 80v, Ibja, = -26 ma and @»9= 20°C the maximum deviation from Top = 90 sec did not exceed 2 seconds (i.e., 


4, Variants of the Pulse Generator Network 





) The described pulse generatcr network is not the only possible one, Figure 5 shows other variants of a con- 


tactless pulse generator network, 


In the networks shown in Figs. 5,a and5,banindirectly heated thermistor is used; in the networks shown in 
Figs.5,c,5,dand5, e adirectly heated thermistor is used, 


The network in Fig.5,b differs from the network in Fig.5a because the thermistor shunts the control winding 





of the CMR and a multiplier resistor R,, is connected in the control circuit, 





stepwise under the same conditions, 


Input 


+o— 





Wherease in the networks shown in 5, aand5 a stepwise variation of the heating current occurs for opera- 
tion or release of the relay, inFigs. 5,c,5,d and 5,e itis the resistance in the control circuit winding which varies 
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However, whereas in Figs, 5,c and5,d semiconductor diodes are used to vary the input resistances, the net- 


workinFig.5,c uses a transistor as a variable resistance. 


The advantage of the networks shownin 5, a and 5, bresides in the fact that there is no galvanic coupling 
between the input and output of the CMR. 
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The operating principle of an automatic optimizer is described and its optimum parameters are 
chosen, The results obtained in testing it with an electronic mock-up of the controlled device 
are given, 


Automatic optimizers consists of devices which are specially constructed for tacking the maximum (or mini- 
mum) of the parameter to be optimized, However, a simple and reliable automatic optimizer can be constructed by 
using well-known mass-produced components, 


An automatic optimizer, which has been developed at the hydropneumatic automation laboratory of the Insti- 
tute of Automation and Remote Control, Academy of Sciences, USSR, is described below. It is to be constructed by 
using individual components of the pneumatic standardized unit system (SUS), which is produced by the Tizpribor 
plant, and of the electropneumatic command device KEP-12U, which is produced by the Kalugapribor plant. The 
possibility of constructing such an automatic optimizer by using mass-produced devices can greatly speed up the 
work on introducing extremal control in the practice of industrial automation, 


The present article describes some of the results obtained in the experimental investigation of an automatic 
optimizer, which operated in combination with an electronic mock-up of the system to be controlled. 
I, Operating Principle of the Automatic Optimizer 

The block diagram of the automatic optimizer (Fig, 1) contains four basic units, 


1, The integrator* 1 with two inputs, In the figure, the inputs are denoted by y. and y, in dependence on 
the sign of the integral that is obtained when the signal is supplied only to one of the inputs, 


2. The summation device 2, which adds the integrator output x,j and the trial oscillations x,. 





3. The floating reading system 3, which consists of an inertial element consisting of the resistor r, the capaci- 
tor c, and the power amplifier A, The output y of the system to be controlled represents the input of the floating 
* Instead of the integrator, a more complex linear device with an optimum weighting function can be provided; this 
optimal device can also be nonlinear (see remark at the end of this article), 








reading system, while its output is the signal y¢), which is obtained by smoothing-out the signal y(t); if the value of 
the output y of the system to be controlled is constant, ym = y. 


4, The control unit 4. Figure 1 shows the contacts k;,..., kg of the control unit as well as the device M 
which controls them, The generator of square trial oscillations x, with the period T, has been organized by means 
of contact kg and a certain constant signal x». The magnitude of the constant signal x9 is equal to the double value 
of the square trial oscillation amplitude, i.e., Xo = 2Xgr, Where Xg_, + Xemax, 


The switching of contacts k;,..., kg alters the scheme of supplying the signals y and y¢, to the integrator out- 
puts yand y. The k;,..., kg contacts are switched simultaneously, These switchings have the same period T,, and 
consequently they are performed in step with the oscillations x,. During a single half-period, the output y of the sys- 
tem to be controlled is supplied, for instance, to the integrator input y, (Fig. 1). At the same time, the output yp 
of the floating reading system is fed to the input y.. During the next half-period, the reverse process takes place, 
etc, In other words, this means that the difference y-y,, during one half-period is integrated with the plus sign, 
while it is integrated with the minus sign, during the next half-period, etc. 


5. The limiter 5 which limits the integrator output value, For a simpler explanation of the operating princi- 
ple of the optimizer we shall assume that the limiter 5 and the floating reading system 3 are absent, i.e., ym = 0; 
we shall also neglect a number of other secondary factors. 


























Consider the operation of an automatic optimizer in conjunction with an inertialess system whose positional 
response is shown in Fig, 2. For the sake of determinancy, let us assume that the system, which consists of the opti- 
mizer and the system to be controlled, is positioned to the left of the extremum (points 1 or 2), 


The operation of contacts k;,..., kg consists in alternately supplying the output y of the system to be con- 
trolled to the two integrator inputs, Thus, the output y of the system to be controlled is integrated with the plus 
sign over one-half of the period T,; during the next half-period, it is integrated with the minus sign, etc. At the 
same time, square oscillations x,, which transfer the system from position 1 to position 2 and vice versa are supplied 
to the input of the system to be controlled (Fig. 2). If the system is to the left of the extremum, the output y of the 
system to be controlled is larger in position 1 than in position 2, Therefore, the integrals over adjacent half-periods 
are correspondingly not equal to each other with respect to absolute magnitude, The sum of these integrals, which 
have opposite signs, yields over one period a positive increment at the optimizer output, which brings the system 
closer to the extremum, 


Similarly, the system approaches the extremum if it was initially positioned to the right instead of to the left 
of the extremum, 


For tracking the maximum instead of the minimum, the phase relationship between the oscillations of con- 
tacts k, and ky,..+.+» Kg should be reversed, 


Let us clarify the purpose of the floating reading system. In the absence of the floating reading system, the 
integrator output x4 performs near-sawtooth oscillations during the scanning and tracking process, which is connected 
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with the fact that the quantity y is alternately integrated with the plus and minus signs, These additional (besides 
x») sawtooth oscillations lower the quality of the optimizer, The greater the output value of the system to be con~ 
trolled (and also the longer the period T and the smaller the integrator's time constant), the greater the oscillation 
amplitude and, consequently, the lower the performance quality, In order to eliminate sawtooth oscillations, it is 
necessary to eliminate the large constant component of the output y of the system to be controlled, For this, it is 
necessary to integrate not the entire quantity y, but only a portion of it, which is measured from a certain variable 
level yp. Thus, the floating reading system filters out the low-frequency oscillation components of the quantity y; 
in this, the x, oscillations are transmitted after they are slightly smoothed out (see sections III and IV). 


Due tothe presence of the floating reading system, 
sawtooth oscillations are absent if the extremum is kept 
stationary, since y-yf) = 0 (Fig. 2). 


In the transient process, the basic component x, of 
the movement toward the extremum (for jump-like varia- 
Ym it} tions of its location) is monotonic, Sawtooth oscillations 
1 can arise only for a short time immediately after a large 
jump in the magnitude of location of the extremum occurs, 


Moreover, due to the presence of the floating reading 
(Zm,Y¥m) system, the effect of unequal time constants in integrating 
| the positive and negative values of the y-yr, difference is 
+— Zz reduced; in the device described here, the above time con- 
| stants differ by 10%, 





Limiter 5 serves for limiting (at the top) the region 
in which the extremal system's output varies, Limiting is 
necessary either for technical reasons or for preventing the 

Zt, (4) "jamming ,"of the system as the control unit stops when the 
extremal position passes beyond the stops, In the pneuma- 
tic optimizer under consideration, the limiting of the range 
from below takes place spontaneously: the pressure xj can 

t not be less than zero, The lower limit of the range is ad- 

Fig. 2, Extremal response of the controlled system; justed to the summation device 2, 

Xm and yr are the extremum coordinates, 





The extremum tracking process is described here only 
in general terms, Actually, the processes in an extremal 
system are somewhat more complicated due to the interation of phenomena when the magnitudes of y, yp, etc. 
change, The particulars of these phenomena can be readily explained by considering individual idealized operating 
regimes of the system in greater detail, Some aspects of the tracking process will be considered in dealing with the 
choice of optimum parameters. 


Il. Description of the Optimizer’s Pneumatic System 





It is assumed that the system to be controlled has a pneumatic output y and a pneumatic input (for instance, 
pneumatic command in an ordinary-type regulator), 


The layout of the pneumatic optimizer is shown in Fig. 3, 
We shall describe the pneumatic equipment used in the optimizer. 


1, The reset control unit 5RB-9B of the standardized unit system is used as the integrator 1, If the 5RB-9B 
unit is used in this way, the pressure in the repeater of the reset control unit serves as its output, In order to bring 
out this pressure, a connecting pipe, u, is provided, while the regular outlet is plugged. In this case, the 5RB-9B 
unit supplies the pressure x,j, Which is proportional to the integral of the difference between the pressures fed to its 
two inputs y. and y,. 


The pressure regulator of the 5RB-9B unit is used as the limiter 5 which limits the maximum pressure at the in- 
tegrator output, For this, the following two channels are covered by means of amembrane that is located under the 
pressure-setting device; the channel through which the feed pressure is supplied to the setting device and the channel 
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through which the setting pressure produced by the setting device are brought out.With these modifications, the con- 
necting pipe “setting” of the 5RB-9B unit serves as the input of the limiter of maximum pressure x,j- (It should be 
noted that limiter 5 can be provided as a separate device.) 


Since the “setting” connecting pipe now serves a different purpose, one of the inputs of the 5RB-9B unit is un- 
available, In order to restore the integrator input, a new connecting pipe, which is marked by y. in Fig. 3, is in- 
stalled in the vertical channel for the setting pressure below the membrane that covers this channel. 
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Fig. 3. c) Pneumatic capacitor; r) pneumatic resistor (choke). 





It is obvious that this modification of the 5RB-9B unit, where two channels are covered and two connecting 
pipes are installed, will hardly affect its reliability. 


The divider D, whichis based on pneumatic resistors (chokes [1]) with attenuation factors that are close to 
unity, is intended for compensating the difference between the effective areas of the y. and y, inputs. 


2. The summation unit 2RB-25A of the SUS system for four quantities to be summed (two with the plus sign 
and two with the minus sign) is used as the summation device 2, 


The integrator output x,; is fed to the two “positive” chambers of the summation device, as a result of which 
the gain with respect to x;j is equal to 2.4, The operating range in which the integrator output x,j varies is thereby 
restricted; this range is of the order of 0,5 atm if x, varies in the 1 atm range. 


The oscillations x, are fed to one of the “negative* chambers of the summation device, The summation 
device's gain with respect to x, is equal to -1.2. 


3, The intermediate relay PR-14 of the SUS system is used as the power amplifier 3 in the floating reading 
system. The power amplifier is necessary because of the large leakages in the pneumatic contact: (pneumatic valves) 
and because the volume c is commensurable with the capacity of the integrator's chambers and piping. 


4. A KEP-12U electropneumatic command device is used as the control unit. Methods for connecting this 
device to the extremal — are shown in Figs. 3 and 4, In this, the outlet to atmosphere in valves k;,..., kg is 


plugged. 
The pressure x» is produced by means of the divider D, (Fig. 3), which is similar to the D, divider, but has an 
attenuation factor of the order of 0,03, 


In constructing the automatic optimizer, the absolute values of the cross-sectional areas of passages in the 
chokes are determined by taking into account the leakages and the capacities of the chambers and piping. 


The KEP-12U device can be replaced by a control unit which is based on pneumatic relays [1]. It contains a 
simple pneumatic oscillator. Since the ratio of half-period durations is unstable, a trigger at the output of which the 
oscillation half-periods are considerably closer to each other is installed beyond the generator, It should be noted 
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that a high oscillation period stability is not required from the generator, The square trial oscillations x, are ob- 
tained from the trigger output through divider D,. The trigger output is controlled by means of four pneumatic 
valves, which are connected in the same way as the KEP -12U valves in Fig. 4. It should be noted that a two-nozzle 
relay can be used instead of each valve pair [1]. 
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The described optimizer can be constructed by using small-size SUS units, In this, the 2RB-25A unit is re- 
placed by a BS-34A unit, and the PR-14 power amplifer is replaced by PR-14M. These units can be used without 
any further modification, The 5RB-9B unit is replaced by a 4RB-32A unit, where the modification consists only in 
installing a connecting pipe in the repeater chamber of the reset control device. Since this unit does not contain a 
setting device, the limiter 5 of maximum pressure x4; is constructed by using a P -280 pressure reducer, where a 
nozzle is installed instead of the valve. This reducer is connected in the same manner as the setting device of the 
5RB-9B unit. Instead of using the 5RB-9B or the 4RB-32A unit, the integrator can be constructed on the basis of the 
layout described in [2]. One of the above control units can be used as the control unit in a small-size optimizer. 


Figure 5 shows an optimizer variant which is based on elements that have been developed at the IAT AN SSSR 
hydropneumatic automation laboratory [1], the mass-production of which has been organized at the Tizpribor plant. 


Il, Adjustment of the Automatic Optimizer 





In adjusting the system, the integrator zero drift is minimized in the first place. As a rule, this adjustment is 
performed only once, when the optimizer is constructed; later, it is necessary only to check and correct the adjust- 
ment if the characteristics of the automatic optimizer equipment have changed. 


In adjusting the optimizer for another system to be controlled of if the extremum tracking conditions have 
changed, the following parameters must be adjusted: the period and the phase of trial oscillations, the oscillation 
amplitude, the floating reading system's time constant, the range in which the extremal system's output pressure x 
changes, and the integrator time constant. 


1. Minimization of the integrator zero drift. The integrator zero drift causes the extremal system to settle in 
an equilibrium position which is at a certain distance from the extremum, This distance is the greater, the greater 
the integrator zero drift and the smaller the coefficient of the parabola that approximates the extremal characteris- 
tics of the system to be controlled. Thus, if a fixed extremum is maintained, a position error arises and, consequent - 
ly, the systems performance deteriorates, In order to reduce this error, it is necessary to minimize the integrator 
zero drift. It should be noted that this error constitutes the main component of the error in maintaining the extre- 
mum in the absence of disturbances and displacements of the extremum position. In the case of intensive x,,,(t) and 
Ym (t) oscillations, the dynamic error and the error due to disturbances can be considerably larger than the position 
error, which will in this case be negligible. 








It is better to minimize the system's zero drift for the optimizer as a whole instead of separately mini mizing 
it at the integrator, The minimization of the optimizer's zero drift can be performed during the short period of 
time when the optimizer is disconnected from the control process without taking it off the system to be controlled, 


The zero drift (its constant component) is minimizedby means of a screw, which determines the nozzle posi- 
tion in the 5RB-9B reset control unit. 


However, when x and y vary in the range from 0 to 1 atm exc, the zero drift changes considerably. Thus, the 
zero drift can have a large average value. Therefore, the indicated changes must be minimized. 


The main reason for linear changes in the zero drift magnitude is the fact that, in making the 5RB-9B units, 
the effective areas of the membranes of inputs y_ and y, as well as of positive and negative feedback are far from 
being equal. Divider D, was provided for eliminating the linear dependence of the drift on the value of y. 


In order to reduce nonlinear changes in the drift with changes in xj, the operating range of changes in x,; was 
restricted; for this, the gain of the summation device 2 with respect to xj was set at 2.4, 


Generally, both methods can be used in compensating the drift with respect to x,; and with respect to y. 


In the local region of changes in x, and y, the zero drift can be compensated more accurately. This region is 
to be found during the optimizer's operation with an actual controlled system. 


2. Fluctuations of the extremum value y,,(t) are manifested as disturbances in tracking the extremum position 
Xm(t). The oscillations of y,,(t) contain high-frequency as well as low-frequency components (the so-called upward 
and downwardextremum “creep"). Theoretical investigations show that, in the sense of stability with respect to low- 
frequency disturbances, the period T, of the square trial oscillations must be as small as possible. The degree to 
which the period can be reduced is limited only by the inertness of the system to be controlled and also by changes 
in the inertness and lag times in the system to be controlled, 


Actually, the smoothing-out of oscillations x, while they are passing through an inertial controlled system 
must not be exaggerated; otherwise, their efficiency would drop, and the error in tracking the extremum would 
begin to increase due to high-frequency disturbances (the high-frequency components of the y,,(t) oscillations). 


In the case of an inertial system to be controlled, the oscillations of the valve k; (Figs. 1 and 4) must lead the 
oscillations of valves ky,..., k, with respect to phase. The lead must be adjusted in such a manner as to compen- 
sate the lag of the trial oscillations x, that is caused by the inertness and lags in the system to be controlled. The 
phase shift between oscillations x, that are supplied to the input of the system to be controlled and the corresponding 
oscillations at the system's output must be as close to zero as possible; the phase shift can also be equal to 2mn or 
wn for a sufficiently small integer n, The reduction in the T, value (and, corresponding, the increase in n) is limited 
by variations in the lag and inertness of the system to be controlled. In taking into account changes in the lag, the 
value of T, must be chosen in such a manner that the oscillation amplitude of the phase shift be considerably less 
than 47. In the opposite case, the efficiency of trial oscillations would drop, and the extremum tracking error due 
to high-frequency disturbances would increase. 


For a large n, it is assumed that the smoothing -out of oscillations xg by the inertness of the system to be con- 
trolled is still not sufficient. 


The adjustment of T, is effected by adjusting the KEP -12U device, 


3. The optimum amplitude of oscillations x2, is determined in such a manner as to provide a compromising 
solution of the contradiction which arises in striving simultaneously to reduce the extremum tracking error caused by 
trial oscillations and the error caused by the integrator zero drift and the high-frequency disturbances at the output 
of the system to be controlled. 


If the disturbances (oscillations of the extremum value y,_(t) ) are not intensive and the position x;,)(t) of the 
extremum changes very slowly, the main cause of the error in maintaining the extremum is the integrator zero drift. 
Although the zero drift was minimized, we did not succed in making it equal to zero throughout the entire range of 
changes in y and x. Then, with a decrease in x) = 2xXgm, the stable equilibrium position of the extremal system, 
which is located to the side of the extremum, becomes increasingly further removed from. the extremum. 


For x» = 0, as a result of the zero drift, the system moves to the stop and the error in maintaining the extre- 
mum will be very large. The error will also be large for exceedingly large values of x9, since the system, whose av- 
erage position is at the extremum, will perform about this position oscillations with an exceedingly large amplitude. 




















It is necessary to find such an optimum value Xp» = 2xgm that the average of the y values in positions 1 and 2 
(Fig.2,a) is at a minimum, 


In the presence of intensive y,(t) oscillations, the error due to the zero drift can be relatively unimportant. 
An exact optimum adjustment of the x value under disturbance conditions is more complicated than the one des- 
cribed above. The optimum % value is chosen in such a manner as to provide a compromising solution of the con- 
tradiction between the error due to disturbances (due to the y,,(t) oscillations) and the error due to the trial oscilla- 
tions, However, the error due to disturbances is difficult to measure. It should be noted that, if any autocorrelation 
function of y,,(t) [and xm(t)}] can be found, the optimum X2Mopt value can be calculated theoretically [3 and 4]. 


4, The time constant of the floating reading system is so chosen as to be equal to one or two periods T, of 
oscillations xp. This is determined by the fact that there is no necessity for a constant output of the floating reading 
system during a single half-period. It is sufficient if the output changes almost linearly. Actually, the area of trial 
oscillations does not change in this case, which means that the efficiency of the trial oscillations also remains un- 
changed. If the time constant of the floating reading system is still further reduced, the nonlinearity (the exponen - 
tial character) of changes in yp during the half-period T,/2 begins to increase considerably. In this, the area of 
trial oscillations x, decreases, and consequently, the efficiency of the trial oscillations also decreases, In connection 
with this, the tracking error due to the y,,(t) disturbances increases, 


However, it is desirable to reduce the time constant, since in this case, the constant component of y is compen- 
sated more quickly and consequently, the sawtooth oscillations are more rapidly eliminated, 


Thus, the magnitude of the floating reading system's time constant is determined by a compromising solution, 
The adjustment of this time constant is effected by varying the values of the pneumatic capacitance c and the 
pneumatic resistance r, 


5. The range of changes in the output x, of the extremal system is adjusted by means of the spring screw of 
the 2RC-25A summation device and the screw in the pressure setting device of the 5RB-9B reset control unit. The 
adjustment should be started with regulating the lower limit of changes in pressure x, since the adjustment of the 
limiter determining the upper limit of changes in x would be otherwise disturbed. 


6. The integrator time constant should be determined in such a manner as to provide a compromising solution 
of a contradiction similar to the contradiction arising in developing optimal servosystems. The system as a whole 
must smooth out sufficiently well the high-frequency components of oscillations of the extremum y,,(t), and it must 
not overly smooth out changes in the position x,,(t) of the extremum that is to be tracked. The low-frequency com- 
ponents of oscillations of the extremum value are filtered out in the process of measuring the systems remoteness 
from the extremum, The indicated low-frequency and high-frequency oscillations of y,,(t) are manifested as dis- 
turbances in tracking the extremum position. The adjustment of the integrator time constant is brought about by ad- 
justing the chokes DD and AP of the 5RB-9B unit. If the DDchoke is completely closed and the AP choke is open, 
the time constant is equal to 20 sec, An exact optimum adjustment of the time constant under the above conditions 
is complicated; the optimum time constant can be calculated theoretically by using the autocorrelation functions 
of the ym(t) and xp(t) quantities [3]. 


In order to improve the quality of extremal control, the gain of the converter at the output of the extremal 
system must be as large as possible, without, however, saturating the converter output during the extremum tracking 
process. In this, the error due to the integrator zero drift decreases, while the error due to oscillations of the extre- 
mum value, i.e., the error due to disturbances, remains the same as before. 


IV. Experimental Investigation 





The pneumatic extremal system was tested in conjunction with an electronic system, which contained a quad- 
trator with inertial elements ahead of it and beyond it, a ball noise generator, and converters for transforming pres- 
sure into electrical voltage and vice versa (EDMU-1 and EPD-32). 


The extremum tracking process under conditions where changes in the extremum positions have the character 
of random switchings is recorded on the oscillograms shown in Fig. 6, 


For low disturbance intensities and small changes in the position of the extremum, the error in maintaining 
the extremum was caused only bythe integrator zero drift. The accuracy in maintaining the extremum for a fixed 
extremum position in the absence of disturbances was equal to 1-1.5% for a 1 kg/em* range of changes in pressure y. 








The accuracy can be increased by a factor of 2 to 3 by more carefully adjusting the divider D, and the screw which 
determines the nozzle position in the 5RB-9B unit, especially if the minimization of the zero drift is performed for 
the local region of changes in the coordinates of the extremum. 





At the present time, industrial tests of the optimizer are in progress in one of the synthetic alcohol plants. 
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Fig. 6. a) Inertialess controlled system without disturbances; b) inertial controlled 
system without disturbances and with 1-sec time constants of the inertial elements; 
c) inertialess controlled system in the presence of disturbances, 1) output y of the 
controlled system; 2) extremum position xp(t); 3) output yp of the floating reading 
system; 4) output x of the automatic optimizer; 5) trial oscillations x,. 


Remark, We shall now determine the position of the described automatic optimizer among other optimizers, 


By comparing the known automatic optimizers [3], it was established that the optimizer described in [3 and 4] 
is apparently more efficient in tracking the extremum in the presence of intensive steady-state random disturbances 
with Gaussian distributions for a controlled system which has at least a little inertia. 


The optimal characteristics of the extremal system described in [3 and 4] were determined in further theoreti- 
cal investigations, On the basis of the analogy with optimal linear servosystems, it can be said that, under the above 
conditions, the optimizer [4] provides the optimum solution of the problem of extremum tracking in a large class of 
linear automatic optimizers, 


It appeared that the indicated optimum characteristics make it possible to simplify the block diagram of the 
extremal system and the construction of its equipment without impairing its performance. The described pneumatic 
realization of the extremal system was developed by using these results. It differs from the system [4] by the com- 
bination of several functions in a single device, 


The integrator represents the continuous element in the described optimizer, In this, the weighting function 
of the entire system is represented by an exponential, which is an optimum function under conditions where the sys- 
tem to be controlled has little inertia, the variations of the extremum position have the character of random switch- 
ings, and the oscillations of the extremum value (disturbance) constitute white noise. Under different conditions, 
for instance, if the extremum position changes linearly over a long period of time, a different weighting function, 
corresponding to a system where a more complicated device is used instead of the integrator, would yield better 
results (see, for instance [4]). 


In the foregoing, we have discussed optimal systems in the class of linear systems. It should be noted that, in 
the class of all possible linear systems, the optimal linear systems are absolutely optimal only in the case of a 
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Gaussian distribution of random steady-state changes in the magnitude and position of the extremum. Otherwise, non- 
linear optimizers can be absolutely optimal [5 and 6]. Thus, if extremum tracking is to be further improved, more 
complex linear and nonlinear unsteady-state devices must be introduced instead of the integrator. It should be noted 
that the advantage of nonlinear systems is strongly pronounced only if the noise-to-useful signal ratio is small (6). 


The optimum adjustment of this extremal system as well as of the system described in [4] can be performed 
automatically by means of another similar extremal system. 
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The present article is concerned with the determination of the noise stability of remote control 
systems with frequency modulation of the carrier and the separation of commands with respect to 
the subcarriers for arbitrary levels of fluctuating noise, The probability of command suppression 
and the probability of spurious command formation are used as the estimation criterion. The in- 
fluence of the structure and the parameters of the receiving device on its noise stability is inves- 
tigated, 


Introduction 





Frequency remote control systems which have a high degree of noise stability are used to an ever increasing 
extent in remote control devices, 


The theoretical calculation of such systems has not yet been treated to a sufficient extent in the literature. 


The present article describes an investigation of the noise stability of remote control systems with frequency 
modulation of the carrier by channel (subcarrier) frequencies, In the receiving device (Fig. 1), the carrier frequency 
is demodulated, while the channel frequencies are divided by means of frequency selectors and are then detected. A 
threshold device and a final device are connected to the detector output in each channel, At each instant of time, 
the carrier frequency is usually modulated by one of the subcarriers, since only a single command is transmitted at a 
time, while an additional command, which is used for protection purposes (for disconnecting the receiver if the trans- 
mitting device breaks down), is transmitted in the intervals between the operating commands, 


The probability of command suppression and the probability of spurious command formation are used as the 
noise stability criterion in command remote control. 


*See English translation, 
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Figure 1 shows a typical block diagram of a receiving device for frequency remote control systems. The com- 
plete circuit of the device comprises band filters 4, $2, and $3, a frequency discriminator, FD, detectors D, and 
limiters, 1 and 2, at the discriminator's input and output. 


In order to find the probability of command suppression in the remote control channel (the channel which trans- 
mits the signal), it is necessary to know separately the magnitudes of noise and of the sinusoidal signal at the channel 
filter $3 output as well as the distribution of the total oscillation envelope. In order to determine the probability of 
spurious command formation in the channel that does not contain the signal, it is necessary to know the noise magni- 
tude and the distribution of its envelope. In calculations, the signal and noise magnitudes at the output of filter $, 
are determined first, and these values are then used for determining the corresponding values at the output of filter $s. 


The following assumptions are made in analysis, 

1. Limiters 1 and 2 are close to idealized limiters. 

2, Filters @ and $3 have a Gaussian characteristic, 

3. Detectors D are linear and inertialess, and they have a transmission factor equal to unity. 


4, The threshold devices at the outputs of the detectors are inertialess with zero time of return to the initial 
state, 
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5. The discriminator output voltage is determined by the expression 
% =kA(o— @2) 


where k is a constant factor, A is the voltage amplitude at the discriminator output, and w, is the frequency of the 
discriminator's “zero” point. 


The above conditions are not very strict, and they are usually fulfilled in actual devices. 


1, Signal and Noise at the Filter 4, Output 





Let us determine the signal powers Qg and the noise power Qn, at the output of filter g, in the presence and 
absence of limiter 1. 


ith Limiter 1 
The signal power at the outputs of the discriminator and the filter $,, whose transmission factor in its transmis- 
sion band is K( a) = 1, will be equal to [1] 


wo (1 — e~*)? 
5 RCT eh (1) 


where D = kA = 4cigk/w, tg is the limiting level of the symmetric limiter 1, z = (Us/U,)}} is the signal power-to- 
noise power ratio at the receiver input (at the output of filter $4), and wy is the frequency deviation amplitude in the 
given channel; wy = 2m fy. 


In Eq. (1), the term (1-e~)* takes into account the suppression of signals by noise. The noise power at the out- 
put of filter $, whose transmission band is Af, = f, -f,, is equal to 
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Cn 2.0 = \W(/) df 
(2) 
where W(f)istheenergy spectrum of noise at the discriminator output. According to [5 and 6], in the presence of 


limiter 1, the noise spectrum for the entire range of z values (0 = z< ~) and the initial spectrum portion, when 
2fmaxMf, = 0.2 can be represented as 
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where Af, is the complete band of output filter 4. 


The first term of expression (3) defines the noise spectrum component that is independent of the frequency, and 


the second term defines the component that is quadratically dependent on frequency. If z = 3, the second term in 
expression (3) can be neglected as a small quantity. 


Expression (3) is valid under the condition that the carrier frequency w, is equal to the center frequency of the 
output filter [$, (w»)] and the frequency of the discriminator’s “zero” point (ws = Wp = wz), i.e., that the modulation of 
the discrimination of the carrier is absent. It was shown in [1] that, in the case of a low limiting level in the presence 
of carrier frequency modulation, the low-frequency components of noise at the discriminator output (2f,,. /Af, = 0.2) 
vary to a slight extent (the high-frequency components of the noise spectrum are subject to the greatest changes). 


Consequently, within the indicated limits, expression (3) can be considered as approximately valid even in the pre- 
sence of modulation. 
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On the basis of (1), (2) and (3), for z = 3, the signal power-to-noise power ratio at the output of filter g, will 
be equal to 





(7) = ee 


(4) 
where 


2ha _ Ah 
"= Th’ Ti = dAf" 


575 











Without Limiter 1 
The signal power at the outputs of the discriminator and filter $2 will be in this case equal to [1] 


Gav= FP Wohagse*[14(5) + (5) 


where I,(x) and 1,(x) are modified zero and first-order Bessel functions of the first kind, 


(5) 


The noise power at the output of filter g, is determined by using Eq. (2), where the spectrum W(f)p for slow 
sinusoidal modulation (fmoqg * Afi/2, 2fmax/Af, = 0.2 and 0 = z< @ is equal to (see Appendix) 


W (fh,=4a°k*(") (th 3 [w (0); mot z (x4) 0.405 Bi On) (6) 
where 


Ww (0); =o =0,405e-°8 for QO<¢2<1.5 


and 





WO); ao= for 1.5<2<oo. 
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On the basis of (1), (2), and (3), for z = 3, the ratio of signal power to noise power at the output of filterggwill 
be given by 





(2 ) ia aste*[1o(5) +4 (5)] wn 
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Figure 2 shows the dependences of the signal and the noise powers at the filter $2 output on the signal-to-noise 
ratio (Z) at the receiver input in the presence (1) and absence (2) of limiter 1, which were calculated by using Eqs. 


(1), (2), and (5); 


yt 


b=K(Y or. d= D* Sa, 





Figure 3 shows the dependences of the signal-to-noise ratios at the 
output of filter $2 on z (1 — with limiter 1; 2 — without limiter 1). 


It must be mentioned that, without limiter 1 and in the presence of 
sinusoidal modulation of the carrier, the intensity of the noise energy spec- 





ae trum at the discriminator output (for z = const) changes in correspondence 
2 with the modulation frequency [1 and 7] (the spectrum intensity increases 
with an increase in the frequency deviation from wy), i.e., in the presence 

771 of the subcarrier, the noise is correlated with the signal. 





2, Signal and Noise at the Outputs of Channel Filters 





In order to calculate the signal and noise powers at the outputs of 

ij channel filters $3, it is generally necessary to consider the passage of noise 
that is correlated with the signal through limiter 2 and filter $3. For such a 
case, exact calculations are very complicated. It can be shown that, in a 
wide range of z, this correlation can be neglected regardless whether limiter 
1 is present or absent, As was shown in [1 and 7], the correlation of noise 
with the signal in the absence of limiter 1 is considerable; however, the noise 
0 ES PAE action at the output of limiter 2 will mainly affect only the shift of the in- 
stant of time when the sinusoidal voltage passes through zero. In correspond- 
ence with this, the calculations can be performed for uncorrelated noise with 
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n = 0. In the case of powerful noise, when the noise voltage that is superimposed on the signal can attain the limit- 
ing level of limiter 2, the correlation between noise and the signal becomes so weak that the difference between noise 
intensities for fy = 0 and fy = max is small. In this case, the calculations can also be performed for uncorrelated noise 
that corresponds to n = 0, 


Until the present time, the signal and noise passage 

through the limiter with a Gaussian distribution has been 

eat a considered in the literature. The noise at the filter 
I 3). a, se) output differs from Gaussian-distribution noise; however, 


EN ss if the band relationship is given by Af, <<Af, it can be 
( ix (232) considered as Gaussian-distribution noise, It should be 
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oa hc noted that the above relationship between the bands is 
> XQ usually valid for industrial remote control systems, Under 
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these assumptions, we shall determine the magnitudes of 
the signal and noise at the outputs of channel filters in the 
(2 ) presence and absence of limiter 2. 


Qo! With Limiter 2 
a0! y U 100 The signal power Qgy and the noise power Q,s (in the 


Bie. 4 frequency band near the input signal first harmonic) at the 
8 output of the idealized symmetrical limiter 2 with the 
characteristic 
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are expressed in terms of their ratio at the output (Q;/Qn), by the following relationships [4} 


a= F(e),4[5 2- (4), ® 
ag = a ePt[$-t-— (St) J +-$Ft[ 4-4. (32) ] + o 
+ art [F-4—(¢) J+] 


where ,F,(8, y , ~x) is a degenerate hypergeometric function, 


2 
Figure 4 shows the dependences of Qg3/ dap and Qns3/ 2a on (Qs/Qn)p. which were plotted by using Eqs. (8) 
and (9) [4]. The figure also shows the broken dashed lines which can be used for approximating these functions, They 
are given by the following simple expressions: 


a8), on 3(s) = cr 
es BH and a) for (z),>1. (11) 
“\Qn/s 


It is obvious from Fig. 4 that the approximation of functions (8) and (9) by the broken curves (10) and (11) 
yields a considerable error in the (Q,/Q,), © 1 region. 


The noise spectrum at the limiter output depends on the limiting level. For a low limiting level, the effective 
width of the output noise spectrum (in the band near the input signal first harmonic) is approximately equal to the 
transmission band of an idealized filter that is connected ahead of the limiter, while the spectrum intensity is uni- 
form [2]. 





In this case, the noise power at the channel filter output in the i-thchannel is given by 


o~ Afs 
Ons a CQ ns Bhs . (12) 


Here, Af, is the energy transmission band of channel filter $3 with the center frequency fj, and Cy is the filter 
transmission factor for the center frequency. 


The signal power at the channel filter output will be obviously equal to 


Qss = CiQs.- (13) 
and the sinusoidal signal amplitude will be 
Ey = V 2Qq- (14) 


Wi er2 
The signal power at the output of channel filter $3 is 


Qse = COQ2 (15) 
where Qg is determined from Eq. (1) or (5), while the noise power is given by 


A 
Om = COW (fi) Mfs = C2Qn2 Fe (16) 
W(fj) is given by Eq. (3) or (6), while the approximate equation is valid for z = 3 in the presence of filter 4. 


3, Probability of Command Suppression 





Let us consider the process of command suppression. 


As can be seen from Eqs. (1) and (5) and Fig. 4, the 
amplitude of the sinusoidal signal E, (14) at the output of 
channel filter $3 and consequently, the signal level at the 





































































































Us output of detector D continuously decrease with an increase 
A ya fs in the noise level at the receiver input. In this, as follows 
. i tae from Eq. (2), the noise power at these points simultaneously 
thr increases, The above phenomenon first leads to partial 
| command suppression and, for a high noise level, to com- 
t plete suppression of the commands, For the sake of simpli- 
Uthr city, we shall consider that the commands transmitted 
: Ae lm through the ith channel under consideration follow each 
& other without pause. In this, the beginnings of commands 
are distributed at random in time. 





Figure 5, a and b, shows the shape of command sig- 
nals at the detector output (the envelope of the total sinus- 
; oidal signal and noise voltage at the output of channel fil- 
Fig. 5 ter $3) for the two limiting cases corresponding to the be- 
ginning of command suppression (a) and to the almost com- 
plete suppression of commands (b) (the rise and decay times 
of the command signal fronts are assumed to be short in comparison with the command duration, which is obviously 
true for Afy >> 1/T%). In Fig. 5, a and b, U,, designates the operating level of the threshold device and E, denotes 
the signal level with respect to which the noise acts; this level is equal to the sinusoidal signal amplitude at the 
linear detector input. It is obvious that the signal level can occupy any of the intermediate positions with respect 
to Uthr, 


If a certain constant level is taken for the signal E, level and the level of Ugh; is assumed to be variable, the 
cases given in Fig. 5, a and b can be represented in a different manner (Fig. 6). 














It is obvious from Fig. 5, a and b, and Fig. 6 that the suppression of command signals will occur either due to 
negative overshoots of the detector output voltage that intersect the Upp, level (a) or due to the intervals between the 
positive overshoots of this voltage that exceed Uh, (b). Suppression will occur in the case where the command signal 
voltage is lower than the operating level U,,,, over the entire period 7, during which the command is mg (Fig. 
6). If the duration % of overshoots at the level U,,, (a) is less than tes suppression will not occur. If % =%,, only 
those commands whose beginning lies in the Ay = %- % interval will be suppressed. If we consider that the com- 
mand signals arrive at random, we find that the probability of signal suppression by a single overshoot with duration 


%) is equal to 
Rus Une) = 2S 


where T is the time interval under consideration, which is of much greater magnitude than the fluctuation process 
correlation time (T>> 7). 


By considering all overshoots % =% or the intervals 















































bs 1% =T between them and by averaging over the time T, we 
-$——— |} ——_+ find that the probability of command suppression at the 
foo hi —/}\ devek. Uthr level for the case a or b will be given by 
si \y ¢ 
4 * thr ™>%) —t 
Ea vs \ L Pte Ru(Gne yee = ne (Unr » T>%) 
a (17) 
iT) 4 where T (T =1%) is the average duration of overshoots (for 
Uthr < E,) oF of the intervals between them (for Ugh; < E) 
t atthe Uhr level of all F = (T = T, for the case a and 
= % for the case b), and a, (Urn. T =T,) is the average 
Fig. 6 ond. of downward overshoots of all T St, during the 


time T. 


The values of T(T =%,} and ny(Uppre 7= %) are determined by the following expression: 


{ w (t) dt 
tc 
T(t >) = Sw (x) dt 
*e 
ny ( (Chr - tT >t) = AU, )TP(t>e) =nUm)T (w(nyar (18) 


t 
c 


where w(T ) is the distribution density of the durations of overshoots or intervals between them at the Uj), level, 
n(Urhr) is the average number of downward overshoots at the Urhy level per unit time, and P (T = 1,) is the probability 
that the overshoots (or the intervals between them) will have a duration T = 7, at the Up, level. 


By substituting (18) in (17), we obtain the expression for the probability of command suppression. 
After normalization, we obtain the equation for the probability of suppression at the level a: 
Py = na) Ir (a. Tat) —td PG, Tt.) = 


= n(a) [ 5 tw (t)dt—t, { w(x) dr] a%) 
where a is the relative limiting level *e *c 
ae Ghr — 
ae (20) 


mMa) is the average number of overshoots of the envelope of the total signal and noise voltage at the a level per unit 
time, and w(T ) = w(a, 7) is the distribution density of the durations of overshoots (or the intervals between them) at 
the level a, 
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The distribution of the linear detector output voltage in the channel with the signal (the envelope of the sinus- 
oidal signal and noise) is expressed by the generalized Rayleigh law. For large signal-to-noise ratios(E,/Un, > 4) this 
distribution approaches the Gaussian distribution. 


For the Gaussian as well as for the Rayleigh distribution, an exact analytical expression for the distribution of 
durations w(T ) has not yet been found, The experimental data provided in a number of papers [8 and 9] indicate 
that the distribution of the durations of overshoots and intervals between them for a Gaussian and a Rayleigh process 
can be coarsely approximated by the Rayleigh distribution for all limiting levels. 


If we use such approximation, i.e., if we assume that 





expression (19) will be given by 
Peake): - 
P =n(a)\e ***’ dr= 2¢n (a)(1 — F (1,25-£)| 

su T 

: (21) 
c 

where T=T(a) isthe average duration % of overshoots (for Uyp,,<Es) or the average duration % of the intervals 
between them (for Uppy > E,) of all 720 at the level a and F(x) is the Laplace function 


x 
chee (e-ma 
F(t) =F Se dt 
It is obvious from Eq. (21) that, for determining the probability that the command (with the duration %) will 

be suppressed, it is necessary to know the relative limiting level a (20), the average duration T of overshoots (or in- 
tervals between them) at this level, and their average number n(a) per unit time. As follows from (1), (5), (8), and 
(13), the limiting level a depends on the signal-to-noise ratio (z) at the input of the receiving device and on its 
structure and parameters, 


Assuming that the detector output voltage in the control channel has a Gaussian distribution (Ez /Un, > 4), we 
have (3): 





n (a) ag bf,e—*”, (22) 
om ~ i RM SES® .. F (a) 
n (a) ofc"? (23) 


where 6 f, is the average width of the fluctuation spectrum; for a Gaussian frequency response of the channel filters 
#3 with the energy transmission band Afg, it is equal to 6 fy ¥0.4 Afy. 


By substituting (22) and (23) in (21), we obtain an expression for the probability of commandsuppression for an 
arbitrary level a. 


In the general case, this expression will have a complex form. Simpler expressions for P,, are obtained for 
particular a values, 


Beginning of Suppression (Fig 5. a). 
In this case, the value of z is negative, and the average overshoot duration is % «7. If | al> 2, we have: 
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Nearly Complete Suppression (Fig. 5, b). 
In this case, the value of a is positive, and the average duration of the intervals between overshoots is Tj >» T,. 
If a> 2, we have: 
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Then 


Puy = 1 — ft, e-2%, (25) 


By expressing the values of E?(14) and Up, in terms of the signal power Qg, and the noise power Q), at the out- 
puts of channel filters, we obtain the equation for the magnitude of a: 





_ OB, max — Fs [ s4.max ou | 
is Uns =¥el6 One Qn (26) 


where o& = Uthr/Es,n4, is the relative operating level of the threshold device, while Es,,,, = Es for z-> =, 


For the case where limiter 2 is present, Qs, and Qn, are found from Eqs. (12) and (13) by using the signal power 
Qyg and the noise power Qns at the output of limiter 2, which are determined from the curve ziven in Fig. 4 or from 
the approximate analytical expressions (10) and (11). 


In the latter case, the expressions for a are given by 
a= VIV $F a—0.88)/ (),] for (G).<! 
a= 2 Sh(a—1) V ®), és (@),>1 (21) 


In the absence of limiter 2, the values of Qy and Q,, are determined by using Eqs. (15) and (16). The signal 
power Qg at the output of filter , and the noise energy spectrum W(fj) at the discriminator output can be found by 
using Eqs, (1) and (3) if the limiter 1 is present or by using Eqs. (5) and (6) in the case where the limiter is absent. 


4, Probability of Spurious Operation 





The formation of spurious commands in the frequency system under consideration can take place due to the 
noise voltage that acts at the outputs of the channel filters in the absence of the signal corresponding to the given 
filter (the case of channels without the signal), The noise voltage level depends on the signal-to-noise ratio at the 
receiver input and is determined by Eq (12) or (16). Let us determine the probability that a command will appear in 
the kth channel during the period of time { in which a command is transmitted through the ith channel. It can be 
shown that this probability can be calculated on the basis of the equations derived for the suppression probability. The 
curve giveninFig.5,c shows the detector output voltage in a channel without the signal. We shall assume that the 
moments of arrival of command signals in the ith channel are distributed at random. Spurious operation will occur 
if the fluctuating voltage exceeds the U,,, level at least at a single point during the time interval % that is equal 
to the command duration. There will be no spurious operation if the time interval ™% falls within the interval be- 
tween overshoots, i.e., if the fluctuating voltage does not exceed Up, at any point during the time %. It can be 
readily seen that the pattern of the absence of spurious operation is similar to the above considered pattern of the 
presence of suppression, where the probability of the fluctuating voltage being everywhere lower than U,,,, during 
the time interval % was determined, 


The probability of spurious operation at the level ay) = Upp, / Up, can be determined by taking into account 
expressions (19) and (21), since 


is i one 1 —n (a) [%; (do, Ty >) —T,] P(e. >t) = 
- tT 
= 1 — 2%, n (a) E —F (1.25 2] (28) 
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where T, is the duration of the interval between positive (with respect to the zero level) fluctuating voltage over- 
shoots at the a) level, n(a9) is the average number of voltage overshoots at the a) level per unit time, 1% (a, 7] =%) 
is the average interval duration at the a level of all% 2 %, 7% is the average interval duration at the a) level of 
all % 20, and P(a», % =7,) is the probability that the intervals at the a) level have the duration % =T,. 


The detector output voltage in the channel without the signal (the noise envelope) has a Rayleigh: distribution, 
In this case, the values of n(a,) and i, which determine Psp are equal to [3] 


n (as) = V2xd/,a,e°8?. |, = P<) ite 
1 n (a) ‘se ate 
V 22 Sfsaye 





We shall determine the expression for Rp at the beginning of spurious operations. In this case > To, and 


F (x)= 5+ + (e= 1.254). 


Then 
P,p= 1 —n (ag) (t; — te) 


If we take into account that the average duration % of overshoots and the average duration % of the intervals 
between them are related by the expression 


T) + % == n (ao) 


we obtain the following equation for the beginning of spurious operations; 


R= N (ao) (t) + Te). (29) 
If ag-+0, then 7 « 7% and as follows from (28), Psy 1. 


It should be noted that the equations obtained here for the probability Psy of command suppression and the pro- 
bability Psp of spurious operation can also be used in the case of signals whose shape at the detector output differs 
from a rectangular shape. In this case, the signal should be approximated by a trapezium, and the command dura- 
tion % should be understood as the duration of the upper flat portion of the trapezium. The proposed method for 
analyzing the noise stability of frequency remote control systems can obviously be used for determining the noise 
stability of other systems where the effect of fluctuating noise on the pulse signal reception is investigated. 


5. Effect of the Receiving Device'’s Structure and Parameters on its Noise Stability 





Let us consider the effect of limiters at the discriminator input (1) and output (2) on the frequency receiver 
noise stability, We shall make a comparison between three receiving device structures; 1) limiter 1 is connected 
and limiter 2 is absent; IL) limiter 2 is connected and limiter 1 is absent; IIl) both limiters 1 and 2 are connected. 
Figure 7 shows the dependence of the signal power and the noise power at the outputs of channel filters $3 on the 
signal-to-noise ratio (z) at the receiver input in the case of equality of the output signal power for z + #, The cal- 
culations were performed by using the formulas derived above and the curves for 

2 
ip pin =5 and 20 n = =74 = 1 (¢= <2). 
It is obvious from Fig. 7 that, for z < 3 and for the receiving device structure II, the signal power has a larger, and 
the noise power a smaller, value than is the case for structure III, i.e., structure II secures a greater noise stability 
than structure III does, In comparison with structure II, structure I provides a lesser probability of spurious operation, 
but a considerably greater probability of command suppression,as is shown by calculations, 


From the equations for Psy and Psp, it is obvious that, by varying the relative operating level a of the threshold 
device, these probabilities can be redistributed. Since it is difficult to compare systems with different Ry and Psp, it 
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is necessary first to compare the magnitudes of one of the probabilities, for instance Py: It can be shown that, in this 
case, the noise stability of the systems is determined by the signal-to-noise ratio at the output of the channel filter $,. 
Figure 8 shows the dependence of the signal-to-noise ratio &;/Up,¢y on the z value for the three receiving device 
structures, It is obvious from Fig. 8 that, for z< 3, structure II yields the best signal-to-noise ratio at the output, and 
consequently, the greatest noise stability of all the three structures, Thus, for weak signals (z< 3), which case is of 
the greatest practical interest for remote control systems, it is advisable to use a receiving device with limiter 2 at 
the discriminator output and without limiter 1 ahead of the discriminator, It should be noted that the conclusion that 
a receiver having only the limiter 2 has the greatest noise stability is valid only for those remote control systems 
whose operation is highly reliable for a signal-to-noise ratio z< 3 at the input. As a rule, telemetering systems re- 
quire large signal-to-noise ratios for reliable transformation of information, and therefore, the above conclusion may 
not be valid for them. 


Besides selecting the device's structure, the correct choice of its parameters is also an important task. 
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Fig. 7 


The receiver under consideration has three filters with the transmission bands Af,, Af,, Afy. Let us consider 
the influence which the band magnitudes of these filters exert on noise stability. It can be readily shown that the 
noise stability can always be improved by narrowing the band of output filter $3; however, this narrowing is limited 
by the system's speed, As can be seen from the complete equation for (Qs/Qn), (which is valid for 0 =z < ), the 
system's noise stability can be improved by restricting the input filter @, band to the optimum value Afjopr. The 
value of Afjopr that corresponds to the maximum of (Qs /Qn), can be determined by using the complete equation for 
(Q3/Qn)a. For 4 <5 the value obtained for Afjopt will be approximate, since the equation for the spectrum W(f) 
yields an accurate value only for », =5. 


Let us consider the influence of the Af, band in the presence of limiter 2 (structure II), 


It is obvious from Fig. 7 that the narrowing of the Af, band (the increase in y,) leads, for the assigned relative 
operating level a of the threshold device, to an increase in the noise power Qn, and consequently to a greater proba- 
bility of spurious operation. However, in this, the signal power Qg, increases and, as calculations, the probability P,, 
of suppression decreases, Thus, for a = const and in the presence of limiter 2, the Psy and Pep values can be redistri- 
buted by varying the band Af,. By increasing the transmission band Af,, the magnitude of Pg, can be reduced to very 
small values regardless of the presence or absence of the carrier. Calculations show that, for intance, for Upp; = 0.5 
Esimax* Sfs = 20 cps, and % = 0.1 sec, it is necessary that the condition Af, >50Af, be satisfied in order to secure 
Pap < 0.001%. 


As was mentioned above, the other method of redistributing Ps, and P,,. consists of changing the operating level 
of the threshold device, By taking into account both of the above-described possibilities of redistributing P,,, and 
tp it can be shown that a reduction of the band Af, would result in greater noise stability, This is confirmed by 


g. 8, from which it is clear that an increase in y, leads to an increase in the signal-to-noise ratio at the output of 
filter $s. 





It should also be mentioned that, moreover, the introduction of limiter 2, leads to the equalization of the noise 
stabilities of all channels in the device. In the absence of limiter 2, the channel with the highest subcarrier has the 
worst noise stability (for z> 1). 
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Fig. 9 


6, Experimental results 





The experimental investigations of the noise stability of the frequency system (Fig. 1) were performed on a ra- 
dio remote control system with a portable transmitting device, which has been developed at IAT AN SSSR. 











Figure 9 shows the dependences of P,, and Ps, on (Us/Un) = 4% (for T, = 0.1 sec), of Peay on %, (for (Us/U,), =0.6 
and & = 0.7), and of Psy and Py, on @ (for U,/U, = 0.6 and % = 0.1 sec), which were obtained by experiment (1) and 
by calculations (2) for the case where limiter 1 is absent and limiter 2 is present (the transmission bands of the input 
filter, the second filter, and the channel filter were Af, = 20 kc, Af, = 1 kc, and Af, = 20 cps, respectively), The Rp 
vs (Us/Upn), Curve was plotted for a = 0.35, since, for « = 0.7, the value Of Psp is almost equal to zero for any noise 
level at the input. 


Figure 9 shows that the experimental and theoretical data are in fairly good agreement. 


SUMMARY 


1, Spectral methods which were formerly used for frequency telemetering systems [5 and 6}, make it possible 
to derive theoretical equations for estimating the noise stability of remote control systems in the presence of arbitrary- 
level fluctuating noise, 


2, The theoretical and experimental results indicate that frequency remote control systems have high noise sta- 
bility and that they can operate with rather low signal-to-noise ratios at the receiver input (up to(Us/Up), = 0.6 - 0.7, 
for y; = Afy/2Af, = 10, y2 = Afy/2Af, = 50 and % = 0.1 sec), 


3. In remote control systems which provide highly reliable operation for a signal-to-noise ratio (Us/Up), at the 
receiver input that is less than 43, it is advisable to use a receiving device with a limiter at the discriminator output 
and without a limiter ahead of the discriminator. 


APPENDIX 


Noise Energy Spectrum at the Discriminator Output in the Absence of Limiter 1 in the 





Case of Carrier Sinusoidal Modulation 





In the absence of carrier modulation, when its frequency uw is equal to the input filter center frequency w) and 
the discriminator "zero" frequency wz, the noise energy spectrum at the discriminator output at the zero frequency 
f = 0 [1] can be represented with a good approximation by the relationship 

Af . 
W (0) ag = 40"! (3) (U, 2W (Oj mo 


where 
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If the signal frequency deviates from wp), the spectral intensity increases approximately following a quadratic 
dependence and attaining, for fy = Af,/, an intensity that corresponds to the action of noise alone without the signal: 
W (0)8(W(0)8 = W(O)P, = 0 for z = 0 [1 and 7]. 


For an arbitrary value of {,, the maximum energy spectrum value W (0)P max will be equal to 


. . 2f . . . 
WN, nay = Op ot (<4) [W O)n—W (0); pol 
In the case of slow sinusoidal modulation (fmoq* Af,/ 2), the noise energy spectrum for f = 0 can be deter- 


mined by using the equation 


nari A “ i ( 2fa\* ° 
WO = dnt *f) OW Ofame + + (we) [0.405 —W (O);gue]} 


For sinusoidal modulation of the carrier frequency and if limiter 1 is not provided ahead of the discriminator, 
the noise spectrum intensity W(f) is approximately uniform for all low-frequency noise components (2fm,/Af, = 0.2) 
[1], on account of which one can consider that W(f) © W(0) within the indicated limits, 
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